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Abstract  of  Dissertation 

Wall  Effects  Observed  in  Tissue  Equivalent  Proportional  Counters 


Tissue-equivalent  proportional  counters  (TEPC)  have  been  used  to  measure  energy 
deposition  in  simulated  volumes  of  tissue  ranging  in  diameter  from  0. 1  to  10  pim.  There 
has  been  some  concern  that  the  wall  used  to  define  the  volume  of  interest  could  influence 
energy  deposition  within  the  sensitive  volume.  These  “wall  effects”  occur  because  the  wall 
has  a  density  significantly  greater  than  the  cavity  gas.  Energy  deposition  measurements 
were  made  for  1  GeV/nucleon  Fe  ions  in  a  TEPC  simulating  1, 2,  and  3  pim  diameter 
spheres  of  tissue.  The  detector  was  nested  within  a  particle  spectrometer  that  provided 
identification  and  flight  path  of  individual  particles.  Energy  deposition  was  studied  as  a 
function  of  pathlength  through  the  detector.  Approximately  25  %  of  the  energy  transfer 
along  trajectories  through  the  center  of  the  detector  escapes  the  sensitive  volume.  The 
response  of  the  detector,  for  trajectories  through  the  detector,  is  always  larger  than 
calculations  for  energy  loss  in  a  homogenous  medium.  This  enhancement  is  greatest  for 
trajectories  near  the  cavity/wall  interface.  An  integration  of  the  response  indicates  that 
charged  particle  equilibrium  is  essentially  achieved  for  a  wall  thickness  of  2.54  mm. 
However,  estimates  of  LET  and  quality  factor  are  influenced  by  these  wall  effects. 

Detector  response  for  fragment  particles  through  the  detector  was  found  to  scale  closely 
to  Z2  for  ions  of  18  <  Z  <  26,  as  expected  by  theory  for  ions  of  high  velocity. 

The  mean  detector  response  to  particle  trajectories  where  a  charge-changing  nuclear 
interaction  occurred  in  the  detector  wall  was  lower  than  that  of  the  response  to  particle 
trajectories  not  involved  in  charge-changing  nuclear  interactions. 

Steven  Edward  Rademacher 
Radiological  Health  Sciences 
Colorado  State  University 
Fort  Collins,  CO  80523 
Fall  1997 
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Chapter  1:  Introduction 

Zirkle  et  al.  ( 1952)  first  introduced  the  concept  of  linear  energy  transfer  (LET)  to 
explain  specific  variations  in  biological  response  from  different  radiation  types.  In  part,  the 
physical  basis  of  relative  biological  effectiveness  (RBE)  is  believed  to  be  due  to  differences 
in  the  spatial  distribution  of  ionization  along  the  track  of  charged  particles  (Katz  et  al. 

1972).  Rossi  and  Rosenzweig  (1955)  developed  the  first  device  to  measure  energy 
deposition  in  simulated  volumes  of  tissue  with  diameters  on  the  order  of  1  micron  (/« n). 
These  devices  are  referred  to  as  tissue-equivalent  proportional  counters  (TEPC). 
Microdosimetric  techniques  have  been  used  in  radiation  protection,  radiobiology,  medical 
radiotherapy,  and  microelectronics .  The  technique  has  been  used  by  the  National 
Aeronautics  and  Space  Administration  (NASA)  for  characterizing  radiation  exposure  to 
personnel  on  manned  space  flights  (Badhwar  et  al  1992)  where  the  radiation  environment 

is  a  complex  combination  of  electrons,  y-rays,  neutrons,  protons,  helium  and  heavier  ions. 

Of  these,  particular  attention  has  been  given  to  the  high-energy  particles  with  Z  greater  than 
2  (HZE)  (NCRP  1989). 

The  concept  of  quality  factor,  Q,  has  been  introduced  to  modify  the  absorbed  dose  and 
define  an  equivalent  dose  that  is  related  to  risk  associated  with  the  incident  radiation.  The 
International  Commission  on  Radiological  Protection  (ICRP)  defines  Q  as  a  function  of  the 
unrestricted  linear  energy  transfer,  L  (ICRP  1977).  The  International  Commission  on 

Radiation  Units  and  Measurements  (ICRU)  has  defined  Q  in  terms  of  lineal  energy 

where  e  is  the  energy  deposited  to  the  matter  in  a  volume  by  a  single  energy  deposition- 

event  and  l  is  the  mean  chord  length  in  that  volume  (ICRU  1983).  Because  lineal  energy 
is  directly  related  to  discrete  energy  deposition,  from  a  biological  standpoint  it  may  be  more 


meaningful  than  L.  One  objective  of  this  study  is  to  determine  how  energy  deposition  in  a 
TEPC  by  HZE  particles  with  known  L  is  influenced  by  the  composition  of  the  detector. 

A  conventional  TEPC  consists  of  a  solid  wall  surrounding  a  gas-filled  cavity  (sensitive 
volume).  Wall  effects  caused  by  the  difference  in  density  between  the  two  media  were  first 
described  by  Failla  and  Failla  (1960)  and  later  by  Rossi  (1968).  Oldenburg  and  Booz 
(1970)  classified  these  effect  into  different  categories.  Kellerer  (1971a  and  b)  provided  a 

detailed  theoretical  analysis  for  heavy  charged  particles,  where  the  V-effect  and  6-ray  effect 

are  most  important.  The  V-effect  is  the  distortion  created  by  the  exaggeration  in  number  of 
nuclear  collisions  that  occur  in  the  wall  as  compared  to  the  homogenous  case  (when  the 
cavity  has  the  same  density  as  the  wall). 

The  6-ray  effect  is  divided  into  one  of  two  cases:  direct  events  for  when  the  primary 

charged  particle  trajectory  crosses  the  cavity  and  indirect  events  for  when  the  primary 
trajectory  passes  outside  of  the  cavity.  For  direct  events,  energy  deposition  in  the  TEPC  is 

enhanced  by  the  wall  because  energetic  6-rays  created  in  the  wall  can  deposit  energy  in  the 

cavity.  For  indirect  events  that  penetrate  the  wall,  event  size  is  enhanced  provided  the 
trajectory  of  the  primary  is  sufficiently  close  to  the  cavity/wall  interface.  For  trajectories 

distant  from  the  interface,  6-rays  are  stopped  in  the  high  density  wall  before  reaching  the 
cavity.  Indirect  events  from  particles  that  pass  outside  of  the  wall  will  have  energy 
deposition  that  is  suppressed  due  to  absorption  of  6-rays  in  the  wall  of  the  detector.  If  the 

wall  is  thicker  than  the  penetration  range  of  the  maximum  energy  6-ray,  deposition  of 

energy  in  the  cavity  by  electrons  will  not  occur. 

Experimental  measurements  have  been  performed  to  characterize  wall  effects  of  HZE 
particles.  Rodgers  et  al.  (1973)  compared  the  response  of  walled  and  nearly  wall -less 
TEPCs  in  the  field  of  279  MeV/nucleon  nitrogen  ions,  B(v/c)  =  0.64.  Data  with  minimum 
absorber  material  upstream  of  the  detectors  indicated  that  the  dose  mean  lineal  energy  of  the 
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walled  detector  was  about  30  %  higher  than  that  of  the  nearly  wall-less  one.  The  difference 
in  dose  mean  lineal  energy  between  the  two  detectors  decreased  when  the  velocity  of  the 
beam  was  reduced  by  slowing  the  incident  particles  in  absorber  material.  This  observation 
would  be  expected  because  for  secondary  electrons  from  low  velocity  ion  interactions,  the 
site  of  energy  transfer  and  deposition  are  in  close  spatial  proximity.  However,  for  high 
velocity  incident  ions,  energetic  secondary  electrons  can  deposit  energy  at  large  distances 
from  the  location  of  initial  energy  transfer. 

Dicello  et  al.  (1991)  found  similar  results  in  the  comparison  of  walled  and  nearly  wall¬ 
less  detectors  in  beams  of  Fe,  Ar,  and  Ne  with  respective  energies  of  535, 570,  and  557 

MeV/nucleon  (p:  0.772  -  0.784). 

The  distributions  of  energy  deposition  collected  with  walled  and  wall-less  detectors  will 
differ  for  high  velocity  ions.  However,  because  of  the  Fano  Theorem  (1954),  the  mean 
absorbed  dose  should  remain  unchanged,  provided  the  frequency  of  nuclear  interactions  in 
the  wall  are  negligible.  This  principle  must  be  satisfied  for  conventional  dosimetry  using 
ionization  chambers. 

Although  the  mean  dose  measured  by  the  two  types  of  TEPC  will  remain  the  same, 
estimates  of  Q  based  on  frequency  distributions  of  y  can  be  significantly  different  because 
of  wall  effects,  differences  among  the  recommendations  of  scientific  commissions,  and  the 
assumptions  that  must  be  made  in  the  calculations. 

NASA  has  assumed  that  y  is  equivalent  to  L  for  energy  deposition  events  measured 
during  manned  space  flights.  Because  energy  deposition  is  influenced  by  the  pathlength  of 
the  ion  through  the  detector,  clearly  this  is  not  true  on  an  event  by  event  basis.  The 

implications  of  this  assumption  for  determining  Q  for  microdosimetric  single  event  spectra 
has  not  been  studied. 

In  this  study,  microdosimetric  measurements  were  collected  with  two  spherical  walled 
TEPCs  using  1  GeV/nucleon  56Fe  ions.  The  wall  thicknesses  were  1.27  and  2.54  mm. 
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The  TEPCs  were  placed  within  a  particle  spectrometer  system  that  provided  identity  and 
trajectory  through  the  TEPC  for  each  incident  particle.  TEPC  response  was  characterized 
for  Fe  particles  through  the  TEPC,  fragments  of  Fe  through  the  TEPC,  and  for  events 
where  the  Fe  was  involved  in  a  charge-changing  nuclear  interaction  in  the  detector  wall. 

This  study  of  the  wall  effect  for  heavy  ions  is  original  and  different  in  comparison  to 
earlier  work  in  many  respects.  This  study  identified  the  particle  trajectory  through  the 
TEPC  for  each  event.  This  work  was  for  ions  of  energy  1  GeV/nucleon,  while  the  highest 
energy  ions  used  previously  were  600  MeV/nucleon.  Because  particle  spectrometer 
systems  were  not  used  in  earlier  work,  the  beams  were  comprised  of  primary  particles  as 
well  as  fragmentation  products  and  assumptions  were  made  concerning  the  spatial 

uniformity  of  the  beam.  Characterization  of  the  6-ray  effect  and  calculations  of  Q  were  not 

made.  In  previous  work,  very  high  energy  deposition  events  in  the  spectra  collected  with 
walled  TEPCs  were  absent  in  the  spectra  of  walless  TEPCs.  Authors  could  only  speculate 
as  to  the  cause.  Because  a  particle  spectrometer  system  was  used  in  this  work,  the 
following  tasks  could  be  accomplished: 

1.  characterization  of  the  6-ray  effect  to  a  monoenergetic  beam  of  Fe  particles, 

2.  calculations  of  Q  for  Fe  particles  (based  on  a  known  value  of  L  for  each  event  and 
the  relationships  between  Q  and  L  of  ICRP),  and 

3.  more  concrete  evidence  as  to  mechanism  for  production  of  very  high  energy 
deposition  events. 

Also,  in  addition  to  analysis  of  detector  response  to  S6Fe  ions,  preliminary  evaluations  are 
provided  for  detector  response  to: 

4.  ions  of  lower  Z  than  Fe  (with  the  same  velocity  as  Fe  particles)  and 

5.  incident  Fe  particles  where  the  particle  is  involved  in  a  charge-reducing  nuclear 
interaction  in  the  wall  of  die  detector. 
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Both  of  these  evaluations  are  considered  preliminary  because  of  the  limited  number  of 
events  available  for  analysis. 
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Chapter  2:  Literature  Search  and  Background 

Ionizing  Radiation  and  Relevance  to  Biological  Systems 

Ionizing  radiation  transfers  energy  to  materials  through  ionizations  and  excitations  of 
electrons,  and  nuclear  interactions.  Biological  effects  from  exposure  to  ionizing  radiation 
are  determined  in  part  by  several  characteristics  of  the  radiation  field.  Among  the 
characteristics,  absorbed  dose  is  perhaps  the  most  important.  Absorbed  dose  (D),  as 
defined  by  the  ICRU  (1980),  is  the  quotient  of  ds  by  dm,  where  ds  is  the  mean  energy 
imparted  by  ionizing  radiation  to  matter  of  mass  dm.  Secondly,  deposition  rate  is  important 
due  to  the  ability  of  living  organisms  to  repair  damage.  Repair  capability  may  be  limited  in 
part  by  dose  rate  and  the  type  of  ionizing  radiation.  The  spatial  distribution  of  energy 
deposition  at  the  microscopic  level  is  another  important  characteristic  that  will  be  discussed 
later. 

Space  Radiation  Environment 

The  space  radiation  environment  is  different  from  that  encountered  on  earth  in  both 
intensity  and  types  of  radiation.  The  earth  has  the  protective  magnetic  field  that  shields  the 
surface  from  charged  particle  radiations  as  well  as  the  attenuation  provided  by  the 
atmosphere.  Space  radiation  consists  of  three  components:  trapped  particle,  solar  particle, 
and  galactic  cosmic  radiations  (GCR). 

The  electron  belts  surrounding  the  earth  contain  particles  trapped  in  the  earth’s  magnetic 
field.  The  belts  consists  of  two  distinct  regions:  an  inner  region  that  extends  out  to  about 
1.8  x  104  kilometers  (km)  and  an  outer  region  beyond  that  of  the  inner  to  7.6  x  104  km 
(NCRP  1989).  Electron  intensities  in  the  outer  region  are  typically  about  10  times  those  in 
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the  inner  region,  and  within  both  regions  intensities  can  vary  drastically  over  short  periods 
of  time  (NCRP 1989).  For  the  inner  region,  exposure  to  protons  is  of  greater  importance 
than  electrons  for  manned  space  missions  (NCRP  1989).  At  450  km  altitude  (28.5° 
inclination),  the  majority  of  the  proton  fluence  has  energies  between  5  and  100  MeV 
(NCRP  1989).  There  can  be  considerable  variation  in  proton  fluence  with  inclination  angle 
to  the  earth’s  axis;  the  most  intense  region  being  the  South  Atlantic  Anomaly  that  is  located 
between  Africa  and  South  America 

Outside  the  protection  of  the  earth’s  magnetic  field,  the  environment  consists  of  solar 
and  galactic  particles.  The  galactic  cosmic  spectrum  is  present  isotropically  in  space  and  is 
from  source(s)  outside  the  solar  system.  Ninety-eight  percent  of  total  galactic  particle 
fluence  is  protons  and  heavier  ions,  with  2  %  electrons  and  positrons  (NCRP  1989).  For 
the  barions,  the  bulk  of  the  particle  fluence  have  energies  between  0. 1  to  10  GeV/nucleon 
and  among  these  the  majority  is  between  0. 1  to  1  GeV/nucleon  (NCRP  1989).  Table 
2-1  lists  the  fractional  abundance  (FA)  of  the  most  important  ions  and  their  FA  multiplied 

by  the  Z2  of  the  ion  (NCRP  1989).  For  galactic  cosmic  rays,  FA  •  Z2  is  a  good  indicator 

of  the  relative  dose  among  the  various  ions  as  collisional  stopping  power  is  closely 
proportional  to  Z2  and  the  energy  spectra  among  the  ions  is  similar  in  except  for  scaling 
factors.  Thus,  though  the  FA  of  Fe  is  2300  times  lower  than  that  of  protons  and  one-tenth 
that  of  C  and  O,  its  contribution  to  dose  is  greater  than  that  of  both  O  and  C,  and  about 


Table  2-1.  Fluence  Weighted  Relative  Fractional  Abundance  (FA) 
of  Galactic  Cosmic  Barions  and  FA  •  Z2  (NCRP  1989). 


Barion 

Fractional  Abundance 

FA  •  Z2 

Hydrogen 

0.87 

0.87 

Helium 

0.12 

0.48 

Carbon 

0.004 

0.13 

Oxygen 

0.004 

0.20 

Iron 

0.0004 

0.26 

7 


30  %  of  that  from  protons.  Ions  that  are  heavier  than  helium  are  referred  to  as  HZE 
particles  (high  Z  and  energy).  The  influence  of  the  sun  to  the  radiation  environment  outside 
the  earth’s  magnetosphere  is  highly  dependent  on  the  activity  of  the  sun.  Solar  winds  are 
created  by  ionized  gas  emissions  from  the  sun  and  carry  magnetic  fields.  These  magnetic 
fields  decrease  the  intensity  of  the  GCR  in  interplanetary  space.  Thus  in  years  of  lower 
solar  activity  (solar  minimum),  the  GCR  reaches  a  maximum,  while  GCR  minimums  are 
achieved  during  periods  higher  solar  activity.  In  addition  to  the  influence  of  solar  winds, 
the  sun  has  large  emissions  of  charged  particles  in  solar  particle  events  (SPE)  that  occur 
during  solar  flare  activity.  The  emissions  consist  of  protons,  helium,  and  heavier  nuclei, 
but  have  energy  spectra  lower  than  that  of  the  GCR.  Solar  particle  events  can  be  grouped 
into  two  categories:  ordinary  and  anomalously  large  events.  One  anomalously  large  event 
occurred  from  4-7  August,  1972  and  had  a  one-week  integrated  fluence  many  orders  of 
magnitude  greater  than  that  of  the  galactic  fluence.  These  events  are  rare  (NCRP 1989). 

In  addition  to  the  primary  particles  described,  secondary  particles  are  produced  through 
interactions  with  spacecraft  materials  and  personnel.  Secondary  particles  consist  primarily 
of  charged  fragmentation  and  spallation  products,  and  neutrons  from  nuclear  interactions . 
Spatial  Patterns  of  Energy  Deposition 

Introduction 

For  different  types  of  radiations  that  provide  similar  dose  and  dose  rates  to  biological 
systems,  the  primary  difference  among  the  radiations  in  reaching  many  biological  outcomes 
is  related  to  the  spatial  pattern  of  energy  deposition. 

In  1952,  Zirkle  et  al.  (1952)  introduced  the  concept  of  linear  energy  transfer  (LET) 
which  referred  to  the  linear  density  of  all  forms  of  energy  transfers  from  ionizing  radiation. 

Zirkle  introduced  the  concept  to  explain  differences  in  survival  of  spores  exposed  to  a- 

particles  and  x-radiations.  In  his  experiment,  the  high-LET  a-particles  exposed  spores  had 
a  lower  survival  than  those  exposed  to  the  x-radiation  (low-LET)  at  the  same  dose.  In 
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1962,  the  ICRU  (1962)  defined  linear  energy  transfer  (LET)  for  charged  particles  in  a 
medium  as  the  quotient  of  dEJdl,  where  dE^  is  the  average  energy  locally  imparted  to  the 
medium  by  a  charged  particle  of  specified  energy  in  traversing  a  distance  of  d /.  Within  this 
definition,  it  was  recognized  that  primary  charged  particles  could  transfer  significant  energy 
to  secondary  electrons  such  that  the  energy  carried  by  these  secondaries  would  impart  the 
majority  of  their  kinetic  energy  outside  of  the  location  of  the  primary  particle  track.  The 
term  “locally  imparted”  referred  either  to  a  maximum  radial  distance  from  the  particle  track 
or  a  maximum  discrete  kinetic  energy  transfer  to  a  secondary  electron  beyond  which  are  not 

included  in  LET.  Secondaries  with  high  kinetic  energy  are  generally  referred  to  as  8-rays 

and  are  created  by  hard  or  “knock-on”  collisions  with  the  primary  where  the  classical 
impact  parameter  (of  the  primary  and  atom  that  possesses  the  secondary  electron)  is  on  the 
order  of  atomic  dimensions.  Subsequent  ICRU  Reports  (1970,1980)  defined  LET 

(restricted  linear  collision  stopping  power  -  of  charged  particles  in  a  medium  as  the 
quotient  of  dEA  by  dl,  where  dEA  is  the  mean  energy  loss  due  to  collisions,  with  kinetic 

energy  transfers  less  than  some  specified  value  A,  and  dl  is  the  distance  transversed  by  the 

particle.  For  the  remainder  of  this  document,  LET  (and  LETJ  will  interchangeably  refer  to 
the  total  (or  unrestricted)  LET  that  accounts  for  all  collisional  losses  regardless  of  energy. 
For  simplicity  an  L  will  be  used  in  place  of  LET. 

Many  researchers  over  the  years  have  investigated  the  spatial  patterns  of  energy 

deposition  for  protons,  a-particles,  and  heavy  ions.  The  work  was  prompted  primarily  by 
concern  over  radiation  hazards  from  neutrons  encountered  in  the  commercial  nuclear  power 
industry,  a-particles  from  natural  and  man-made  sources,  ions  encountered  in  space 

operations,  and  ions  used  in  medicine.  The  investigations  have  consisted  of: 

a)  measurements  of  radial  dose  distributions  and  formulation  of  empirically  based 
models. 
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b)  Monte  Carlo  simulations  of  energy  deposition,  and 

c)  microdosimetric  measurements. 

Track  Structure  Models 
Katz  Model 

Katz  and  co-workers  (Butts  and  Katz  1967,  Katz  eal.  1971)  provided  one  of  the  first 
track  structure  models  for  heavy  ions  based  on  a  6-ray  distribution  formula  and  an 

extrapolated  range-energy  relation  for  electrons.  The  radial  distribution  of  dose,  calculated 
from  these  relations,  follows  an  approximate  (1/r2)  relationship,  where  r  is  the  radial 
distance  from  the  primary  track.  In  biological  studies  with  dry  enzymes  and  viruses,  Katz 
found  their  response  to  heavy  ions  irradiation  to  be  predicted  solely  on  the  distance  from 
the  primary  track. 

Chunxiang  etal.  (1985)  compared  the  Katz  model  of  radial  distribution  of  dose  to  four 
sets  of  experimental  data  collected  by  Varma  et  cd.  ( 1975, 1976, 1980a,  and  1980b)  and 
another  model  of  radial  dose  distribution  proposed  by  Fain  et  al.  ( 1974).  Comparison  of 
the  model  to  experimental  data  provided  close  agreement  for  helium  (Z  =  2),  oxygen  (Z  = 

8),  neon  (Z  =  10),  and  bromine  (Z  =  35)  ions  over  a  velocity  range  of:  0.072  <  (3  <  0.7. 
For  comparison  with  the  Fain  model,  there  was  close  agreement  for  carbon  (Z  =  6  and 
P  =  0.065)  and  iron  (Z  =  26  and  p  =  0.41)  ions. 

The  Katz  model  was  used  by  NASA  (Cucinotta  1991)  tb  predict  biological  damage  to 
mammalian  cell  cultures  from  galactic  cosmic  rays.  Cucinotta  etal  (1995)  have  modified 
input  parameters  to  the  Katz  model  to  improve  accuracy  of  die  model  both  near  and  far 
from  the  ion  trajectory  where  the  model  had  the  greatest  uncertainties. 
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Chatteriee  Model 

Chatteijee  and  co-workers  (Chatteijee  etal.,  1973;  Chatteijee  and  Shaefer,  1976;  and 
Chatteijee,  1993)  developed  a  model  similar  to  Katz’s  in  the  (1/r2)  relationship  of  dose  to 
radial  distance  from  the  primary  track.  The  model  is  simpler  than  Katz’s  model  having 
only  three  input  parameters:  L,  penumbra  radius,  and  core  radius.  The  model  divides 
radial  dose  in  two  distinct  regions:  an  inner  region  of  very  dense  ionization  called  the  core 
that  arises  primarily  from  deposition  of  energy  from  electrons  created  from  soft  collisions, 

and  an  outer  region  that  encompasses  the  range  of  6-rays  produced  along  the  track  that  is 

called  the  penumbra.  Energy  deposition  is  assumed  not  to  occur  outside  the  penumbra. 

The  model  is  independent  of  particle  charge.  The  penumbra  and  core  radii  are  the  same  for 
particles  of  the  same  total  kinetic  energy  per  nucleon  and  are  dependent  only  on  velocity. 
Like  the  Katz’s  model,  Chatteijee’s  model  closely  agreed  with  microdosimetric 
measurements  of  600  MeV/amu  Fe  ions  that  were  stripped  of  their  electrons  (Metting  etal. 
1986).  Toburen  etal.  (1990)  provided  microdosimetric  measurement  of  partially  stripped 

U  and  Ge  ions  that  had  a  maximum  velocity  of  0  =  0. 189.  The  mean  energy  deposited 

was  in  close  agreement  with  the  Chatteijee  model,  but  had  energy  deposition  beyond  the 
penumbra  radius  of  the  model.  One  advantage  of  the  Chatteijee  model  over  the  Katz’s 
model  is  that  it  incorporates  energy  conservation  and  can  be  used  to  calculate  the  fraction  of 
L  absorbed  in  volumes  surrounding  the  particle  track.  For  this  reason,  this  model  will  be 
used  for  calculations  in  this  report.  Though  some  researchers  measured  energy  deposition 
beyond  the  penumbra  radii  specified  in  the  model,  this  inaccuracy  of  the  model  should  have 
no  influence  on  the  calculations  used  in  this  report  because  they  were  made  for  radii  well 
within  the  penumbra  Appendix  A  summarizes  the  calculations  performed  with  this  model. 
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Monte  Carlo  Simulations 


Several  computer  codes  based  on  the  Monte  Carlo  method  are  currently  used  to 
simulate  the  transport  of  charged  particles  through  media  (Kraft  etal.,  1992  and  Wilson 
etal.,  1994).  The  codes  generally  accomplish  the  same  task,  but  may  be  subtly  different  in 
the  extent  to  which  tracks  are  followed  (Rossi  and  Zaider,  1996).  Ottolenghi  etal.  (1995) 
have  adapted  a  Monte  Caiio  track  structure  code  to  simulate  strand  breaks  produced  in 

DNA  from  a-particles  and  protons. 

The  Monte  Carlo  method,  being  based  on  the  stochastic  characteristics  of  energy 
deposition,  more  closely  approximate  physical  reality  than  track  structure  models  that 
provide  only  mean  energy  deposition.  Under  some  conditions,  mean  values  may  be 
sufficient  in  describing  energy  deposition.  For  large  sites  in  and  near  the  track  of  charged 
particles,  or  for  multiply  charged  ions  with  high  ionization,  the  variance  in  energy 
deposition  may  be  small.  But,  for  sensitive  volumes  distant  from  the  track,  small  volume 
sites,  or  singly  charged  ions,  the  stochastic  nature  of  energy  deposition  becomes  more 
important,  making  mean  values  less  meaningful. 

Microdoslmetrv 

1.  Micrometer  Measurements.  Rossi  and  coworkers  (1955a,  1955b)  developed  the 
first  device  to  measure  energy  deposition  in  simulated  volumes  of  tissue  with  diameters  on 
the  order  of  1  ptm.  These  devices  are  referred  to  as  tissue-equivalent  proportional  counters 
(TEPC).  Microdosimetric  techniques  have  been  used  in  radiation  protection,  radiobiology, 
medical  therapy,  and  microelectronics.  Volume  size  is  chosen  relevent  to  structure  of 
biological  interest,  for  example,  cell  nuclei  of  diameters  typically  between  1  and  10  jn n. 

The  technique  has  been  used  for  characterizing  radiation  exposure  to  NASA  personnel  on 
manned  space  flights,  with  simulated  diameters  in  the  micron  range  (Badhwar  etal.  1992). 

2.  Conventional  TEPCs.  A  conventional  TEPC  consists  of  a  solid  wall 
surrounding  a  gas-filled  cavity  (sensitive  volume).  The  most  commonly  utilized  shape  of 
the  sensitive  volume  has  been  the  sphere  since  it  possesses  complete  symmetry  with  respect 
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to  isotropic  and  directional  radiation  sources  (ICRU  1983).  Another  shape  commonly 
implemented  is  the  cylinder  and  is  the  type  used  for  dosimetry  on  manned  space  flights. 
Cylindrical  devices  are  easier  to  construct  and  generally  are  structurally  more  rugged, 
which  makes  them  more  attractive  for  use  on  spacecraft  Detectors  constructed  of  solid 
tissue-equivalent  plastic  walls  to  define  the  sensitive  volume  are  easy  to  construct  and 
rugged.  But,  wall-less  detectors  (where  the  medium  surrounding  the  sensitive  volume  is  of 
the  same  density  and  composition),  are  the  most  desirable  because  radiation  equilibrium  is 
preserved  as  generally  is  the  case  in  tissue.  Distortions  created  by  the  solid  wall,  termed 
“wall  effects,”  are  the  primary  focus  of  this  research. 

3.  Limitations  of  Microdosimetrv.  The  measurement  of  energy  deposition  on  the 
micrometer  scale  has  been  useful  for  predicting  a  wide  range  of  biological  experimental  data 
(Kellerer  and  Rossi,  1972).  Limitations  in  the  predictive  ability  of  microdosimetry  was 
apparent  for  very  low-energy  X-rays  (Cox  et  al.  1977,  Goodhead  et  al.  1978).  For  very 
low-energy  X-rays,  ionizations  are  closely  spaced  and  resultant  biological  effects  could  not 
be  predicted  by  energy  deposition  on  the  micrometer  scale.  Experiments  of  this  type  stress 
the  importance  of  energy  deposition  on  the  nanometer  scale  (in  addition  to  that  on  the 
micrometer  scale)  in  predicting  biological  effects  (Kliauga  1994). 

4.  Nanodosimeters.  Kliauga  (1990)  evaluated  the  response  of  a  ultraminiature 
walled  cylindrical  TEPC  to  photons  and  neutrons  down  to  site  sizes  of  10  nm.  More 
recently  Kliauga  (1994)  evaluated  the  response  of  a  nearly  wall-less  miniature  cylindrical 
TEPC  to  heavy  ions.  Acceptable  performance  was  achieved  for  simulated  site  sizes  down 
to  20  nm. 

5.  Dosimetry  on  NASA  Missions.  Based  on  the  above  discussion,  the  ideal  TEPCs 
would  be  of  nearly  wall-less  designs  and  incorporate  devices  that  simulated  nanometer  and 
micrometer  site  sizes.  For  many  reasons  this  has  yet  to  occur.  Wall-less  devices  are  not  as 
rugged  as  walled  devices  and  may  not  be  able  to  withstand  G-stresses  encountered  during 
lift-off.  To  properly  create  radiation  equilibrium  with  a  wall-less  detector,  the  detector 
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container  would  have  to  be  orders  of  magnitude  larger  than  detector  sensitive  volume 
(Metting  et  al.  1988).  This  is  not  practical  for  spacecraft.  Miniature  detectors  that  are 
necessary  to  simulate  nanometer  site  sizes  are  difficult  to  construct  and  not  as  rugged  as 
devices  that  simulate  micrometer  site  sizes.  Distortions  created  by  the  wall  effect  (to  be 
discussed  in  detail  later  in  this  Chapter)  will  be  more  pronounced  for  nanodosimeters  than 
microdosimeters  in  the  field  of  high  velocity  ions. 

6.  Quantities  and  Units.  The  measured  ionization  produced  in  a  TEPC  is 

proportional  to  the  energy  deposited  in  the  cavity,  s.  In  the  field  of  microdosimetry,  the 
concept  of  lineal  energy  (y)  is  more  commonly  used.  ICRU  ( 1983)  defined  lineal  energy 
as  the  quotient  of  £  by  l  ,  where  l  is  the  mean  chord  length  of  /^-random  particle  traversals 

through  that  volume.  Lineal  energy  is  a  random  or  stochastic  quantity  and  is  normally 
expressed  in  units  of  keV  /an"1.  By  recording  the  number  of  events  within  event  size 
(energy)  bins  of  some  defined  value  and  dividing  each  by  the  total  number  of  events  among 
all  bins,  a  discrete  probability  density  function  (pdf),  f(y),  in  lineal  energy  can  be  defined. 
These  distributions  or  functions  are  called  single  event  spectra.  The  first  moment  of  f(y), 

yp ,  is  called  the  frequency  mean  lineal  energy,  is  a  nonstochastic  property  of  f(y),  and  is 
defined  as: 

__  00 

yf  -Jyf(y)dy.  (2-i) 

o 

The  second  moment  of  f(y),  yD,  is  called  the  dose  mean  lineal  energy,  is  a  nonstochastic 
property  of  f(y),  and  is  defined  as: 

00 

/  y2  f(y)  dy  ^ 

yD  =  - — = - =  /  y  d(y)  dy ,  (2-2) 

yF  o 

where  d(y)  =  ■=¥—  f(y). 
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Radiation  Dosimetry  and  Radiation  Protection 

Relative  biological  effectiveness  (RBE)  is  a  concept  used  in  radiobiology  to  compare 
the  dose  of  a  test  radiation  with  the  dose  of  250  kVp  x-rays  required  to  achieve  a  similar 
biological  outcome.  RBE  is  defined  as  the  quotient  of  the  absorbed  dose  of  the  test 
radiation  by  that  of  250  kVp  x-rays.  While  the  concept  of  RBE  may  appear  to  be 
straightforward,  it  is  difficult  to  apply  directly  to  radiation  protection  because  RBE  changes 
with  dose,  dose  rate,  mixed  radiation  fields,  and  specific  biological  outcomes. 

For  the  purpose  of  radiation  protection,  quality  factor  (Q),  has  been  assigned  to 
radiations  of  different  types  to  account  for  the  varying  ability  to  cause  detrimental  radiation 
effects.  Risks  for  detrimental  effects  from  radiation  exposure  are  related  to  the  dose  of 
radiation  received  modified  by  the  Q  of  the  radiation.  The  International  Commission  on 
Radiological  Protection  in  Publication  26  (ICRP  1977)  provided  a  recommended 
relationship  between  L  and  Q.  The  recommendation  was  changed  in  Publication  60  (ICRP 
1990)  where  Q  was  replaced  by  w„  the  radiation  weighting  factor.  In  this  report,  Q  will  be 
used  interchangeably  for  both  set  of  recommendations.  The  recommendations  of  both 
publications  is  given  in  Fig.  2-1. 

For  charged  particles  of  similar  L,  but  with  different  charge  and  velocity,  there  can  be 
large  differences  is  the  spatial  pattern  of  energy  deposition  and  therefore  biological  effects. 
Also,  L  is  a  quantity  that  is  difficult  to  measure.  ICRU  (1986)  has  recommended  a 
relationship  between  Q  and  y.  The  relationship  is  provided  in  Fig.  2-2  for  1  pern  diameter 
spheres.  Because  lineal  energy  represents  discrete  energy  deposition,  from  a  biological 
standpoint  it  may  be  more  meaningful  than  L.  A  TEPC  can  be  designed  to  measure  y 
directly  and  therefore  provide  direct  measurements  of  Q. 

NCRP  ( 1989)  recommended  dose  equivalent  limits  for  space  travel.  In  the  document, 
it  was  recognized  that  the  ICRU  (1986)  recommended  Q  values  for  neutrons  should  be 
implemented  for  radiation  exposures  in  space.  But,  the  report  did  not  specifically 

recommend  determination  of  Q  directly  from  measured  microdosimetric  spectra.  NASA  is 


15 


Figure  2-1.  ICRP  26  and  60  Quality  Factor  Relationships  to  L 


Figure  2-2.  ICRU  40  Quality  Factor  Relationship  to  y. 
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currently  determining  Q  under  the  assumption  that  y  =  L  and  with  the  ICRP  recommend¬ 
ed  relationship  between  Q  and  L.  Plastic  detectors  CR-39,  cellulose  nitrate,  and  Lexan 
polycarbonate  have  been  used  by  NASA  on  past  missions  to  measure  L  (NCRP  1989). 
NASA  has  collected  microdosimetric  spectra  with  walled  cylindrical  TEPCs  on  manned 

space  flights.  Calculations  of  Q  have  been  made  from  these  spectra  under  the  assumption 
that  y  is  equivalent  to  L  for  discrete  energy  deposition  events.  Because  energy  deposition 
is  influenced  by  the  pathlength  of  the  ion  through  the  detector,  clearly  the  equivalence  is  not 
true  on  an  individual  event  basis.  NASA  has  been  criticized  for  this  practice  (Dicello  et  al. 

1991).  However,  it  has  not  been  determined  whether  calculations  of  Q  are  significantly 
different  from  the  known  values  for  high  velocity  ions  because  the  wall  effects  for  these 
detectors  have  not  been  firmly  established. 

Ricourt  et  al.  (1981)  recommended  that  dose  equivalent,  H,  be  determined  from 
measurements  of  D(y)  using  the  following  relationship: 

H  =  /Q(L  =  %y)  D(y)  dy  (2-3) 

‘The  substitution  of  8/9  y  for  L  is  derived  from  the  moments  analysis  of  y  in  the  limit  when 
the  simulated  tissue  diameter  approaches  zero,”  according  to  Borak  and  Stinchcomb 
(1982).  Borak  and  Stinchcomb  (1982)  found  that  the  substitution  provided  an  almost  exact 
estimate  of  the  actual  Q  for  neutrons  of  energy  from  thermal  to  4  MeV  if  the  simulated 
diameter  of  the  tissue  is  zero.  But,  as  the  simulated  diameter  increased,  calculated  values  of 
Q  decreased. 

Radiation  risk  coefficients  based  on  a  fluence  based  system  were  first  discussed  by 
Curtis,  Dye,  and  Sheldon  (Wilson  et  al.  1990).  The  system  used  risk  coefficients  for  a 
particular  endpoint  based  on  the  fluence  of  individual  particles  in  the  radiation  field. 
Dosimetry  under  this  system  would  require  particle  identification  and  velocity.  Practically, 
this  could  be  accomplished  with  a  particle  telescope  that  contained  a  series  of  silicon 
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detectors.  Unfortunately,  a  system  of  this  type  could  not  record  dose  from  indirectly 
ionizing  radiations  that  are  present  in  the  space  environment. 

The  purpose  of  this  work  is  to  characterize  the  wall  effects  from  56Fe  incident  on  walled 

TEPCs.  This  data  can  be  used  to  determine  an  appropriate  method  for  calculation  of  Q . 

Wall  Effects 
General 

Wall  effects  caused  by  the  difference  in  density  between  two  media  were  first  described 
by  Failla  and  Failla  ( 1960)  and  later  by  Rossi  ( 1968).  Oldenburg  and  Booz  ( 1970) 
classified  these  effects  into  different  categories.  Kellerer  (1971a,  1971b)  provided  detailed 

theoretical  analysis  for  heavy  charged  particles,  where  the  V-effect  and  6-ray  effect  are 

most  important. 

6-Rav  Effect 

The  effect  is  most  important  for  high  velocity  charged  particles  where  6-ray  secondaries 

can  have  significant  energy  and  range.  The  effect  is  dependent  on  the  trajectory  of  the 
charged  particle  with  respect  to  the  cavity  of  the  detector.  Three  types  of  trajectories  will  be 
described. 

Figure  2-3  illustrates  the  6-ray  effect  for  a  walled  cavity  as  compared  to  energy 

deposition  in  homogenous  media.  For  both  cases,  the  primary  charged  particle  traverses 
the  sensitive  volume  (cavity).  These  are  referred  to  as  “direct  events.”  In  the  case  of  the 
walled  cavity,  the  energy  deposited  will  be  the  sum  of  that  from  the  primary  particle  and 

from  the  6-ray(s)  created  in  the  wall.  However,  in  homogenous  media,  the  distance 

between  die  location  of  6-ray  creation  and  the  entrance  point  of  die  detector  may  be 

sufficiently  large  such  that  only  one  of  the  particles  may  traverse  the  sensitive  volume.  As 
shown  in  the  figure  for  homogenous  media,  the  equivalent  amount  of  energy  that  was 
deposited  in  the  walled  cavity,  is  represented  by  two  discrete  energy  deposition  events. 
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For  highly  charged  particles,  a  single  event  in  a  walled  cavity  may  be  represented  by  the 
sum  of  numerous  discrete  energy  deposition  events. 


Figure  2-3.  Diagram  Illustrating  the  6-Ray  Effect  in  a  Walled  Cavity  as  Compared 
to  Equivalent  Energy  Deposition  in  Homogenous  Media.  [In  the  walled  cavity,  energy 
deposition  is  represented  by  a  single  event,  but  for  homogenous  media  the  same  energy 
deposition  is  represented  by  two  discrete  events.]  [Adapted  from  Kellerer  (1971b)]. 


a)  walled  cavity  b)  homogenous  media 

Figure  2-4  provides  another  illustration  of  this  effect  for  a  primary  charged  particle 
traversing  a  spherical  cavity  along  the  longest  chord  length.  Evaluation  of  the  difference  in 
energy  deposition  between  a  walled  cavity  and  homogenous  media  for  these  types  of 
trajectories  have  not  been  accomplished. 

Figure  2-4.  Diagram  Illustrating  the  Difference  in  Energy  Deposition  in  a  Walled  Cavity 
and  Homogenous  Media  for  Direct  Events  [Adapted  from  Kellerer  (1971a)]. 


a)  walled  cavity  b)  homogenous  media 

Indirect  events  are  those  where  the  primary  charged  particle  passes  outside  of  the 
sensitive  volume.  Indirect  events  for  the  walled  cavity  can  either  enhance  or  suppress 


energy  deposition  as  compared  to  that  in  homogenous  media  dependent  on  the  location  of 
the  particle  trajectory. 

Figure  2-5a  illustrates  an  indirect  event  where  the  trajectory  of  the  primary  particle 
traverses  outside  of  the  cavity  but  in  the  wall  of  the  detector.  Energy  deposition  for  these 
events  will  be  much  greater  as  compared  to  that  in  homogenous  media  (Fig.  2-5b). 


Figure  2-5.  Diagram  Illustrating  Energy  Deposition  in  a  Walled  Cavity  and 
Homogenous  Media  for  Indirect  Events  Traversing  Outside  the  Sensitive  Volume. 


a)  walled  cavity 


b)  homogenous  media 


The  frequency  and  magnitude  of  these  events  depends  on  the  range  of  the  6-rays,  proximity 
of  the  primary  particle  to  the  cavity/wall  interface,  simulated  site  volume,  and  the  Z  of  the 
primary  particle.  The  range  of  the  6-rays  is  directly  related  to  the  velocity  of  the  primary 

particle.  Kellerer  (1971a)  calculated  the  absorbed  dose  in  homogenous  media  for  proton 
recoils  and  heavy  charged  particle  (HCP)  of  initial  energies:  1  to  30  MeV/nucleon  in 
spherical  sites  of  diameter  from  0.01  to  10  ptm.  Absorbed  dose  was  determined  separately 
for  direct  and  indirect  events.  For  a  30  MeV/nucleon  HCP,  the  fraction  of  absorbed  dose 
from  indirect  events  was  as  high  as  0.35  for  0.01  ptm  diameter  sites  and  lower  than  0.09 
for  10  Jim  diameter  sites.  Theoretical  calculations  haven’t  been  published  for 
1  GeV/  nucleon  HCP. 

Figure  2-6a  illustrates  an  indirect  event  where  the  trajectory  of  the  primary  is  outside  of 
the  cavity  wall.  If  the  wall  is  thicker  than  the  range  of  the  maximum  energy  6-ray,  then 
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there  will  be  no  deposition  of  energy  in  the  cavity  from  this  events,  neglecting  x-rays  from 
Bremsstrahlung.  If  the  cavity  wall  is  thinner  than  the  penetration  range  of  the  maximum 

energy  6-ray,  then  there  is  the  potential  for  energy  deposition  in  the  cavity  from  this  event, 

but  the  net  result  of  this  type  of  indirect  event  is  a  suppression  of  energy  deposition  as 
compared  to  that  in  homogenous  media  (Fig.  2-6b). 

Figure  2-6.  Diagram  Illustrating  Energy  Deposition  in  a  Walled  Cavity  and  Homogenous 
Media  for  Indirect  Events  with  Trajectories  Outside  the  Wall  of  the  Walled  Cavity. 


a)  walled  cavity  b)  homogenous  media 

A  description  of  three  trajectories  have  been  given  to  characterize  the  6-ray  effect  for 

charged  particles.  For  two  of  the  descriptions,  energy  deposition  was  enhanced  over  the 
homogenous  case,  while  for  the  other  it  was  suppressed.  For  a  uniformly  irradiated  walled 

cavity  TEPC,  the  net  result  in  single  event  spectra  is  an  increase  in  yp  and  yD  as 
compared  to  homogenous  media. 

Reentry  Effect 

Figure  2-7a  illustrates  the  reentry  effect  for  an  electron  in  a  walled  cavity.  From  the 
diagram,  the  electron  traverses  the  cavity,  enters  the  wall,  and  reenters  the  cavity  a  second 
time.  For  homogenous  media  (Fig.  2-7b),  the  electron  will  traverse  the  cavity  only  once. 
For  the  event  depicted,  the  net  effect  in  the  walled  cavity  is  an  increase  in  energy  deposition 
compared  to  that  of  homogenous  media.  The  reentry  effect  is  only  significant  for  electrons 
because  they  are  the  only  charged  particle  with  significant  curvature  of  path.  Because 
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photons  deposit  energy  through  transfer  of  energy  to  electrons,  it  is  applicable  to  photons 
as  well. 


Figure  2-7.  Diagram  Illustrating  the  Reentry  Effect  in  a  Walled  Cavity. 
[Energy  deposition  in  a  walled  cavity  is  represented  by  the  electron  traversing 
the  walled  cavity  twice,  while  in  homogenous  media  it  is  represented  as  two 
discrete  energy  deposition  events.]  [Adapted  from  Kellerer  (1971a)]. 


i  homogenous  media 

Scattering  Effect 

The  third  categoiy  of  wall  effect  is  called  the  scattering  effect.  The  effect  is  applicable 
to  indirect  ionizing  radiations  (i.e.  neutrons  and  photons).  The  effect  is  illustrated  for  a 
walled  cavity  in  Fig.  2-8a,  where  two  scattering  events  occur  within  close  proximity  of 
each  other  and  the  scattered  charged  particles  both  deposit  energy  in  the  cavity. 


a)  walled  cavity 


Figure  2-8.  Diagram  Illustrating  the  Scattering  Effect  in  a  Walled  Cavity. 
[Energy  deposition  in  the  walled  cavity  is  represented  by  the  sum  of  energy  deposited 
by  the  two  charged  particle  recoils,  but  in  homogenous  media  it  would  be  represented 
by  two  discrete  energy  deposition  events.]  [Adapted  from  Kellerer  (1971a)]. 
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In  homogenous  media  however,  deposition  of  energy  by  both  secondaries  in  the  same 
active  volume  is  extremely  improbable.  The  effect  is  dependent  on  the  relative  densities  of 
the  cavity  and  wall  materials,  thickness  of  cavity  wall,  and  interaction  cross-sections  of  the 
primaries.  The  net  result  of  this  effect  is  enhanced  energy  deposition  as  compared  to  that  in 
homogenous  media. 

V-Effect 

The  fourth  type  of  wall  effect  is  the  V-effect  and  can  occur  in  both  elastic  and  inelastic 
(absorption)  nuclear  interactions.  The  effect  for  a  walled  cavity  is  illustrated  in  Fig.  2-9a 

Figure  2-9.  Diagram  Illustrating  the  V-Effect  Effect  in  a  Walled  Cavity. 

[Energy  deposition  in  the  walled  cavity  is  the  sum  of  that  from  both  charged 
particles  trajectories,  but  in  homogenous  media  it  would  be  represented  by 
two  discrete  energy  deposition  events.]  [Adapted  from  Kellerer  (1971a)]. 


€C 

a)  walled  cavity 

where  the  nuclear  interaction  occurs  in  close  proximity  of  the  cavity.  The  “V”  formed 
illustrates  the  trajectory  of  both  the  primary  and  a  secondary  nuclei  where  both  particles 
contribute  to  energy  deposition  in  the  cavity.  In  homogenous  media  (Fig.  2-9b),  the  two 
events  would  not  likely  deposit  energy  in  the  same  site.  The  net  effect  in  a  walled  cavity  is 
enhanced  energy  deposition  compared  to  homogenous  media  where  the  energy  deposition 
would  be  represented  by  to  discrete  energy  deposition  events.  The  example  provided  here 
is  the  simplest  case  with  only  one  secondary  nucleus  produced.  In  violent  nuclear 
collisions,  especially  seen  in  HZE  particle  collisions,  there  can  be  many  secondary  nuclei 
produced.  Appendix  C  provides  details  of  nuclear  interactions  for  HZE  particles  with 
specific  data  for  S6Fe  ions. 


b)  homogenous  media 
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Wall  Effects  Characterized  by  Experimental  Measurements 


General 

A  number  of  investigators  experimentally  evaluated  the  difference  in  response  of  walled 
and  wall-less  TEPC  designs.  The  bulk  of  the  data  was  collected  in  the  1970s.  Because  the 
focus  of  this  thesis  is  on  the  wall  effects  of  HZE  particles  a  more  complete  summary  of 
relevant  work  will  be  provided,  while  few  other  radiation  types,  only  a  small  fraction  of  the 
relevant  work  will  be  summarized. 

Photons 

Photons  were  the  primary  source  of  radiation  in  radiation  biology  and  medical  therapy 
during  the  early  years  of  these  disciplines.  As  a  consequence,  a  considerable  amount  of 
experimental  microdosimetry  was  accomplished  on  these  fields.  Braby  and  Ellett  (1971) 
compared  single  event  spectra  from  spherical  solid-  and  grid-walled  TEPCs  in  the  field  of 
250  kVp  x-rays  and  a  60Co  source.  The  devices  simulated  volumes  of  tissue  with 
diameters  from  1  to  6  jim.  The  result  of  the  work  is  summarized  in  Fig.  2-10  in  terms  of 
the  mean  event  size  of  the  single  event  spectra.  From  the  plot,  mean  event  size  decreases  as 
the  simulated  site  diameter  increases  for  both  radiation  and  detector  types.  The  difference 
in  the  mean  event  size  for  the  two  detectors  in  the  field  of  the  250  kVp  x-rays  is 
considerable  with  the  maximum  ratio  of  1.52  at  2  pim  and  minimum  ratio  of  1.36  at  6  ptm. 
Few  the  irradiations  with  60Co,  the  difference  between  the  two  detectors  was  considerably 
less,  with  the  maximum  ratio  being  1.2.  The  authors  attributed  the  majority  of  the  wall 

effects  to  the  reentry  effect,  but  noted  that  for  the  60Co  y-rays,  a  significant  fraction  of  the 
effect  may  be  due  to  the  6-ray  effect. 
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Figure  2-10.  Frequency  Mean  Event  Size  from  Single  Event  Spectra  of  60Co  and  Filtered 
250  kVp  X-Rays  Collected  with  Solid-  and  Grid-Walled  TEPCs  (Braby  and  Ellett  1971). 
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Neutrons 

Oldenburg  and  Booz  (1970)  provided  theoretical  calculations  and  measurements  of 
6  MeV  neutrons  for  simulated  spherical  volumes  6.5  ptm  in  diameter.  For  events  with 
energy  deposition  below  20  keV,  the  calculated  spectrum  for  a  walled  cavity  decreased  with 
decreasing  neutron  energy,  while  the  measured  spectrum  increased.  This  difference  was 

believed  to  be  due  to  the  neglect  of  6-ray  production.  For  events  with  medium  energy 


deposition,  the  measured  and  calculated  spectrums  were  nearly  equivalent.  But  for  events 
with  high  energy  deposition,  event  frequency  in  the  calculated  spectrum  was  higher  than 
that  of  the  measured  spectrum.  The  calculated  single  event  spectra  for  the  walled  and  wall¬ 
less  TEPCs  were  in  close  agreement.  For  events  with  energy  deposition  from  20  to  55  keV 
(y:  4.6  - 12.7  keV/^m),  there  was  a  2  %  difference  in  event  frequency,  whereas  for  those 
with  energy  deposition  from  55  to  370  keV  (y:  12.7  -  85.4  keV/ptm)  there  was  no 
measurable  difference.  For  events  with  energy  deposition  greater  than  370  keV,  the  walled 
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TEPC  had  an  event  frequency  9  %  greater  than  the  walless  TEPC.  Since  theoretical 
calculations  had  produced  only  17  V-effect  events  out  of  237,000  runs  (0.007  %),  the 
authors  assumed  that  most  of  the  observed  wall  effects  must  be  from  multiple  scattering. 
Since  the  single  event  spectrum  for  the  walled  TEPC  had  exaggerated  the  frequency  of  high 
energy  deposition  events  in  comparison  to  that  measured,  it  was  concluded  that  the 
calculated  response  for  high  energy  events  was  biased  high  and  therefore,  the  expected 
response  of  walled  and  wall-less  detectors  would  be  nearly  equivalent.  Kellerer  (1971b) 
concluded  that  the  fraction  of  energy  deposition  involved  in  wall  effects  was  small  and 
would  have  an  insignificant  effect  on  the  frequency,  f(y),  and  dose,  D(y),  single  event 
spectra  for  neutron  energies  below  10  MeV. 

Heavy  Charged  Particles 

Theoretical  microdosimetric  spectra  from  high  velocity  charged  particles  are  not  well 
established  because:  1)  distributions  of  highly  energetic  6-rays  are  not  known  to  a  high 

degree  of  certainty  (ICRU  1995)  and  2)  energy  deposition  from  recoil  nuclei  from  nuclear 
interactions  have  not  been  characterized.  There  is  limited  experimental  work  on  the  wall 
effects  from  HZE  particles  because  of  limited  access  to  and  the  expense  of  operating 
particle  accelerators. 

Glass  and  Braby  (1969)  measured  a  broad  beam  of  5.3  MeV  a-particles  with  a  wall¬ 
less  spherical  counter.  Theoretical  calculations  were  made  under  the  assumption  that  all 
energy  transfers  were  locally  deposited.  Comparison  of  the  two  spectra  clearly 

demonstrated  the  importance  of  6-ray  contributions  to  the  single  event  spectra  and  the  need 
for  wall-less  detectors  in  the  fields  of  heavy  charged  particles. 

Gross  etal.  (1970)  measured  a  broad  beam  of  a-particles  from  an  241Am  source  with  a 

spherical  wall-less  TEPC  simulating  sites  0.5, 1,  and  2  ptm  in  diameter.  In  comparison  of 
the  single  event  spectra,  the  authors  noted  that  the  most  probable  value  decreased  as  the 
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simulated  site  diameter  became  smaller.  The  most  probable  value  of  the  spectra 
corresponds  to  a-particles  traversals  with  the  longest  pathlength  through  the  sensitive 
volume.  The  decrease  in  the  most  probable  value  is  due  to  a  decrease  in  the  amount  of 
energy  deposited  by  5-rays.  Similarly,  as  the  simulated  site  diameter  decreased,  the 
fraction  of  absorbed  dose  from  indirect  events  increased.  The  authors  compared  the 
detector  response  from  5.44  (P  =  0.053)  to  8.8  MeV  (P  =  0.069)  a-particles  and  found  the 

relative  contribution  to  absorbed  dose  from  6-rays  created  by  indirect  events  to  be  higher 
for  the  a-particles  of  higher  energy. 

Rodgers  etal.  (1973)  compared  the  response  of  spherical  walled  and  wall-less  (92  % 
transparent)  TEPCs  simulating  sites  2  pim  in  diameter  in  the  field  of  279  MeV/nucleon  14N 

ions  (P  =  0.638).  The  theoretical  maximum  6-ray  energy  is  700  keV,  significantly  higher 

than  that  of  5.44  MeV  a-particles  (2.9  keV).  The  authors  collected  spectra  for  the  two 

detectors  at  various  depths  in  a  phantom.  Data  with  minimum  absorber  material  upstream 
of  the  detectors  indicated  that  the  mean  lineal  energy  of  the  walled  detector  was  about  30  % 
higher  than  the  wall-less  detector.  The  difference  between  the  response  of  the  two 
detectors  became  less  apparent  when  the  velocity  of  the  particles  was  reduced  through 
additional  absorber  material.  Near  the  Bragg  peak,  the  difference  in  the  dose  mean  lineal 
energy  between  the  two  spectra  was  less  than  10  %.  ICRU  (1983)  discussed  the  findings 
of  this  work  and  compared  the  single  event  dose  distributions  for  each  detector  without 
absorber  material.  The  most  probable  value  of  each  distribution  was  nearly  the  same.  The 
walled  counter  was  deplete  in  low-energy  deposition  events  that  were  present  for  the  wall¬ 
less  detector.  But,  the  walled  detector  had  a  significantly  higher  number  of  high-energy 
deposition  events  as  compared  to  the  wall-less.  High-energy  deposition  events  were 
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attributed  to  “spallation  and  a  contribution  from  6-rays  at  low  lineal  energies”  (ICRU 
1983). 

Dicello  etal.  (1991)  found  similar  results  in  the  comparison  of  spectra  collected  with 
walled  and  wall-less  detectors  in  beams  of  bare  Fe,  Ar,  and  Ne  with  respective  energies  of 

535, 570,  and  557  MeV/nucleon  (f3:  0.772  -  0.784).  In  this  work,  direct  comparisons  of 

frequency  and  dose  mean  lineal  energy  values  were  not  made  because  the  geometry  of  the 
detectors  was  different. 

Kliauga  etal.  (1978)  measured  the  response  of  two  spherical  wall-less  TEPCs  in  the 
field  of  450  MeV/nucleon  (p  =  0.738)  40Ar  and  400  MeV/nucleon  (p  =  0.715)  12C.  One 

detector  had  a  diameter  of  2.54  cm  and  was  92  %  transparent  while  the  other  had  a  diameter 
of  0.64  cm  and  was  88  %  transparent.  Both  detectors  were  enclosed  in  a  10. 15  cm 
diameter  container.  The  ratio  of  container  to  cavity  diameters  for  the  two  detectors  were  4 
and  16,  respectively.  For  the  experiments,  0.25, 1,  and  4  ptm  diameter  sites  were 
simulated  in  both  detectors.  For  the  40Ar  ions  evaluated  at  the  plateau,  the  ratio  of  the  dose 
mean  lineal  energy  of  the  physically  larger  diameter  cavity  to  smaller  was  1.45  when  both 
detectors  were  simulating  1  pim.  At  the  plateau,  particle  kinetic  energy  was  429 

MeV/nucleon  (P  =  0.73),  L  was  about  95  keV//<m,  and  maximum  6-ray  energy  was  1.15 

MeV.  Under  similar  conditions,  measurements  made  with  10.5  cm  of  water  absorber 
yielded  a  ratio  in  dose  mean  lineal  energy  between  the  detectors  of  0.88.  Under  these 

conditions,  L  was  about  210  keV/^m,  p  =  0.46,  the  maximum  6-ray  energy  was  0.27 

MeV,  and  the  wall  effect  was  considerably  less.  The  data  collected  with  absorber  material 
is  difficult  to  interpret  because  the  particle  fluence  had  a  significant  contribution  from 
fragment  products  and  had  considerable  variance  in  particle  velocity. 
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Contributions  from  this  Work 


This  study  of  the  wall  effects  from  heavy  ions  extends  the  work  done  previously.  An 
original  and  novel  method  of  evaluating  the  wall  effects  of  heavy  ions  in  walled  TEPCs 
was  developed.  The  method  involved  the  use  of  a  particle  spectrometer  system  that  tagged 
every  particle  with  identification  of  the  charge  and  trajectory  through  the  experimental  set¬ 
up.  Having  knowledge  of  particle  trajectory,  the  analysis  could  be  tailored  to  evaluation  of 
the  response  of  the  detector  to  a  spatially  uniform  particle  fluence.  In  previous  work,  beam 
spreaders  (filters)  were  used  to  provide  a  more  spatially  uniform  particle  fluence  than  that 
of  an  unfiltered  beam.  Because  the  earlier  work  used  beam  spreaders,  assumptions  were 
made  concerning  the  spatial  uniformity  of  the  beam.  Furthermore,  filters  provided 
additional  fragmentation  and  increased  the  variance  in  velocity  of  the  primary  beam.  Both 
of  these  characteristics  are  undesirable  because  they  diminish  the  validity  of  the  assumption 
of  a  homogenous  beam  of  particles.  Thus,  conclusions  from  these  works  were  concerning 
qualitative  in  nature  because  the  beams  were  not  homogenous.  The  spectrometer  system 
use  in  this  work  allowed  evaluation  of  particles  that  had  similar  velocity  and  the  same 
charge  by  selection  on  the  energy  deposition  in  a  series  of  silicon  detectors. 

This  work  evaluated  the  wall  effects  for  1  GeV/nucleon  Fe  particles,  higher  energy  than 
previous  work.  The  highest  energy/nucleon  ions  used  in  earlier  work  was  600 
MeV/nucleon.  The  energy  spectra  of  galactic  cosmic  HZE  particles  extends  from  about  0. 1 
to  10  GeV/nucleon,  with  the  majority  of  the  fluence  between  0. 1  and  1  GeV/nucleon. 

Thus,  the  wall  effects  characterized  in  this  work  represents  a  practical  upper  bound,  as  the 

distortions  created  by  the  6-ray  effect  increase  for  increases  in  ion  velocity. 

The  primary  focus  of  this  work  was  the: 

1.  characterization  of  the  6-ray  effect  to  a  monoenergetic  beam  of  Fe  particles  based  on 
particle  trajectory  through  the  detector. 
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2.  calculation  of  Q  for  Fe  particles  based  on  measured  energy  deposition  in  a  TEPC 
and  under  assumptions  that  y  =  L  and  8/9y  =  L,  and 

3.  comparison  of  known  Q  factors  for  Fe  to  those  calculated. 

This  work  was  possible  because  a  particle  spectrometer  system  was  used.  These  details 
were  absent  from  earlier  work  on  heavy  ions. 

In  the  work  by  Rodgers  etal.  (1973),  the  ICRU  (1983)  speculated  on  the  cause  of  the 
very  high  energy  deposition  events  observed  in  the  single  event  spectra  collected  by  a 
walled  TEPC.  Because  a  particle  spectrometer  system  was  used  in  this  work,  evidence  as 
to  the  mechanism  responsible  for  these  events  could  be  provided. 

Also,  in  addition  to  analysis  of  detector  response  to  Fe  ions,  preliminary  evaluations  are 
provided  for  detector  response  to: 

1.  ions  with  Z  lower  than  Fe  (with  the  same  velocity  as  Fe  particles)  and 

2.  events  where  the  incident  Fe  particle  is  involved  in  a  charge-reducing  nuclear 
interaction  in  the  wall  of  the  detector. 

Both  of  these  evaluations  are  considered  preliminary  because  the  number  of  events 
analyzed  was  limited.  More  complete  evaluations  would  require  significantly  longer  beam 
time.  For  these  experiments,  a  6.35  mm  thick  copper  target  was  used.  Increased 
fragmentation  could  be  produced  through  use  of  thicker  targets,  but  increased  variance  in 
the  velocity  of  primary  and  secondary  particles  will  be  introduced. 
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Chapter  3:  Materials  and  Methods 

Beam-Line  Experimental  Set-Up  and  Equipment  Description 

Experiments  were  conducted  at  the  Alternating  Gradient  Synchrotron  (AGS)  at 
Brookhaven  National  Laboratory  (BNL)  using  S6Fe  ions  with  an  on-target  energy  of 
1.053  +  0.005  GeV/nucleon  (6  =  0.883).  Figure  3-1  shows  a  schematic  diagram  of  the 
experimental  arrangement .  The  TEPC  was  about  30  cm  from  each  PSD  pair. 

Figure  3-1.  Experimental  Arrangement  (Not  to  Scale). 


Triggers  PSD1  PSD2 


The  two  trigger  detectors,  Tj  and  T2,  were  330  ym  thick  silicon  detectors  with  active  areas 
of  300  mm2.  The  output  signal  from  each  was  sent  to  a  locally-mounted  charge  -sensitive 
preamplifier,  a  variable-gain  shaping  amplifier,  and  a  peak- voltage  sensing  analog-to-digital 
converter  (ADC).  Coincidence  between  T,  and  T2  defined  the  trigger  and  gates  for  the 
experiment  A  pair  of  position-sensitive  detectors  (PSD1)  was  placed  between  the  triggers 
and  TEPC,  while  another  pair  (PSD2)  was  placed  downstream  of  the  TEPC.  Each  pair 
was  oriented  to  provide  horizontal  and  vertical  information,  designated  as  (xt,  yt)  and  (x2, 
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y2)  in  Fig.  3-1.  Each  position-sensitive  detector  also  provided  a  signal  that  was 
proportional  to  the  total  energy  deposited  (i.e.,  AE  yj,  AE  xx,  .  . .).  In  context,  AE  is  a 

measure  of  stopping  power,  dE/dx,  without  normalizing  the  signal  by  the  thickness  of  the 
silicon  detector.  A  complete  description  of  the  PSDs  is  provided  by  Wong  etal .  (1990). 
One  experiment  was  conducted  with  a  6.35  mm  copper  target  while  all  others  were  without 
a  target 

Tissue-Equivalent  Proportional  Counter  Spherical  proportional  counters1  with 
internal  cavity  diameters  of  12.7  mm  (0.5  in)  were  used.  Figure  3-2  shows  cross-sectional 
diagrams  of  the  thin-walled  detector,  while  Figure  3-3  provides  cross-sectional  diagrams  of 
the  thick-walled  one.  The  left  diagram  of  both  figures  is  a  cross-sectional  view  on  the  same 
plane  as  the  experimental  set-up  given  in  Fig.  3-1.  The  right  diagram  is  a  cross-sectional 
view  along  the  beam  axis.  The  A-150  plastic  wall  had  a  thicknesses  of  1.27  and  2.54  mm, 
respectively  for  the  two  detectors.  A  stainless  steel  central  anode  wire  was  coaxially 
surrounded  by  a  stainless  steel  helix  1  mm  in  diameter.  The  potential  between  the  helix  and 
anode  was  used  to  create  a  symmetrical  electric  field  along  the  length  of  the  anode  (ICRU 
1983).  A  hole  in  the  wall  was  introduced  for  gas  filling.  Since  this  hole  was  located  on  the 
upstream  side  of  the  thin-walled  detector,  particles  intersecting  this  region  of  the  detector 
were  rejected  from  analysis.  The  vent  hole  for  the  thick-walled  detector  was  on  the 
downstream  side  of  the  detector.  Each  detector  was  encapsulated  within  a  0.018  mm  shell 
of  aluminum  that  maintained  the  desired  vacuum.  Figure  3-4  shows  “pathlengths”  through 
the  cavity  and  wall  with  respect  to  impact  parameter  (distance  from  the  detector  center)  for 
the  thin  and  thick-walled  detectors. 

The  proportional  counters  were  filled  with  propane-based  tissue  equivalent  gas 
consisting  of  C02,  N2,  propane,  and  a  trace  amount  of  isobutane  in  mass  fractions  shown 
in  Table  C-l.  The  ideal  gas  law: 


1  Far-West  Technologies,  Goleta  CA 
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Figure  3-2.  Cross-Sectional  Diagrams  of  Thin-Walled  TEPC 
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Cross-Sectional  View  Cross-Sectional  View 

Perpendicular  to  Beam  Axis  Along  Beam  Axis 


Figure  3-3.  Cross-Sectional  Diagrams  of  Thick-Walled  TEPC 
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Cross-Sectional  View  Cross-Sectional  View 

Perpendicular  to  Beam  Axis  Along  Beam  Axis 


P  V  =  n  R, 


(3-1) 


was  used  to  determine  pressure,  P,  of  the  gas  required  to  simulate  the  required  diameters 
where  V  is  detector  volume,  n  is  number  of  moles  of  the  fill  gas,  and  R  =  24.1  liter  atm 
mole'1  at  21°  C. 


0  3  6  9 


Impact  Parameter  (mm) 

Table  C-l  summarizes  the  calculated  parameters  for  simulation  of  a  3  pi n  diameter  site  at 
21°C  and  pressure  of  13,332  Pa  (100  mmHg).  Other  diameters  were  simulated  by  directly 
scaling  the  pressure.  For  each  use  of  the  detector,  a  minimum  of  three  flushings  of  the 
detector  were  performed  prior  to  use. 

The  ratio  of  densities  between  the  wall  and  fill  gas  was  about  14,400,  9,600,  and  4800, 
respectively  for  simulated  unit  density  diameters  of  1, 2,  and  3  pm. 
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Table  3-1.  Calculated  Parameters  for  Propane-Based  Tissue-Equivalent  Gas  at 
21  °C,  Volume  =  1.0725  cm3.  Pressure  =  13,332  Pa,  and  Simulated  Diameter  =  3  /*m. 


Molecule 

Mass  Fraction 

Moles 

Molecular  Weight 
(grams) 

Mass  (grams) 

co2 

0.3939 

2.30  x  10-6 

44.01 

1.013  x  10-4 

n2 

0.0564 

3.30  x  lO'7 

28.02 

9.237  x  10‘6 

c3h8 

0.5485 

3.21  x  lO"6 

44.09 

1.414  x  10-4 

c4h10 

0.0012 

7.01  x  10'9 

58.12 

4.077  x  10'7 

Total 

1 

5.847  x  10'6 

— 

2.523  x  10-4 

Detector  bias  potentials  were  selected  to  provide  reasonably  high  gas  amplification  with 
acceptable  proportionality.  A  summary  of  bias  potentials  used  for  each  detector  and 
simulated  site  diameter  is  listed  in  Table  3-2. 


Table  3-2.  TEPC  Bias  Potentials  (Volts). 


Detector 

1  yira 

2  pim 

3  pi n 

Thin  (SN  1373) 

600 

618 

658 

Thick  (SN  1376) 

598 

606 

644 

Calibrations  were  performed  using  the  most  probable  event  size  value  (unscattered  a- 

particles)  from  a  built-in  gravity  controlled  244Cm  a-radiation  source.  Using  published 

values  of  stopping  power  values  (ICRU  1993)  and  applying  a  continuous  slowing  down 
approximation  to  account  for  energy  loss  in  the  detector,  the  mean  energy  transferred  from 

unscattered  a-particles  crossing  the  cavity  diameter  was  84. 1  keV//rm  for  the  thin- walled 
detector  at  1  pim.  Calibrations  were  also  performed  using  the  proton  edge  from  an 
unmoderated  239PuBe  neutron  source  and  a-edge  from  the  built-in  a-source  by  the 
technique  described  by  Schrewe  etal.  (1993).  When  empirical  functions  of  Goodman  and 
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Coyne  (1980)  for  mean  energy  per  ion  pair  formed  (W)  were  applied,  all  three  methods 
agreed  to  within  +  2.4  %.  Details  of  the  calibrations  are  in  Appendix  B. 

The  output  signal  from  the  detector  was  sent  to  a  locally-mounted  charge-sensitive 
preamplifier,  split,  and  sent  through  two  identical  variable-gain  shaping  amplifiers  with  a 
ratio  of  gains  equal  to  16. 

Data  Analysis 

Particle  Identification 

The  first  step  in  the  analysis  consisted  of  evaluating  energy  deposited  in  the  sequence  of 
trigger  detectors  and  positional  devices.  Since  data  were  recorded  on  an  event  by  event 
basis,  the  silicon  devices  could  be  used  to  determine  particle  species  along  the  experimental 
beamline. 

Signals  in  the  trigger  detectors  were  normally  distributed  about  a  mean  value 
corresponding  to  the  energy  loss,  AE,  of  Fe.  Events  with  AE  beyond  +  3a  of  the  mean  in 

either  T,  or  T2  were  rejected  from  the  analysis  to  insure  particles  entering  the  target  or 
PSD1  were  iron. 

Similarly,  for  the  experiments  conducted  without  a  target,  an  event  was  rejected  if  AE 

in  PSD1  and  PSD2  was  beyond  +  2.5a  from  the  mean  corresponding  to  Fe.  For  the 
experiment  conducted  with  the  Cu  target,  individual  particles  species  could  be  identified 
and  selected  for  analysis  based  on  the  AE  signal  in  the  PSDs.  The  relative  location  of  the 

peaks  in  the  PSD  AE  distributions  were  roughly  proportional  to  the  Z2  of  the  particle 

responsible  for  each  peak. 

Particle  Trajectories 

An  analysis  of  y,  vs.  y2  and  x2  vs.  x2  indicated  that  the  vertical  and  horizontal  position 
of  individual  particles  through  the  two  PSD  pairs  had  good  correlation  indicating  that  the 
beam  was  highly  parallel.  Events  without  good  correlation  were  believed  to  be  due  to 
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nuclear  collisions  at  some  location  along  the  beam  axis.  Because  the  purpose  of  the 
experiments  conducted  without  the  Cu  target  were  to  characterize  the  response  of  the  TEPC 
to  Fe  particles  not  involved  in  nuclear  interactions,  events  without  good  correlation  were 
rejected  from  analysis.  The  experiment  with  the  Cu  target  was  not  screened  for  trajectory 
correlation. 

The  exact  location  of  the  TEPC  in  reference  to  the  PSD’s  was  accomplished  using  an 
iterative  process  that  imaged  the  shape  of  the  detector  in  each  PSD  pair.  This  was  done  by 
selecting  either  very  small  energy  deposition  events  (i.e.,  particles  that  passed  outside  of 
the  wall)  or  very  high  energy  deposition  events  (i.e.,  particles  that  passed  along  the  cavity/ 
wall  interface).  The  final  uncertainty  in  the  horizontal  and  vertical  location  of  the  detector 
center  was  estimated  to  be  less  than  0. 1  mm. 

After  the  geometry  was  established,  the  trajectory  of  each  particle  that  satisfied  the  AE 

selection  criteria  was  reconstructed  under  the  assumption  that  the  particle  had  a  straight 
trajectory  through  the  system.  The  impact  parameter  (i.e.  perpendicular  distance  from  the 
center  of  the  detector)  and  chord  length  through  the  detector  were  computed  and  placed  in  a 
data  file  with  other  information  associated  with  that  event.  Details  of  these  calculations  are 
provided  in  Appendix  D.  Lineal  energy  was  computed  for  the  two  TEPC  response  signals 

by  dividing  the  energy  deposited,  e,  by  the  mean  chord  length  of  the  simulated  sphere.  For 

simulation  of  a  1  pa n  diameter  sphere,  the  mean  chord  length  for  /(-random  particle 
traversals  was  %  pim.  The  lineal  energy  from  the  two  amplified  TEPC  signals  overlapped 
at  about  20  keV  pan1.  The  two  signals  were  combined  into  one  in  for  data  analysis. 

TEPC  Response  vs  Impact  Parameter  for  Experiments  without  Cu  Target 

Detector  response  was  evaluated  in  terms  of  lineal  energy  as  a  function  of  the  impact 
parameter  of  the  trajectory.  This  was  done  using  two  methodologies.  The  first  method 
sorted  events  into  sequential  intervals  of  impact  parameter  from  0  to  9.6  mm.  There  were 
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about  1000  events  in  each  interval  and  the  mean  value  of  energy  deposition  associated  with 
that  impact  parameter  was  determined  by  fitting  the  data  to  a  normal  distribution. 

Uncertainties  in  positional  information  can  have  a  large  effect  on  observed  energy 
deposition  for  impact  parameters  near  the  cavity/wall  interface  where  the  TEPC  response 
varied  rapidly  with  small  changes  in  impact  parameter.  For  these  regions  of  the  detector, 
the  selection  process  was  reversed.  Events  were  sorted  into  sequential  intervals  of  energy 
deposition  and  an  analysis  of  the  impact  parameter  was  performed  each  interval.  Both 
methods  converged  to  the  same  result  for  intermediate  regions. 

Experiment  with  CU  Target 

The  data  from  this  experiment  was  split  into  three  separate  data  files  based  on  the 
particle  identified  in  PSD1  and  PSD2. 

1.  Fe  in  PSD1  and  PSD2  Events  in  this  data  file  had  a  AE  signal  in  all  four 
PSD  devices  that  was  within  ±  2.5 o  of  the  mean  corresponding  to  Fe. 

2.  Fe  in  PSD1.  but  Fragment  in  PSD2  Events  in  this  data  file  had  a  AE 

signal  in  both  PSD1  detectors,  but  had  a  AE  signal  at  least  2.5a  less  that  the  mean 

corresponding  to  Fe  in  both  PSD2  detectors.  These  events  were  believed  to  be  involved  in 
charge-changing  nuclear  interactions  somewhere  between  the  two  PSD  pairs,  with  the 
likeliest  location  being  in  the  wall  of  the  TEPC.  Because  the  Fe  particles  were  involved  in  a 
nuclear  interaction,  the  assumption  that  particle  had  a  straight  trajectory  between  the  two 
PSD  pairs  is  invalid.  The  impact  parameter  for  these  particles  was  based  only  on  the 
position  data  in  PSD1  and  will  have  greater  uncertainties  than  determination  of  impact 
parameter  using  both  pair  of  PSDs.  Appendix  D  provides  estimates  of  the  uncertainties  in 
impact  parameter.  Since  the  response  of  unfragmented  Fe  particles  (above)  are  compared 
to  this  data,  impact  parameter  was  calculated  for  the  Fe  events  in  this  manner  as  well  so 
bias  is  not  introduced. 
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3.  Fragment  in  PSD1  and  PSD2  Events  in  this  data  file  had  a  AE  signal  at 


least  2.5a  less  than  the  mean  corresponding  to  Fe  in  all  four  PSD  detectors.  These  events 

were  believed  to  be  involved  in  charge-changing  nuclear  interactions  between  T2  and  PSD 
y15  with  the  likeliest  location  being  in  the  copper  target.  Due  to  the  low  number  of 
fragments  produced  in  the  copper  target,  the  analysis  was  restricted  to  evaluation  of  energy 
deposition  of  fragments  with  impact  parameter  less  than  3  mm. 

Model  for  Energy  Deposition 

The  Chatteijee  model  was  used  to  calculate  energy  deposition  in  a  1  ptm  diameter  sphere 
for  a  range  of  impact  parameters  using  an  LET  of  149  keV//un  in  propane-based  TE  gas 
(ICRU  1993).  Details  of  the  calculations  are  provided  in  Appendix  A. 
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Chapter  4:  Results 


Energy  Deposition  Measurements  in  the  Position  Sensitive  Devices 

Figure  4-1  is  a  plot  of  energy  deposition  in  PSD  y2  for  the  experiment  with  the  thin- 
walled  detector  and  copper  target  in  place.  The  prominent  peak  near  270  MeV  corresponds 

to  Fe  particles  and  the  peaks  with  lower  values  of  AE  correspond  to  fragments.  The  plot  is 

truncated  to  display  the  peaks  for  energy  deposition  of  fragments  of  Z  =  17  (Cl)  to 
unfragmented  Fe.  For  charge-changing  interactions  at  this  velocity,  the  fragments  have 

velocities  close  to  the  velocity  of  unfragmented  Fe.  Accordingly,  AE  values  for  individual 

peaks  scaled  closely  to  the  Z2  of  the  ion  as  expected  from  the  Bethe-Bloch  stopping  power 
formula  (Attix  1986). 

Normal  distributions  were  fit  to  the  Fe  peaks  in  the  energy  spectra  of  the  four  PSDs. 
The  coefficients  of  variation  (CV)  (a/p)  were  less  than  2  %.  Discrimination  of  Fe  events 

from  fragments  using  paired  PSD  AE  signals  (i.e.  AE  yt  &  AE  xj  was  accomplished  with 

less  than  1  %  probability  of  falsely  selecting  a  fragment,  based  on  the  mean  energy 
deposition  for  Fe,  Mn,  and  the  CV  of  the  data  in  each  peak. 

Figs.  4-2a  and  b  are  AE  distributions  for  PSDs  xt  and  y,  after  Fe  was  selected  in  T,, 

T2,  and  PSD1.  The  CVs  for  the  two  distributions  were  1.575  +  0.004  %  and  1.454  + 
0.004,  respectively,  based  on  normal  distribution  fits  to  the  data.  Variance  in  the 
distributions  is  caused  by  many  factors  including:  1)  energy  straggling,  2)  variance  in 
energy  deposition  due  to  non-uniformity’s  in  detector  thickness,  3)  random  errors  in  signal 
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Number  of  Events 


Figure  4-1.  AE  in  PSD  Y2  for  Experiment 
with  Copper  Target  (Labels  Correspond 
to  the  AE  of  Individual  Ions) 


AE  y2  (MeV) 
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Figure  4-2.  Energy  Deposition  Distributions  of  Fe  Peaks  in  PSD  XI  and  Y1 
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processing,  4)  variance  in  the  stopping  power  of  the  particles,  and  5)  variance  in  energy 
deposition  due  to  pathlength  variations  through  the  PSDs. 

All  of  the  factors  are  stochastic  in  nature  and  independently  affect  the  energy  deposition 
in  a  single  PSD.  Additionally,  the  first  three  factors  independently  affect  energy  deposition 
between  each  detector  in  a  PSD  pair.  However,  the  last  two  factors  are  covariant  and  will 
affect  the  energy  deposition  in  each  detector  in  a  PSD  pair  in  the  same  manner. 

Per  discussions  with  Dr  Zeitlin  (1997),  the  CV  introduced  by  each  of  the  first  two 
factors  was  about  1  %  and  jointly  account  for  the  majority  of  the  variance  in  the 
distributions  of  Fig.  4-2.  To  assess  the  relative  contribution  of  the  last  factor  on  the 
variance  in  energy  deposition,  an  evaluation  of  particles  trajectory  was  performed.  From 
distributions  of  x2  -  x}  and  y2  -  yt  for  Fe  particles  that  traverse  through  PSD1  and  2  without 
charge-reduction  fragmentation,  distributions  of  pathlength  through  the  PSDs  could  be 
calculated.  The  CV  of  pathlength  through  the  PSDs  was  estimated  at  8  x  10‘4  %.  This 
value  is  very  small  compared  to  the  total  CV  of  energy  deposition  in  the  PSDs. 

A  distribution  of  for  the  events  of  Fig.  4-2  is  plotted  in  Fig.  4-3.  The  data 

AExj 

agrees  well  to  a  normal  distribution.  The  CV  for  this  distribution  is  1.936  +  0.004  %. 
Appendix  E  contains  an  analysis  of  variance  of  the  data  given  here.  The  estimated 

covariance  between  AE  y2  and  AE  x1  was  3. 1  MeV2  and  represents  about  22  %  and  16  % 

of  the  total  variance  in  AE  yt  and  AE  x,,  respectively.  All  factors  contributing  to  the 

covariance  between  the  energy  deposition  in  the  two  detectors  is  not  known.  Variance  in 
stopping  power  of  the  Fe  particles  is  one  factor.  Without  knowledge  of  the  contribution  of 
all  factors,  however,  it  can  be  concluded  that  the  covariance  introduced  by  the  variance  in 
stopping  power  is  limited  by  the  total  covariance.  Therefore,  because  the  covariance 
between  the  two  distributions  is  low,  it  can  be  concluded  that  the  Fe  particles  used  for 
analysis  were  very  homogenous  in  stopping  power. 
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Figure  4-3.  Distribution  of  (AE  yl  /  AE  xl) 
for  Fe  Particles  in  PSD1 
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Positioning  of  the  TEPC 


The  location  of  the  TEPC  in  the  z-dimension  was  known  based  on  measured  distances 
between  the  PSDs  and  TEPC.  However,  fine  positioning  of  the  device  in  the  x-y  plane 
was  based  on  evaluation  of  the  TEPC  response.  Fig.  4-4  is  a  plot  of  detector  response  for 
Fe  particles  in  the  thin- walled  TEPC.  A  1  ptm  diameter  site  was  simulated  without  the 
copper  target . 

Fig.  4.5a  is  a  scatter  plot  of  Xj  vs.  y,  for  events  that  have  lineal  energy  distinctly  greater 
than  those  in  the  peak  of  the  distribution  (y  >  210  keV  /mi'1).  About  one-third  of  these 
events  are  located  in  proximity  to  the  anode  and  helix  wires  and  are  assumed  to  have  hit  a 

wire.  The  enhanced  energy  deposition  is  likely  due  to  energetic  6-rays  produced  in  the 

wire.  The  other  events  are  located  in  proximity  to  the  cavity/wall  interface. 

By  selection  of  events  with  lineal  energy  much  higher  than  210  keV  ptm'1,  wire  hits 
events  could  be  eliminated  .  The  distribution  of  events  with  y  >  255  keV  ptm'1  were  located 
primarily  in  proximity  to  the  cavity/wall  interface.  For  position  determinations,  these 
events  were  assumed  to  arise  from  trajectories  at  the  cavity/wall  interface  (impact  parameter 
=  6.35  mm).  After  many  iterations  of  estimating  scaling  and  position  offset  factors  for 
each  PSD  pair,  coordinates  of  a  line  through  the  TEPC  center  and  the  PSD  pairs  was  made. 
Fig.  4-6  contains  a  histogram  of  impact  parameter  for  events  with  y  >  255  keV  ptm'1.  For 
events  in  the  peak  of  the  distribution,  the  mean  was  estimated  at  6.356  +  0.019  mm.  The 
percent  fractional  error  of  the  this  value  with  the  specified  detector  radius  of  6.35  mm  was 
(-  0.2,  +  0.4  %). 

Figure  4-5b  is  a  scatter  plot  of  x  vs.  y  in  the  detector  plane  for  events  with  lineal  energy 
between  25  and  35  keV  ptm'1.  The  scatter  plot  was  generated  for  corrected  position 
parameters  and  is  reflected  by  (0, 0)  coordinates  of  the  detector  center  and  circular 
symmetry  of  the  particles  with  respect  to  center.  The  deviation  of  the  particles  from  a 
perfect  circular  pattern  arises  from  uncertainties  in  the  positional  devices  and  variance  in  the 
actual  location  of  interaction.  Particles  distant  from  the  circular  pattern  were  assumed  to  be 
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Figure  4-4.  Detector  Response  to  Fe  Particles 
with  a  Simulated  Diameter  of  1  jim  (No  Target) 
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Figure  4-5.  Scatter  Plots  in  Position  for  Detector  Response  to  Fe  Particles  in 
Two  Selected  Intervals  of  Lineal  Energy  [Simulated  Diameter  =  1  um. 
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the  result  of  nuclear  scattering.  The  mean  impact  parameter  of  unscattered  particles  was 
about  6.6  mm  and  corresponds  to  trajectories  though  the  wall  but  not  the  cavity.  Nuclear 
scattering  will  be  discussed  later  in  detail. 

Impact  Parameter  Uncertainties  and  Bias 

Uncertainties  in  impact  parameter  were  estimated  based  on  the  distribution  of  impact 
parameter  for  events  with  y  >  255  keV  /mi'1.  At  the  cavity/wall  interface,  the  estimated 
uncertainty  in  impact  parameter  was  assumed  to  be  approximately  equal  to  the  standard 
deviation  of  the  events  contained  in  the  peak.  For  events  between  5.7  and  7.0  mm  in  Fig. 
4-6,  the  standard  deviation  was  0.23  mm.  Uncertainties  in  impact  parameter  were 
estimated  to  be  close  to  0.23  mm  for  the  entire  range  of  impact  parameter,  based  on  the 
assumption  that  uncertainty  in  impact  parameter  was  caused  by  the  uncertainty  in  position 
signals  from  the  PSDs.  Details  of  the  evaluation  are  in  Appendix  D. 

Uncertainties  in  position  signals  from  the  PSDs  will  propagate  bias  in  calculated  values 
of  impact  parameter.  Estimations  of  bias  in  impact  parameter  were  made  with  Monte  Carlo 
simulations  detailed  in  Appendix  D.  The  bias  was  positive  and  had  the  greatest  values  for 
small  impact  parameters.  The  estimated  bias  was  less  than  2.6  %  for  impact  parameters 
greater  than  1  mm.  The  bias  was  determined  to  have  a  negligible  effect  on  this  work 
because:  1)  the  interval  of  impact  parameters  that  contained  significant  bias  in  impact 
parameter  reflected  a  small  fraction  of  the  total  fluence  and  2)  the  detector  response  in  the 
interval  of  impact  parameter  having  the  greatest  bias  had  the  most  uniform  detector 
response.  A  summary  of  estimated  of  bias  values  are  in  Appendix  D. 

Detector  Response  vs.  Impact  Parameter  for  Fe 

Fig.  4-7  shows  the  distribution  of  lineal  energy  for  Fe  particles  in  the  thin-walled 
detector  for  several  selected  intervals  of  impact  parameter.  Collectively,  these  distributions 
were  used  to  construct  the  detector  response  vs.  impact  parameter. 

Figure  4-7a  Fig.  4-7a  corresponds  to  trajectories  through  the  center  of  the  detector. 
The  solid  line  represents  a  fit  to  the  peak  using  a  normal  distribution.  The  CV  for  the 
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Figure  4-6.  Distribution  of  Impact  Parameter 
for  Fe  Events  with  High  Lineal  Energy 
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Figure  4-7.  Detector  Response  to  Fe  for  Selected  Impact 
Parameter  Intervals  [Simulated  Diameter  =  1  urn] 
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Figure  4-7  cont.  Detector  Response  to  Fe  for  Selected  Impact 
Parameter  Intervals  [Simulated  Diameter  =  1  um] 
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distribution  was  about  8  %.  An  analysis  of  variance  indicated  that  this  was  primarily  due  to 
energy  straggling  and  variations  in  gas  amplification.  There  were  a  few  events  that 
occurred  beyond  the  tails  of  the  normal  distribution.  The  higher  values  were  attributed  to 
trajectories  through  the  anode  and  helix  wires  (wire  hits).  Very  small  values  of  lineal 
energy  were  attributed  to  events  with  nuclear  interactions  in  the  walls  of  the  detector. 

These  particles  were  believed  to  have  impact  parameters  near  the  cavity/wall  interface. 
Scatter  of  the  incident  particle  caused  the  calculation  of  impact  parameter  to  be  incorrect. 

Trajectory  Analysis 

All  of  the  distributions  of  lineal  energy  had  events  of  lineal  energy  below  the  tail  of  the 
normal  distribution.  To  determine  whether  these  events  were  the  result  of  nuclear 
scattering  processes,  distributions  of  (x2  -  x,)  were  created  using  events  with  impact 
parameter  <  3  mm.  Events  with  impact  parameter  in  this  range  were  selected  because  it 
was  broad  enough  to  provide  a  sufficient  number  of  particles  for  analysis,  but  sufficiently 
distant  from  the  wall  such  that  events  with  incorrect  impact  parameter  would  be  due 
primarily  to  nuclear  scattering  rather  than  that  caused  by  PSD  uncertainties  alone. 

Fig  4-8a  contains  a  plot  of  lineal  energy  for  Fe  events  in  the  thin-walled  detector 
simulating  1  pm.  The  distribution  of  (x2  -  xt)  for  a  selection  of  particles  with  lineal  energy 
about  the  mean  of  the  peak  is  shown  in  Fig.  4-8b.  The  mean  of  this  distribution  was  -3.6 
mm  with  a  standard  deviation  of  0.77  mm.  These  events  were  assumed  to  be  unscattered 
through  die  TEPC.  Similar  distributions  were  collected  for  events  with  lineal  energy  in  the 
intervals:  0  -  50, 50  - 100,  >  210  keV  pm1,  and  repeated  for  the  thick-walled  detector.  A 
summary  is  provided  in  Table  4-1.  From  the  table,  the  mean  and  standard  deviations  of  the 
events  in  the  intervals:  0-50  and  50  - 100  keV  pm'1  are  significantly  different  than  that  for 
events  in  the  peak  region  of  the  distribution.  In  particular,  the  standard  deviations  of  these 
distributions  are  considerably  higher.  As  these  particles  are  otherwise  identical  to  the 
particles  assumed  to  be  unscattered,  it  is  logical  to  conclude  that  the  increase  in  standard 
deviation  is  due  to  scattering.  The  events  with  y  >  210  keV  pm1  have  distributions  of 
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Figure  4-8.  Detector  Response,  and  (X2  -  XI)  for  Particles  in 
Peak  Region,  for  Fe  Events  with  Impact  Parameters  3  mm 
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(x2  -  Xj)  with  mean  and  standard  deviation  similar  to  that  of  the  events  in  the  peak.  This 
indicates  that  it  is  unlikely  that  these  events  were  involved  in  nuclear  scattering. 


Table  4-1.  (X2  -  Xt)  Summary  Data  for  Fe  Particles  with  Impact  Parameters 
<  3  mm  in  the  Thin-  and  Thick-Walled  Detectors  [Simulated  Diameter  =  1  jxm] 


Lineal  Energy 
Interval  (keV/j<m) 

Parameter 

Detector 

Thin-Walled 

Thick-Walled 

All 

n 

7273 

6813 

Peak 

n 

3653 

3174 

(Xj  -  X,)  }i  (mm) 

-3.6 

-4.0 

(x2  -  Xj)  a  (mm) 

0.77 

0.78 

0-50 

n 

88 

106 

(x2  -  Xj)  J4  (mm) 

-3.1 

-3.6 

(x2  -  Xj)  a  (mm) 

1.2 

1.2 

50-  100 

n 

18 

12 

(x2  -  X!)  n  (mm) 

-3.2 

-3.6 

(x2  -  Xi)  a  (mm) 

0.92 

1.3 

>210 

n 

62 

148 

(*2  -  x,)  \i  (mm) 

-3.5 

-3.9 

(x2  -  xj  a  (mm) 

0.83 

0.76 

Figure  4-7b  Fig.  4-7b  shows  the  data  for  impact  parameters  between  4.0  and  4.2 
mm.  The  CV  for  the  fitted  distribution  was  12  %.  The  impact  parameter  interval  is  for  a 
region  of  the  detector  where  small  changes  of  impact  parameter  introduce  larger  changes  of 
pathlength  through  the  cavity  due  to  the  curvature  of  the  spherical  wall.  Thus,  uncertainties 
in  the  measured  position  of  the  particle  are  associated  with  a  larger  variance  in  lineal  energy 
compared  with  particles  passing  through  the  center  of  the  detector. 

Figure  4-7c  and  d  Fig.  4-7c  shows  the  data  for  impact  parameters  between  6.25 
and  6.35  mm.  These  should  be  trajectories  just  inside  the  cavity/wall  interface.  However, 
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because  of  uncertainties  in  position  measurements  and  nuclear  scattering,  the  data  also 
contains  particles  that  do  not  traverse  the  cavity.  The  central  region  of  the  abscissa 
(70  <  y  <  170  keV  pim'1)  corresponds  to  trajectories  through  the  cavity.  The  region  of 
small  y  (<  70  keV  //m'1)  corresponds  to  trajectories  through  the  wall-only  and  the  large 
y  (>  170  keV  pim1)  are  from  trajectories  that  graze  the  cavity/wall  interface. 

Fig.  4-7d  shows  the  data  for  impact  parameters  between  6.35  and  6.45  mm.  These 
should  be  trajectories  just  outside  the  cavity/wall  interface.  The  distribution  is  similar  to 
that  of  Fig.  4-7c  for  large  y  events,  but  is  enhanced  in  low  y  events  and  suppressed  in 
medium  y  events. 

Figure  4-9  Figure  4-9a  shows  the  detector  response  for  impact  parameters  between 
7.95  and  8. 1  mm  for  Fe  particles  in  the  thick-walled  detector  simulating  1  ptm.  Distribut¬ 
ions  of  (x2  -  xx)  were  generated  separately  for  events  less  than  25  keV  pim'1  (small  y)  and 
those  greater  than  25  keV  ptm1.  Figure  4-9b  shows  the  distribution  for  events  greater  than 
25  keV  ja n1  with  mean  and  standard  deviations  of  -3.6  and  1.20  mm,  respectively.  The 
distribution  for  the  small  y  events  is  not  displayed,  but  had  a  mean  and  standard  deviation 
of  -4.2  and  0.85  mm,  respectively.  For  the  large  y  events,  the  standard  deviation  is  41  % 
higher  than  the  small  y  events  and  indicates  that  these  events  were  likely  involved  in  nuclear 
interactions.  It  is  reasonable  to  conclude  that  these  events  had  trajectories  through  the 
cavity,  but  had  incorrect  impact  parameter  due  to  nuclear  scattering.  Although,  the  large  y 
events  are  few  in  number,  the  mean  detector  response  will  be  biased  to  higher  values  for 
these  intervals  of  impact  parameter. 

Impact  Parameter  vs,  TEPC  Response  for  Fe 

Impact  parameter  uncertainties  and  incorrect  calculation  of  impact  parameter  due  to 
nuclear  scattering,  made  assessment  of  detector  response  difficult  by  the  method  used 
above  for  trajectories  near  the  cavity/wall  interface  and  in  the  wall-only.  An  alternate 
method  of  analysis  was  used  to  overcome  these  difficulties.  Distributions  of  impact 
parameter  were  created  based  on  selected  intervals  of  detector  response.  Figure  4-10 
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Figure  4-9.  Detector  Response,  and  (X2  -  XI)  for  Particles  with  y  >  25  keV/jim 
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y  (keV//im)  (x2  -  xl) 


Figure  4-10.  Distributions  of  Impact  Parameter  for  Selected  Intervals  of 
Detector  Response  in  the  Thin-Walled  Detector  [Simulated  Diameter  =  1  jim] 
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Impact  Parameter  (mm)  Impact  Parameter  (mm) 


provides  two  examples  of  impact  parameter  distributions  for  selections  of  detector  response 
in  the  thin-walled  detector  simulating  1  // m.  For  each  distribution,  the  solid  line  represents 
a  normal  distribution  fit  to  the  peak  of  the  data.  Based  on  the  mean  of  the  normal 
distributions,  these  events  were  for  trajectories  through  the  wall-only.  The  standard 
deviation  of  the  normal  distributions  were  0.26  and  0.32  mm,  respectively,  for  the  data  of 
Figs.  4-  10a  and  b.  Overall,  for  distributions  created  in  this  manner,  the  standard  deviations 
of  the  normal  distribution  fits  to  the  peak  were  at  a  minimum  at  the  cavity/wall  interface. 
This  trend  in  the  data  is  logical  since  the  incremental  change  in  detector  response  vs.  the 
change  in  impact  parameter  was  maximum  at  the  cavity/wall  interface.  This  trend  was 
consistent  for  the  two  detectors  regardless  of  the  simulated  diameter. 

Composite  Plots  of  TEPC  Response  vs.  Impact  Parameter  for  Fe  Events 

General 

Composite  plots  of  TEPC  response  vs.  impact  parameter  were  generated  by  combining 
the  data  from  both  of  the  methods  described  above.  Figure  4-11  shows  the  overlap  in  data 
for  in  the  intermediate  region  of  impact  parameter.  The  data  from  the  two  methods  overlap 
well  from  4  to  4.7  mm.  For  the  individual  experiments,  plots  of  the  overlap  region  were 
made.  The  interval  of  impact  parameter  with  the  best  overlap  in  the  two  data  sets  was 
selected  qualitatively  based  on  visual  inspection.  The  data  of  detector  response,  based  on 
selection  of  events  within  intervals  of  impact  parameter,  was  truncated  for  impact  parameter 
intervals  greater  than  the  overlap  region.  While,  the  data  of  impact  parameter,  based  on 
selection  of  events  within  intervals  of  detector  response,  was  truncated  for  impact 
parameters  less  than  the  overlap  region. 

1  Mm  Simulated  Diameter  Experiments  -  Fe  Only 

Figure  4-12  is  a  plot  of  mean  lineal  energy  vs.  impact  parameter  through  the  TEPC  for 
the  thin-walled  detector  to  Fe  particles.  The  open  circles  represent  the  data  with  a  smooth 
line  through  the  distribution  of  points.  The  vertical  bars  on  two  of  the  data  points  represent 
the  standard  deviation  of  the  data  at  these  points.  The  solid  line  represents  lineal  energy  in 
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Figure  4-11.  Overlap  in  Data  for  Both  Analysis  Methods  in  the  Intermediate  Region 
of  Impact  Parameter  for  the  Thick-Walled  Detector  for  Fe  Particles  [Simulated 
Diameter  =  3  pim.  Representative  Errors  Bars  Shown  for  Solid  Circles  (•),  Vertical 
Errors  Bars  for  Open  Circles  (O)  were  between  0.5  to  0.7  keV  pim1  (not  shown)] 
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a  sphere  of  tissue  that  is  1  pm  in  diameter  using  the  Chatterjee  model  that  includes  radial 
energy  loss  due  to  high  energy  6-rays.  The  solid  circles  correspond  to  the  distribution  of 
lineal  energy  for  particles  with  L  =  149  keV  //m'1  but  assuming  that  all  energy  transfer  is 
absorbed  locally.  Thus,  there  is  no  loss  of  energy  radially  and  e  =  L  •  pathlength  through 
the  cavity. 

For  trajectories  through  the  center  of  the  gas  cavity,  the  mean  lineal  energy  obtained 
from  the  model  is  70  %  of  the  value  for  particles  with  no  radial  energy  loss.  However,  the 
data  were  higher  than  predicted  by  the  model  that  only  incorporated  radial  energy  loss. 

This  implies  that  there  is  a  significant  reduction  in  energy  deposition  due  to  the  escape  of 

6-rays  and  that  this  is  only  partially  compensated  by  6-ray  production  in  the  front  wall  of 

the  detector. 
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Figure  4-12.  Data  from  Experiment  with  Fe  Particles  in  Thin- Walled  Detector 
Simulating  1  Jim  [Open  Cirles  (O)  is  Data,  Solid  Line  is  Chatterjee  Model,  and 
Solid  Circles  (•)  is  Particles  without  Radial  Loss  of  Energy] 
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For  trajectories  near  the  cavity/wall  interface,  the  data  are  higher  than  the  model  that 
does  not  include  radial  energy  loss.  Thus,  the  reduction  in  energy  deposition  in  the  cavity 

due  to  the  escape  of  6-rays  created  in  the  cavity  is  more  than  compensated  by  the  added 

energy  deposition  from  6-rays  created  in  the  wall. 

Figure  4-13  is  a  plot  of  mean  lineal  energy  vs.  impact  parameter  through  the  TEPC  for 
the  thick-walled  detector.  The  open  circles  represent  the  data  with  a  smooth  line  through 
the  distribution  of  points.  The  solid  line  is  the  distribution  of  lineal  energy  using  the 
Chatteijee  model  as  in  Fig.  4-12. 

1.  Moments  of  the  Distribution  of  Mean  Lineal  Energy.  The  data  in 
Fig.  4-13  were  used  to  determine  yF  and  yD  for  //-random  (ICRU  1983)  Fe  particles  at 
1  GeV/nucleon.  Table  4-2  provides  a  summary  of  the  results  compared  with  computations 
assuming  that  e  =  L  •  pathlength  through  the  cavity.  Calculations  were  performed 

according  to  the  method  in  Appendix  F.  A  marked  difference  exists  in  the  frequency  mean 
lineal  energy,  with  the  data  having  a  value  nearly  half  of  that  for  the  particles  with  no  radial 

energy  loss  from  6-rays.  The  dose  mean  lineal  energy  values  are  closer,  with  the  value  for 

the  data  being  15  %  lower  than  that  of  the  model  compared. 

Table  4-2.  Computed  yF  and  yD  Based  on  Measured  Response  of  Fe  Particles  in  the 
Thick-Walled  Detector  Compared  to  Calculations  for  Particles  with  L  =  149  keV  /mT1 


Data 

(keV  /mi1) 

£  =  L  •  pathlength 
(keV  pm'1) 

yF 

79 

150 

yD 

143 

169 
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pire  4-13.  Data  from  Experiment  with  Fe  Particles  in  Thick- Walled  Detector 
Simulating  1  Jim  [Open  Cirles  (O)  is  Data,  Solid  Line  is  Chatterjee  Model] 
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2.  Comparison  of  Data  to  Homogenous  Case  (Chatteriee  Model) 

Figure  4- 14  is  a  plot  of  the  ratio  of  mean  lineal  energy  measured  in  the  thin-  and  thick- 
walled  detectors  to  that  of  the  Chatterjee  model.  The  ordinate  of  the  primary  plot  is  on  a 
logarithmic  scale  while  the  abscissa  is  based  on  the  impact  parameter  squared,  that 
weighted  the  data  to  a  ^-random  irradiated  detector.  The  abscissa  of  the  small  inset  plot  is 
on  the  same  scale  as  the  primary  plot,  but  has  a  truncated  linear  ordinate  scale.  The  inset 
plot  emphasizes  the  smaller  ratios  for  impact  parameters  near  the  center  of  the  detector. 

The  ratio  of  the  data  to  the  model  is  small  for  intervals  of  impact  parameter  through  the 
center  of  the  cavity,  with  the  thick-walled  detector  being  slightly  higher.  The  enhanced 
response  in  the  thick-walled  detector  over  the  thin-walled  in  this  region  is  likely  due  to  a 

higher  number  of  6-rays  produced  in  the  front  wall  of  the  detector. 

The  ratio  of  the  data  to  the  model  is  the  greatest  at  the  cavity/wall  interface,  with  the 
highest  ratios  due  to  grazing  events.  The  ratio  of  the  data  to  the  model  was  high  for 
trajectories  through  the  wall-only,  with  the  response  of  the  thick-walled  detector  being 
higher  than  that  of  the  thin.  This  result  would  be  expected  since  the  pathlength  through 
wall  material  is  greater  for  these  trajectories  in  the  thick-wall  detector  as  compared  to  the 

thin  (see  Fig.  3-4)  and  would  result  in  a  greater  number  of  6-rays  being  produced. 

3.  Charged  Particle  Equilibrium  (CPE)  in  the  Thick-Wailed  Detector 

“Charged  particle  equilibrium  (CPE)  exists  for  a  volume  if  each  charged  particle  of  a 
given  type  and  energy  leaving  the  volume  are  replaced  by  an  identical  particle  of  the  same 
energy  entering,  in  terms  of  the  expectation  values,”  according  to  Attix  (1986).  For  walled 
cavity  detectors,  CPE  can  be  approximated  if  the  thickness  of  the  wall  is  greater  than  the 
maximum  penetration  range  of  secondary  particles  produced  and  the  wall  material  is  of 
approximately  the  same  atomic  composition  (to  provide  equivalent  mass  stopping  power  in 
both  media).  For  the  detectors  used  in  this  work,  the  wall  material  and  fill  gas  had  very 
similar  atomic  composition.  For  1  GeV/nucleon  particles,  the  maximum  kinetic  energy  of 
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Figure  4-14.  Ratio  of  Measured  Values  of  Mean  Lineal  Energy  to 
Chatterjee  Model  for  the  Thin-  and  Thick-Walled  Detectors  to  Fe 
[Simulated  Diameter  =  1  pim] 
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an  electron  is  about  3  MeV  and  has  an  approximate  CSDA  range  of  1.5  cm  in  water. 
However,  the  fraction  of  energy  transferred  to  the  medium  by  electrons  of  this  energy  is 
very  small.  The  thick-walled  detector  has  a  wall  thickness  of  2.54  mm  and  will  completely 
stop  all  electrons  with  energy  lower  than  0.76  MeV.  It  is  assumed  here  that  approximate 
CPE  can  be  achieved  in  the  thick-walled  detector. 

The  data  in  Fig  4-13  can  be  used  to  determine  if  the  enhancement  in  energy  deposition 

by  6  -rays  created  in  the  wall  is  sufficient  to  create  CPE  for  a  //-random  fluence  of  Fe 

particles.  Calculations  were  performed  according  to  the  method  described  in  Appendix  F. 
The  mean  lineal  energy  deposited  in  the  cavity  for  a  //-random  fluence  corresponded  to  an 
effective  L  of  154  keV  //m'1.  For  CPE,  this  value  should  equal  L.  The  value  obtained  here 
is  3  %  higher  than  expected  for  CPE.  Possible  reasons  for  the  difference  are:  a)  calibration 

uncertainties,  b)  position  uncertainties,  and  c)  a  possible  difference  between  W  for  the  a- 

particle  calibration  source  and  the  Fe  ions. 

4.  Calculated  Mean  Quality  Factors  Mean  quality  factors,  Q ,  for  the 
response  of  the  thick-walled  detector  to  a  //-random  Fe  particles  were  calculated  using: 

a.  ICRU  40  and 

b.  ICRP  26  and  60  with: 

(1)  L  =  y 

(2)  L  =  8/9  y  (Stinchcomb  and  Borak  1982). 

Identical  calculations  were  performed  for  //-randomly  incident  particles  on  the  detector,  but 
without  radial  energy  loss  (i.e.  e  =  L  •  pathlength  through  the  cavity).  Calculations  of  Q 

for  other  ions  found  in  the  galactic  cosmic  radiation  spectrum  were  made.  The  calculations 
were  made  for  the  same  velocity  of  Fe,  but  with  L  and  detector  response  scaled  according 
to  the  Z2  of  the  ion  compared  to  that  of  Fe.  This  scaling  procedure  should  provide 
reasonable  accuracy  because  the  doubly  differential  electron  spectra  of  ions  of  high 
velocity  scales  roughly  according  to  Z2  (ICRU  1996).  Details  of  the  method  of  calculation 
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are  provided  in  Appendix  F.  Tables  4-3  and  -4  provide  a  summary  of  die  results  of  the 
calculations  for  the  data  and  model,  respectively. 

For  the  experimental  data,  ICRU  40  provided  higher  Q  (for  all  ions  but  He)  than  the 
two  calculations  using  ICRP  26.  For  ICRU  40  compared  to  ICRP  60,  some  Q  were 
higher  while  some  were  lower.  For  the  ions  without  any  6-ray  losses,  ICRU  40  provided 

higher  Q  values  than  either  calculation  with  ICRP  26  and  60,  except  He. 

Figures  4-15  and  -16  are  bar  graphs  of  the  percent  discrepancy  in  the  calculated  Q  of 
the  data  from  Tables  4-3  and  -4.  For  the  experimental  data,  the  percent  discrepancy  was 
negative  or  near  zero  for  both  ICRP  report  recommendations.  The  calculations  of  Q  using 
y  =  L  provided  better  agreement  to  the  actual  Q  than  using  L  =  8/9  y.  For  the  ions  without 

any  6-ray  losses,  generally,  using  y  =  8/9L  provided  a  lower  absolute  percent  discrepancy 
than  using  y  =  L. 

2  and  3  «m  Simulated  Diameter  Experiments  -  Fe  Only  -  (Thick-Walled) 

Figure  4-17  contains  the  composite  distribution  of  mean  lineal  energy  from  Fe  for  the 
thick-walled  detector  for  simulated  diameters  of  2  and  3  pim.  The  shape  of  the  two 
distributions  are  similar  to  that  collected  for  a  1  pim  simulated  diameter.  Table  4-5  contains 
summary  statistics  for  the  distributions  and  those  from  the  experiment  with  a  1  jim 
simulated  diameter.  The  integrated  dose  for  the  experiments  with  the  detector  simulating  1 
and  2  ft m  diameter  sites  were  3  %  higher  than  expected  for  CPE,  while  the  calculation  for 
the  experiment  where  the  detector  simulated  a  3  y,m  diameter  site  was  close  to  that  expected 
for  CPE.  Also,  for  the  data  from  the  experiments  with  the  detector  simulating  1  and  2  ftm 
diameter  sites,  there  is  close  agreement  in  the  frequency  mean  and  dose  mean  lineal  energy 
values.  The  experiment  where  the  detector  was  simulating  3  pim  had  frequency  mean  and 
dose  mean  lineal  energy  values  about  3  %  lower  than  those  of  the  other  two  experiments. 

Schrewe  et  al.  (1988)  expected  uncertainties  in  calibrations  with  internal  a-  radiation 
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Table  4-3.  Calculated  Mean  and  True  Quality  Factors  for 
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He  (2)  0.23 


He  (2) 


Figure  4-15.  Percent  Discrepancy  of  Calculated  Mean  Quality  Factors  from  Mean  Lineal  Energy 
Distribution  of  Thick-Walled  Detector  to  Known  Quality  Factors  from  ICRP  26  and  60 
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Figure  4-16.  Percent  Discrepancy  of  Calculated  Mean  Quality  Factors  from  Particles 
Without  Radial  Energy  Loss  to  Known  Quality  Factors  from  ICRP  26  and  60 
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Mean  Lineal  Energy  (keV/^m)  Mean  Lineal  Energy  (keV/^m) 


Figure  4-17.  Data  From  Experiment  with  Fe  Particles  in 
the  Thick-Walled  Detector  Simulating:  a.  2  jim  and  b.  3  }im 
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Table  4-5.  Computed  yp  and  yp>  Based  on  Measured  Response  to  Fe  Particles  in 
the  Thick-Walled  Detector  and  Integrated  Dose  Compared  to  that  Expected  for  CPE 


Simulated  Diameter  (^m) 

Parameter 

1 

2 

3 

y  F  (keV  jtm1) 

79 

79 

76 

yD  (keV  j^m1) 

143 

145 

140 

Ratio:  Data/CPE 

1.03 

1.03 

1.00 

calibration  sources  to  be  on  the  order  of  a  few  percent.  Therefore,  determination  of  the 
differences  between  the  three  experiments  is  beyond  the  inherent  resolution  of  the 
calibration  technique.  It  is  speculated  that  there  exists  a  calibration  bias  between  the  3  pim 
experiment  and  both  the  1  and  2  pim  experiments. 

Figure  4-18  contains  a  plot  of  mean  lineal  energy  vs.  impact  parameter  for  Fe 
trajectories  through  the  cavity  for  the  three  simulated  diameters.  Standard  error  in 
estimation  of  the  mean  bars  for  individual  data  points  were  typically  about  0.5  keV  pim1 
and  were  not  displayed  because  they  were  very  small.  For  trajectories  through  the  center  of 
the  detector,  higher  mean  lineal  energy  was  associated  with  lower  simulated  site  diameter. 
But,  for  trajectories  closer  to  the  cavity/wall  interface,  the  difference  was  less  apparent. 

CV  Among  Thick-Walled  Experiments:  1.  2.  and  3  fim 
The  CV  was  calculated  for  distributions  of  lineal  energy  for  particle  trajectories  through 
the  center  of  the  detector  for  the  three  simulated  diameters  in  the  thick-walled  detector. 

CV’s  were  computed  for  three  separate  intervals  of  impact  parameter:  0-1,0- 1.5,  and 

0  - 1.75  mm  with  the  results  in  the  bar  graph  of  Figure  4-19.  The  error  bars  represent  la 

standard  error  in  the  estimate  of  the  mean.  The  CV  of  the  distributions  are  lower  for  larger 
simulated  site  diameter.  This  result  is  expected  because  the  variance  in  the  distribution  of 
lineal  energy  caused  by  energy  straggling  is  lower  for  higher  mean  energy  deposition. 
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Figure  4-18.  Mean  Lineal  Energy  vs.  Impact  Parameter 
for  Fe  Trajectories  through  the  Cavity  of  Thick-Walled 
Detector  [Simulated  Diameters  of  1, 2,  and  3  jim ; 
Standard  Error  in  the  Data  about  0.5  keV//im  (not  shown)] 
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For  the  experiments  where  the  detector  simulated  diameters  of  2  and  3  pim,  the  CV 
increased  as  the  interval  of  impact  parameter  increased.  This  was  due  to  the  increasing 
influence  of  variance  in  pathlength  for  larger  intervals  of  impact  parameter. 


Figure  4-19.  Coefficient  of  Variation  for  Distributions  of  Mean  Lineal  Energy  for  Three 
Impact  Parameter  Intervals  in  the  Thick- Walled  Detector  (1, 2,  and  3  ]4i n) 

10-, - — - - 


o  -  1  mm  0  -  1.5  mm  0  -  1.75  mm 

Impact  Parameter  Interval  (mm) 

Analysis  of  Events  with  Charge-Changing  Interactions  in  the  TEPC  Wall 


General 

The  results  presented  to  this  point  have  concentrated  on  the  analysis  of  Fe  trajectories 
through  the  TEPC.  By  rejecting  events  that  did  not  have  good  correlation  between  the 
position  values  in  PSD1  and  PSD2,  the  analysis  had  a  reduced  potential  of  including 
particles  that  were  involved  in  nuclear  interactions.  These  interactions  include  elastic 
scattering,  neutron(s)-only  loss  in  the  projectile,  and  target-only  fragmentation.  The 


75 


analysis  provided  here  is  for  events  involved  in  charge-changing  interactions  that  occurred 
between  the  PSD  pairs,  with  the  likeliest  location  being  in  the  wall  of  the  TEPC. 

Because  the  interactions  are  assumed  to  occur  in  the  wall  of  the  TEPC,  particle 
trajectory  information  determined  from  the  combination  of  signals  from  PSD1  and  PSD2 
would  be  incorrect.  Therefore,  impact  parameters  for  this  analysis  were  based  solely  on 
PSD1  signals.  As  such,  uncertainties  in  impact  parameter  are  typically  about  3.5  times  that 
of  impact  parameter  determined  with  PSD1  and  PSD2  (see  Appendix  D).  Due  to  large 
uncertainties  in  impact  parameter,  the  analysis  was  based  on  large  intervals  of  impact 
parameter.  The  response  of  the  detector  for  this  analysis  is  due  to  energy  deposition  from 
recoil  nuclei  and  electrons.  The  analysis  of  the  data  from  die  experiments  with  Fe-only, 
was  based  on  a  relatively  uniform  detector  response  at  a  specific  impact  parameter. 
However,  for  this  analysis,  the  influence  of  recoil  nuclei  on  energy  deposition  is  not  well 
known  for  TEPC  detectors.  Also,  for  identical  trajectories,  detector  response  could  be 
significantly  different  dependent  on  whether  the  nuclear  interaction  occurred  in  the  front  or 
back  wall  of  the  detector.  Evaluation  of  the  detector  response  for  trajectories  near  the 
cavity/ wall  interface  and  wall-only  could  not  be  as  detailed  as  in  the  case  of  the  Fe-only 
experiments.  Because  the  technique  of  determining  mean  impact  parameter  based  on 
selecting  an  interval  of  detector  response  may  provide  ambiguous  results.  Experimental 
results  are  for  the  thin- walled  detector  simulating  1  pi m. 

Analysis 

1*  Scatter  Plot  TEPC  Response  vs.  PSD  vt  Figure  4-20  is  a  plot  of 
energy  deposited  in  PSD  y2,  the  detector  directly  downstream  of  the  TEPC.  The  plot  is  for 
events  with  impact  parameter  less  than  2  mm  (as  estimated  from  PSD1  only).  From  the 
plot,  it  is  apparent  that  the  most  abundant  product  fragment  is  manganese  (Mn),  as  expected 
from  experiment  (see  Appendix  C).  The  mean  energy  deposited  in  PSD  y2  was  170.6 
MeV,  about  60  %  of  the  mean  energy  deposited  by  Fe  particles. 
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Figure  4-20.  AE  Y2  for  Fragments 
Produced  in  the  Wall  of  the  TEPC 
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Figure  4-21  is  a  scatter  plot  of  the  detector  response  vs.  energy  deposited  in  PSD  y2, 
for  events  with  estimated  impact  parameter  less  2  mm.  The  scatter  plot  has  two  obvious 
characteristics  shown  by  the  selected  rectangular  regions. 

The  upper  rectangular  region  is  characterized  by  a  detector  response  independent  of 
energy  deposited  in  PSD  y2.  This  region  corresponds  to  events  that  enter  the  detector  as 
Fe,  have  an  energy  deposition  as  expected  for  Fe  particles,  but  are  involved  in  a  charge¬ 
changing  nuclear  interaction  in  the  back  wall  of  the  detector.  Because  target  fragments  are 
largely  forward  directed,  recoil  nuclei  would  not  be  expected  to  deposit  energy  in  the 
detector,  consistent  with  the  observation  made  here. 

The  diagonally  oriented  rectangle  corresponds  to  events  where  the  response  of  the 
detector  is  generally  proportional  to  the  energy  deposited  in  PSD  y2.  These  events  appear 
to  have  been  involved  in  a  charge-changing  interaction  before  or  in  the  forward  wall  of  the 
detector.  Because  the  response  of  the  detector  is  proportional  to  energy  deposited  in  the 
downstream  silicon  detector,  it  does  not  appear  that  significant  enhancement  of  the  detector 
response  occurs  from  recoil  nuclei. 

There  are  many  particles  that  generally  do  not  fit  into  the  two  rectangular  regions  of 
Fig.  4-21.  The  largest  number  of  particles  have  AE  y2  values  about  250  MeV, 

corresponding  to  Mn.  A  group  of  these  particles  have  a  detector  response  (TEPC)  <  15 
keV  jm'\  while  another  group  has  a  response  from  100  to  140  keV  /mi'1.  Because  the 
detector  response  is  less  than  that  expected  for  Mn  or  Fe  particles  traversing  the  center,  it  is 
speculated  that  these  particles  traversed  the  detector  near  the  cavity/wall  interface  or  in  the 
wall-only.  These  events  could  be  the  result  of  uncertainties  in  impact  parameter.  But  more 
likely,  these  events  are  the  result  of  nuclear  interactions  in  the  downstream  portion  of  PSD 
xt  or  the  air  between  PSD  x,  and  the  detector  (TEPC).  If  either  of  these  occurred,  the 
particle  would  have  scattered,  possibly  causing  its  true  impact  parameter  to  be  significantly 
different  from  that  predicted  by  the  signal  in  PSD1. 
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Figure  4-21.  Detector  Response  vs.  AE  y2  for 
Fragments  Produced  in  the  Wall  of  the  TEPC 
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For  fragmentation  events  that  have  Mn  as  a  product,  there  is  a  reasonably  probability 
that  an  interaction  could  have  occurred  in  the  downstream  portion  of  PSD  without  an 
appreciable  decrease  in  the  energy  deposited,  because  the  L  of  Mn  is  about  92.4  %  of  Fe. 
For  fragmentation  events  that  have  Mn  as  a  product,  and  occur  at  a  penetration  depth  at 
least  53  %  of  the  way  through  the  PSD  x,,  they  will  have  an  energy  deposition  that  is 

within  2.5  o  of  the  mean  of  the  Fe  peak  (if  energy  deposition  by  the  Fe  and  Mn  are  equal  to 

their  expectation  values).  Forty-seven  percent  of  the  thickness  of  a  silicon  PSD  roughly 
corresponds  to  0.8  mm  of  A- 150  plastic.  From  Fig.  3-4,  the  total  wall  thickness 
encountered  for  a  trajectory  through  the  center  of  the  detector  is  2.54  mm.  Thus,  for  these 
trajectories,  about  25  %  of  the  fragmentation  events  that  have  Mn  as  a  product  will  occur  in 
the  downstream  portion  of  PSD  xt  and  falsely  be  attributed  to  a  TEPC  wall  interaction.  For 
trajectories  through  cavity/wall  interface,  about  9  %  of  the  fragmentation  events  that  have 
Mn  as  a  product  will  occur  in  the  downstream  portion  of  PSD  xt  instead  of  the  wall  of  the 
TEPC.  For  fragmentation  events  that  have  Cr  as  a  product,  about  16  and  5  %  of  the 
interactions,  respectively  for  trajectories  through  the  center  and  cavity/wall  interface,  will 
occur  in  the  downstream  portion  of  PSD  xl  rather  than  the  wall  of  the  TEPC. 

2.  TEPC  Response  vs.  Impact  Parameter  Distributions  of  detector 
response  were  created  for: 

a.  events  where  Fe  was  identified  in  PSD1  and  PSD2, 

b.  events  where  Fe  was  identified  in  PSD1  and  but  fragment  in  PSD2,  and 

c.  events  where  Fe  was  identified  in  PSD1  and  but  a  Mn  fragment  was 
identified  in  PSD2. 

The  last  category  was  included  because  Mn  fragments  were  abundant  and  would  be 
expected  to  have  an  energy  deposition  in  the  TEPC  close  to  that  of  Fe  (92.4%). 

Figure  4-22  contains  distributions  of  detector  response  for  impact  parameters  less  than 
2  mm.  Figure  4-22a  is  for  Fe  events,  while  Fig.  4-22b  is  for  all  charge-changing 
fragmentation  events.  The  mean  detector  response  for  the  events  that  have  been  involved  in 
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Figure  4-22.  Detector  Response  for  Estimated  Impact  Parameters  less 
than  2  mm  [a.  Fe  and  b.  All  Charge-Changing  Fragmentation  Events] 
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charge-changing  fragmentation  is  about  72  %  of  the  mean  for  the  distribution  of  Fe  events. 
If  50  %  of  the  charge-changing  interactions  that  occur  in  the  TEPC  do  so  in  the  back  wall 
and  the  detector  response  for  these  events  is  the  same  as  Fe  events,  then  the  mean  response 
of  the  detector  to  events  with  charge-changing  interactions  that  occur  in  the  forward  wall  is 
44  %  of  Fe. 

Figure  4-23  contains  the  distribution  of  detector  response  for  impact  parameters 
between  6  and  7  mm.  Figure  4-23a  is  for  Fe  events,  while  Fig.  4-23b  is  for  all  charge¬ 
changing  fragmentation  events.  The  mean  detector  response  for  the  events  that  have  been 
involved  in  a  charge-changing  fragmentation  is  about  90  %  of  the  mean  for  the  distribution 
of  Fe  events.  The  ratio  of  mean  lineal  energy  for  these  distributions  is  significantly  higher 
than  that  of  the  ratio  for  events  through  the  center  of  the  detector. 

Figure  4-24  is  a  summary  for  all  distributions  of  detector  response  for  impact 
parameters  less  than  8  mm.  The  mean  detector  response  to  fragmentation  events  that  have 
Mn  as  a  product  was  higher  than  the  mean  response  that  included  all  charge-changing 
fragmentation  products.  The  detector  response  to  Fe  particles  was  higher  than  the  mean  for 
Mn  and  all  fragments,  except  for  impact  parameter  intervals  near  the  cavity/wall  interface 
and  in  the  wall-only. 

3.  Events  with  TEPC  Response  Greater  than  190  keV  urn'1 

The  mechanism  responsible  for  production  of  the  events  with  large  values  of  energy 
deposition  was  of  interest  for  the  experiments  were  Fe  had  traversed  the  detector.  It  is 
possible  that  that  these  events  had  unusually  high  energy  deposition  due  to  a  contribution 
from  recoil  nuclei.  The  events  evaluated  here  are  known  to  have  been  involved  in 
fragmentation  interactions  in  the  wall  of  the  detector.  A  comparison  of  the  frequency  of 
production  of  these  high  energy  deposition  events  with  that  for  Fe  traversals  through  the 
detector  is  provided  here.  This  analysis  will  assist  in  determining  the  possibility  that  recoil 
nuclei  are  responsible  for  high  energy  deposition  events  where  Fe  traverses  the  detector. 
Figure  4-25  is  a  histogram  of  the  percent  of  total  events  with  y  >  190  keV  pt m'1  for  Fe 
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Figure  4-23.  Detector  Response  for  Estimated  Impact  Parameters  between 
6  and  7  mm  [a.  Fe  and  b.  A1I  Charge-Changing  Fragmentation  Events] 
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Figure  4-24.  Mean  Lineal  Energy  vs.  Impact  Parameter  for  Fe  Particles  through 
the  Detector,  All  Charge-Changing  Fragmentation  Events  in  the  Wall,  and 
Fragmentation  Events  with  Mn  as  a  Product 
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particles  through  the  detector,  all  charge-changing  fragmentation  events  that  occur  in  the 
wall  of  the  detector,  and  fragmentation  events  that  have  Mn  as  a  product.  Vertical  error  bars 

represent  an  estimate  of  the  la  standard  error  in  the  estimate  of  the  mean  using  Poisson 

statistics. 

For  events  with  impact  parameters  less  than  2  mm,  the  percent  of  total  events  for  the 
charge-changing  fragmentation  events  was  about  half  of  that  for  Fe  events  through  the 
TEPC.  This  naturally  follows  from  the  observations  and  conclusions  drawn  from  the 
scatter  plot  of  Fig.  4-21.  For  the  same  range  of  impact  parameters,  the  Mn  products  had  a 
lower  percent  of  total  events  with  lineal  energy  greater  than  190  keV  ptm'1,  than  that  for  all 
charge-changing  interactions,  but  uncertainties  in  this  value  are  high  with  a  sample  size  of 

only  6  events  and  bias  introduced  by  difficulties  in  interpreting  AE  in  PSD  xx  for  Mn 

products. 


84 


For  events  near  the  cavity/wall  interface  (impact  parameters:  5-6  and  6-7  mm),  the 
percent  of  total  events  with  lineal  energy  greater  than  190  keV  pm -1  was  higher  for  the  Mn 
fragments  than  for  the  Fe  events  through  the  detector. 


Figure  4-25.  Percent  of  Total  Events  with  Lineal  Energy  >  190  keV  pim  1 
[Fe  Particle  through  the  Detector,  All  Charge-Changing  Fragmentation  Events 
that  Occur  in  the  Detector  Wall,  and  Fragmentation  Events  with  Mn  Products] 
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Detector  Response  to  Fragments  Produced  in  the  Copper  Target 
Event  Selection 

Events  with  energy  deposition  in  PSD1  2.5  a  less  than  the  mean  energy  deposition  for 


Fe  were  selected  for  analysis.  Figure  4-26  is  a  plot  of  — —  these  events.  The 

AExj 

distribution  has  a  prominent  peak  and  a  broad  distribution  of  events  with  a  ratio  of  energy 
deposition  lower  than  the  peak.  The  broad  peak  is  believed  to  be  the  result  of  one  of  two 
factors.  First,  some  projectile  fragments  (within  the  same  event)  that  penetrate  the  active 
volume  of  PSD  xx  do  not  penetrate  the  active  volume  of  PSD  y2.  The  magnitude  of  this 
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Figure  4-26.  Quotient  (AE  y2)  by  (AE  xl)  for 
Fragmentation  Events  Produced  in  the  Copper 
Target  [Bars  Show  Events  Selected  for  Analysis] 
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effect  is  expected  to  increase  as  the  value  of  charge-loss  in  the  projectile  increases,  because 
there  will  be  an  increasing  number  of  projectile  fragments.  Second,  some  of  the  projectile 
fragments  may  be  involved  in  additional  fragmentation  events  between  PSD  Xj  and  PSD  y2. 

A  plot  of  the  same  quotient  distribution  was  generated  for  Fe  particles  identified  in 
PSD1  and  PSD2  (not  displayed).  This  data  was  normally  distributed  with  mean  and 
standard  deviation  of  0.957  and  0.024,  respectively.  The  location  of  the  peak  in  Fig.  4-26 
was  very  close  to  that  observed  for  Fe  particles.  The  analysis  was  restricted  to  events  with 
quotient  values  between  the  bars  of  Fig.  4-26.  Because  the  number  of  fragments  produced 
in  die  copper  target  was  small,  the  analysis  of  detector  response  was  limited  to  particles 
with  impact  parameter  less  than  3  mm. 

Detector  Response 

Figure  4-27  is  a  scatter  plot  of  the  detector  response  vs.  the  energy  deposited  in  PSD 
y2.  The  data  has  a  linear  correlation  between  the  response  of  the  two  detectors.  One  event 
had  an  unusually  high  lineal  energy  response  in  the  detector  (TEPC),  near  270  keV  pirn1 
and  was  assumed  to  be  the  result  of  a  wire  hit.  The  plot  also  contains  numerous  events 
with  low  lineal  energy  response  in  the  detector.  As  discussed  earlier  for  evaluation  of  the 
response  of  the  thin-walled  detector  to  Fe: 

1.  these  events  probably  crossed  the  wall-only  of  the  detector  and  have  low 
detector  response,  and 

2.  were  involved  in  a  scattering  event  that  caused  the  impact  parameter  to  be 
incorrect. 

The  number  of  these  events  observed  in  Fig.  4-27  is  statistically  consistent  to  the  number 
observed  in  the  analysis  of  Fe  events. 

Figure  4-28  is  a  scatter  plot  of  the  quotient  of  detector  (TEPC)  response  by  the  energy 
deposited  in  PSD  y2.  The  distribution  of  events  appears  to  be  distributed  about  the  ratio  of 

0.6,  independent  of  AE  y2.  A  linear  regression  on  the  data  in  this  form  was  not  performed 

because  the  coefficients  would  have  been  dominated  by  the  response  of  the  heavier 
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TEPC:  y(keV///m) 


Figure  4-27.  Detector  Response  vs.  (AE  y2) 
for  Fragments  Produced  in  the  Copper  Target 
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Figure  4-28.  Quotient  of  (Detector 
Response  by  AE  y2)  vs.  (AE  y2)  for 
Fragments  Produced  in  the  Copper  Target 
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fragments  of  Mn  and  Cr.  These  fragments  were  greater  in  number  than  all  others 
combined. 

Since  a  linear  regression  of  the  data  in  Fig.  4-28  could  have  produced  misleading 
results.  The  ratio  of  detector  response  to  the  energy  deposited  in  PSD  y2  was  evaluated  for 
groups  of  fragments.  Table  4-6  contains  a  summary  for  the  fragments  and  Fe  particles 
through  the  detector.  The  events  that  were  postulated  to  have  been  involved  in  scattering 
events  in  the  detector  wall  (15)  and  the  one  wire  hit  were  omitted  from  the  analysis. 

Table  4-6.  Quotient  of  Detector  Response  by  AE  y2  for  Various  Groups  of  Fragments 


Particle(s) 

Z 

TEPC  /  AE  y2 
(mean) 

n 

Fe 

26 

0.6015  ±  0.0006 

16605 

Mn  -  Cr 

25-24 

0.604  ±  0.003 

372 

V-Sc 

23-21 

0.601  ±  0.006 

146 

Ca  -  Ar 

20-  18 

0.586  ±  0.008 

111 

T-tests  on  the  mean  of  the  three  fragment  groupings  with  that  of  the  Fe  data  did  not 
demonstrate  a  difference  at  a  significance  level  of  a/2  =  0.05;  but  at  a  significance  level 

of  a/2  =  0.025,  the  Ca  -  Ar  fragment  group  was  different 
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Chapter  5:  Discussion 


Fraction  of  L  Absorbed  in  the  Detector 

The  mean  energy  deposited,  e,  in  the  sensitive  volume  for  Fe  particle  trajectories 

through  the  center  of  the  thin-walled  detector  was  11 1.9  +  0.3  keV  (y  =  167.8  +  0.4 
keV  pa n'1)  for  a  1  pan  simulated  diameter.  This  accounts  for  about  75  %  of  L  for  incident 
Fe  particles.  The  mean  value  for  the  thick-walled  detector  was  1 18.1  +  0.5  keV  (y  =  177.1 
+  0.5  keV  pun'1)  corresponding  to  about  79  %  of  L.  The  model  described  by  Chatterjee 
etal.  (1976)  only  describes  energy  deposition  in  the  radial  direction  for  homogenous 
cylinders.  These  cylinders  were  superimposed  on  a  spherical  geometry  with  the  central 
axis  of  the  cylinder  corresponding  to  various  impact  parameters  (see  Fig.  4-12).  For 
trajectories  through  the  center  of  the  detector,  the  model  predicted  that  70  %  of  L  would  be 
deposited  in  a  1  pan  sphere.  If  one  assumes  that  the  model  is  correct,  then  5  %  of  the 

contribution  to  the  measured  values  of  energy  deposition  is  from  8-rays  produced  for  the 

thin  wall  and  9  %  from  the  thick  wall. 

Metting  etal.  (1986)  used  a  nearly  wall-less  cylindrical  TEPC  that  simulated  a  unit 
density  site  of  1.3  pan  to  collect  measurements  near  the  trajectory  of  600  MeV/nucleon  S6Fe 
ions.  The  authors  reported  a  loss  of  about  23  %  in  L  but  were  not  sure  if  it  was  due  to 

8-ray  losses  or  non-linear  response  of  the  detector.  The  value  of  energy  escaping  is  close 

to  that  predicted  by  the  Chatteijee  model  (1976)  for  600  MeV/nucleon  ions. 

Dicello  et  al.  (1991)  collected  microdosimetric  spectra  with  a  wall-less  cylindrical  TEPC 
in  beams  of  Fe  (535  MeV/nucleon),  Ar  (570  MeV/nucleon),  and  Ne  (557  MeV/nucleon). 
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The  detector  had  a  mean  pathlength  of  0.4  pim.  The  authors  reported  that  for  direct 
traversal  primaries,  about  20  %  of  the  energy  transferred  was  not  absorbed  in  the  detector. 

Large  Energy  Deposition  Events 
General 

The  events  with  the  highest  energy  deposition  for  Fe  particles  through  the  detector  did 
not  correspond  to  those  that  traversed  the  center,  rather  they  were  due  to  those  that  just 
grazed  the  cavity/wall  interface.  It  is  of  interest  whether  energy  deposition  for  these  events 
is  the  result  of  electron  ionization  or  that  from  both  electrons  and  recoil  nuclei. 

ICRU  (1983)  published  spectra  from  walled  and  wall-less  detectors  collected  by 
Rodgers  et  al.  (1973)  for  279  MeV/nucleon  nitrogen  ions.  The  ICRU  presumed  that  the 
very  large  energy  deposition  events  in  the  walled  counter  were  either  from  spallation  in  the 

wall  of  the  detector  or  a  contribution  from  6-rays  at  low  lineal  energies.  The  experimental 

set-up  used  by  Rodgers  et  al.  did  not  have  the  capability  to  separate  the  events  that 
traversed  the  detector  unfragmented  from  those  that  were  involved  in  fragmentation 
interactions  in  the  wall  of  the  detector.  The  analysis  presented  here  allowed  separation  of 
the  two  respective  types  of  events.  The  data  presented  below  provides  convincing  evidence 
for  1  GeV/nucleon  Fe,  that  the  mechanism  of  energy  deposition  for  the  events  with  the 
highest  energy  deposition  is  not  dominated  by  that  from  recoil  nuclei,  but  rather  that  from 
soft  collision  electrons.  Data  presented  was  for  the  thin-walled  detector. 

Nuclear  Interaction  Cross-Sections 

The  first  step  in  the  evaluation  of  these  events  is  a  review  of  nuclear  interaction  cross- 
sections  to  determine  if  the  probability  for  nuclear  interactions  in  the  wall  of  the  detector  is 
sufficiently  high  to  meet  the  frequency  of  these  events  observed  in  the  experimental  data. 

In  tiie  data  analysis  for  Fe  ions  through  the  detector,  events  with  nuclear  interactions  in  the 
wall  of  the  detector  could  not  be  rejected  from  the  analysis  if  the  incident  Fe  particle  did  not 
experience  a  loss  of  charge.  These  include:  a)  elastic  scattering  (Wilson  and  Costner, 
1975),  b)  fragmentation  of  the  incident  particle  with  neutron(s)  loss-only  (Zeitlin  et  al 
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1997),  and  c)  fragmentation  of  the  target  nucleus-only  (Heckman  etal.  1978).  An  analysis 
using  cross-sections  from  the  listed  references  determined  that  the  probability  of  a  nuclear 
interaction  in  the  wall  for  trajectories  near  the  cavity/wall  interface  was  about 
5  %.  This  analysis  assumed  that  only  half  of  the  charged  recoils  from  elastic  scattering  in 
the  wall  could  reach  the  sensitive  volume.  In  the  case  of  target-only  fragmentation,  the 
fraction  measured  for  160  nuclei  was  applied,  because  the  literature  did  not  have 
experimental  data  for  56Fe.  This  estimate  is  likely  to  be  higher  than  that  of  56Fe,  since  160 
ions,  of  lower  Z  and  doubly-magic,  are  likely  to  be  more  stable  (a  trend  evidenced  in  the 
Heckman  data,  see  Appendix  C  for  summary). 

From  evaluation  of  Fig.  4-7c,  5  %  of  the  events  have  y  >  190  keV  /mi'1.  As  noted 
above,  the  probability  of  an  undetected  nuclear  interaction  along  the  longest  pathlength  was 
about  5  %  for  those  that  could  direct  a  recoil  nuclei  into  the  cavity.  Therefore,  based  on  the 
nuclear  interaction  cross-section  data  alone,  it  is  possible  that  these  large  energy  deposition 
events  are  caused  by  recoil  of  target  atom  nuclei.  However,  the  probability  that  a 
significant  fraction  of  the  energy  deposited  for  these  events  is  from  recoil  nuclei  is  unlikely, 
because: 

1.  recoil  nuclei  produced  in  the  downstream  hemisphere  of  the  detector  will  almost 
certainly  not  reach  the  cavity, 

2.  many  recoil  nuclei  produced  in  the  wall  will  be  stopped  in  the  wall,  and 

3.  for  proton  recoils,  the  maximum  energy  transfer  in  a  1  /<m  diameter  site  is 
98  keV  (y  =  147  keV  fim1),  assuming  the  recoil  nuclei  traverses  the  diameter. 

Other  events  in  the  vicinity  of  the  cavity/wall  interface  had  energy  deposition  of  about 
50  keV  (y  =  75  keV  ptm'1).  Energy  deposition  from  these  events  is  believed  to  be  from 
secondary  electrons  only.  Combining  this  with  the  maximum  energy  deposition  from  a 
proton  recoil,  energy  deposition  is  about  150  keV  (y  =  225  keV  /mV1),  a  value  lower  than 
that  of  over  two-thirds  of  the  events  with  y  >  190  keV  ptm'1  from  Fig.  4-7c.  Thus,  events 
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with  proton  recoils  alone  (elastic  scattering)  could  not  be  attributed  to  two-thirds  of  these 
events. 

Trajectory  Analysis 

Another  method  was  used  to  determine  if  events  with  high  energy  deposition  were  from 
electrons  or  the  sum  of  electrons  and  recoil  nuclei.  The  difference  in  position,  x2  -  xt, 
using  the  PSD  detectors  was  plotted  for  events  with  40  <  y  <  60  keV  }i m'1  and  y  >  240 
keV  pim'1  when  Fe  was  identified  both  upstream  and  downstream  of  the  detector.  Events 
with  y  >  240  keV  pim'1  were  chosen  because  it  limited  the  number  of  wire  hits  events  from 
the  analysis.  The  difference  in  position  was  also  evaluated  when  Fe  was  identified  using 

AE  in  PSD1,  but  a  fragment  of  Fe  was  identified  using  AE  in  PSD2.  This  data  is 

summarized  in  Table  5-1. 

The  standard  deviation  for  the  distribution  of  (Xj  -  xt)  was  significantly  higher  in  events 
that  were  known  to  have  been  associated  with  nuclear  interactions  between  PSD1  and 
PSD2.  The  standard  deviation  increased  with  the  magnitude  of  charge  loss  by  the  incident 
Fe.  For  fragmentation  events  with  Mn  products,  the  standard  deviation  was  over  twice  that 
of  the  events  with  y  >  240  keV  pim1.  Fragmentation  events  with  Mn  products  on  average 
are  the  least  violent  of  the  charge-changing  fragmentation  interactions  because  they  have  the 
smallest  overlap  of  projectile  and  target  nuclei.  Fragmentation  events  resulting  in 
neutron(s)  loss-only  and  target-only  fragmentation  have  momentum  transfer  of  nearly  the 
same  magnitude  as  fragmentation  events  with  Mn  products.  Thus,  the  standard  deviation 
of  the  (x2  -  xj  distribution  for  these  events  would  be  expected  to  be  similar  to  that  of 
fragmentation  events  with  Mn  products.  Therefore,  it  appears  unlikely  that  the  majority  of 
events  with  high  energy  deposition  from  Fe  through  the  detector  have  been  involved  in 
nuclear  interactions. 

Soft  Collision  Electron  Deposition  Near  the  Trajectory  of  Fe  Ions 

A  simple  analysis  was  performed  using  the  Chatterjee  Model  (1976)  to  further  assess 
the  mechanism  responsible  for  energy  deposition  from  events  with  the  highest  magnitude 
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Table  5-1.  Summary  of  (x2  -  xt)  Trajectory  Analysis 
as  Measured  Respectively  in  PSD2  and  PSD1 


Data  Analysis 

TEPC  Response 
Range  (keV//im) 

x2  -  Xx 

PSD1 

PSD2 

ji  (mm) 

a  (mm) 

Fe 

Fe 

y  >  240 

l 

-u 

• 

o 

1.8 

Fe 

Fe 

-  3.9 

1.8 

Fe 

Mn 

NA 

-2.9 

3.9 

Fe 

Cr 

NA 

-33 

4.1 

Fe 

y 

NA 

-3.6 

3.9 

Fe 

Ti 

NA 

-3.7 

3.9 

Fe 

Sc 

NA 

-3.7 

4.4  ■ 

Fe 

Ca 

NA 

-33 

4.1 

Fe 

K 

NA 

-3.9 

4.4 

Fe 

Ar 

NA 

-3.6 

4.7 

Fe 

Cl 

NA 

-33 

5.2 
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of  detector  response.  From  the  model  (Fig.  A-l),  the  core  radius  is  9  nm  for  1  GeV/ 
nucleon  particles.  The  model  assumes  that  about  50  %  of  the  energy  transferred  to  matter 

is  from  soft  collision  electrons  and  the  other  from  6-rays.  The  model  assumes  that  the  soft 

collision  electrons  deposit  their  kinetic  energy  in  the  core  region. 

Figure  5-1  is  a  cross-section  of  the  detector  cavity  and  depicts  an  Fe  particle  grazing  the 
cavity/wall  interface.  In  the  diagram,  three  distinct  segments  of  the  particle  trajectory 
through  the  wall  material  have  been  defined.  A  and  c  are  respectively,  the  upstream  and 
downstream  segments  of  the  trajectory  from  the  detector  midpoint.  The  length  of  segment 
b  was  limited  in  length  by  the  criterion  that,  “every  point  on  the  segment  had  a  potion  of  the 
core  extending  into  the  gas  cavity.”  For  a  core  radius  of  9  nm,  segment  b  had  a  length  of 
21.5  }im. 

Figure  5-1.  Fe  Grazing  the  Cavity/Wall  Interface  of  the  Detector.  [Segment  “a”  is 
Upstream  Portion  of  the  Trajectory,  Segment  “b”  is  Portion  of  Trajectory  with  the  Core 
Extending  into  the  Cavity,  and  Segment  “c”  is  Downstream  Portion  of  die  Trajectory] 


For  grazing  events,  the  energy  deposited  in  the  detector  will  be  assumed  to  be 
comprised  of  two  components:  energy  deposition  from  6-rays  produced  in  segments  a  and 

c,  and  from  soft  collision  electrons  in  segment  b.  The  contribution  from  the  6-rays  is 
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assumed  to  be  about  50  keV  (75  keV  /im'1).  This  is  an  energy  deposition  that  is 
characteristic  of  Fe  events  that  have  trajectories  near  the  cavity/wall  interface,  but  in  the 
wall-only  (see  Fig.  4-13). 

The  magnitude  of  energy  deposition  with  respect  to  impact  parameter  is  not  described 
by  the  Chatterjee  Model  for  the  core.  The  model  only  specifies  a  mean  value  of  energy 
deposition.  Therefore,  accurate  estimates  of  energy  deposition  in  the  cavity  from  soft 
collision  electrons  created  by  a  grazing  event  cannot  me  made  from  the  model.  Other 
models  of  energy  deposition  from  heavy  ions  also  suffer  large  uncertainty  in  energy 
deposition  for  impact  parameters  very  close  to  the  trajectory. 

A  rough  estimate  of  energy  deposition  is  made  with  the  following  assumptions: 

1.  75  keV  ptm'1  (one-half  of  L)  is  the  energy  transferred  by  soft  collision  electrons, 

2.  only  soft  collision  electrons  generated  in  segment  b  are  included, 

3.  50%  of  the  soft  collision  electrons  are  directed  away  from  the  cavity,  and 

4.  75  %  of  those  directed  toward  the  cavity  are  absorbed  in  the  wall,  without  any 
energy  deposited  in  the  cavity. 

The  amount  of  energy  transferred  in  segment  b  that  is  absorbed  in  the  wall  (assumption 
4)  cannot  be  predicted  precisely  from  the  Chatteijee  model.  In  the  model,  it  is  assumed  that 
energy  deposited  in  the  core  is  primarily  from  soft  collision  electrons  and  that  the  energy 
deposition  is  uniform.  The  authors  of  the  model  noted  that  the  spatial  distribution  of 
energy  deposition  in  the  core  was  not  known,  but  the  assumption  of  a  core  radius  was 
“introduced  mainly  for  conducting  the  analytical  process.”  Assumption  4  in  the  least  is 
conservative  to  our  analysis.  Because  for  segment  b,  the  average  distance  between  the 
primary  ion  track  and  cavity  is  about  one-half  the  core  radius,  leading  to  an  average  energy 
deposition  of  50  %  according  to  the  Chatteijee  model.  In  our  analysis,  an  average  energy 
deposition  of  only  25  %  is  assumed.  The  estimated  soft  collision  electron  energy 
deposition  in  the  cavity  is  200  keV  (y  =  300  keV  /mi'1). 
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Therefore,  from  this  simple  analysis,  one  may  conclude  that  the  high  energy  deposition 
events  could  be  the  result  of  energy  deposition  from  electrons  alone.  Furthermore,  for 

grazing  events,  soft  collision  electrons  contribute  more  to  energy  deposition  than  6-rays. 

Events  with  Charge-Changing  Interactions  in  the  Detector  Wall 

1 .  General  Based  on  the  foregoing  discussion,  it  appears  unlikely  that  recoil 
nuclei  are  responsible  for  the  high  energy  deposition  events  for  Fe  through  the  detector. 
Charge-changing  fragmentation’s  occur  with  only  a  few  percent  of  the  ^-random  Fe 
particles  incident  on  the  detector.  The  evaluation  of  these  events  also  reinforce  the 
conclusion  that  recoil  nuclei  are  not  responsible  for  the  majority  of  high  energy  deposition 
events  for  Fe  through  the  detector. 

2.  Mn  Data 

For  events  near  the  cavity/wall  interface  (impact  parameters:  5-6  and  6-7  mm),  the 
percent  of  total  events  greater  than  190  keV  yim1  was  greater  for  fragmentation  events  that 
had  Mn  products  than  for  1)  the  Fe  events  through  the  detector  and  2)  all  charge-changing 
fragmentation  events  (see  Fig.  4-25).  The  uncertainty  in  the  percent  value  for  Mn  is 
relatively  high  compared  to  the  mean  since  the  total  number  of  observations  is  low. 
Conclusions  about  the  difference  in  the  mean  of  the  data  for  Fe  events  through  the  detector 
and  that  for  fragmentation  events  that  have  Mn  products  cannot  be  made  with  statistical 
confidence. 

A  99  %  confidence  interval  for  the  percent  of  total  events  greater  than  190  keV  ptm1 
was  calculated  for  the  Mn  data  for  impact  parameters  between  6  and  7  mm.  The  high  side 
value  is  4.6  %.  As  discussed  in  Chapter  5,  about  9  %  percent  of  these  events  may  have 
occurred  upstream  of  the  detector  in  PSD  xr  Independent  of  these  small  errors,  it  appears 
inconceivable  that  a  significant  fraction  of  the  high  energy  deposition  events  for  Fe  through 
the  detector  (that  are  associated  with  a  non-charge-changing  nuclear  interaction  in  the  wall) 
could  be  due  to  recoil  nuclei,  if  at  best  only  5  %  of  the  events  that  are  known  to  have  been 
associated  with  fragmentation  events  that  have  Mn  products  had  high  energy  deposition. 
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3.  Bias 


Some  care  must  be  taken  in  interpreting  the  data  from  the  charge-changing  nuclear 
interactions  in  the  detector  wall.  For  the  Fe  particles  through  the  detector,  there  was  a 
uniform  particle  fluence.  But,  for  the  data  collected  for  charge-changing  fragmentation 
events  that  occurred  in  the  detector  wall,  the  particle  fluence  was  not  uniform  because: 

1.  the  pathlength  through  the  wall  material  of  the  detector  was  not  uniform  (see 
Fig.  3-4)  and 

2.  probability  for  a  nuclear  interaction  in  the  detector  is  nearly  proportional  to  the 
pathlength. 

Large  uncertainties  in  impact  parameter  for  the  data  from  this  experiment  prevented 
resolution  of  the  differences  in  fluence.  Because  the  longest  pathlength  through  wall 
material  occurred  at  the  cavity/wall  interface,  the  data  had  a  higher  particle  fluence  for  these 
trajectories.  This  bias  would  increase  the  fraction  of  events  with  y  >  190  keV  /un'1  higher 
than  that  for  a  uniform  fluence  of  particles.  Therefore,  for  the  evaluation  of  these  events 
provided  in  this  section,  any  bias  introduced  would  be  conservative  to  any  conclusions 
reached. 

Summary 

Absolute  determination  of  the  mechanism  responsible  for  the  events  with  the  highest 
energy  deposition  for  Fe  particles  through  the  detector  is  beyond  the  scope  of  this  work. 
But,  based  on  the  data  and  arguments  presented,  it  appears  that  the  majority  of  high  energy 
deposition  events  are  the  result  of  secondary  electrons  only. 

Implications  of  Wall  Effects  on  Radiation  Protection  Dosimetry 

General 

Walled  detectors  do  not  accurately  represent  single-event  energy  deposition  spectra  that 
would  be  observed  in  homogenous  media  from  high  velocity  ions.  Comparison  of  the 
experimental  data  to  the  Chatteijee  model  in  Fig.  4-12  shows  that  the  experimental  data  had 
higher  energy  deposition  than  the  model  for  trajectories  through  the  cavity  and  wall-only  of 
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the  detector.  The  mean  lineal  energy  for  the  Chatteijee  model  is  not  displayed  for  impact 
parameters  greater  than  the  wall  boundary.  But,  according  to  the  model  (Appendix  A), 
energy  deposition  in  the  sensitive  volume  will  occur  for  particle  trajectories  with  impact 
parameters  out  to  250  pim.  Integration  of  the  energy  deposited  from  these  indirect  events 

will  account  for  the  energy  lost  from  the  direct  events  due  to  the  escape  of  energetic  6-rays. 

In  homogenous  media,  indirect  events  will  account  for  a  large  fraction  of  the  total 
number  of  events  in  a  uniform  radiation  field  of  high  velocity  ions,  but  the  magnitude  of  the 
energy  deposition  for  individual  events  will  be  very  low.  Many  indirect  events  may 
involve  only  a  single  electron  crossing  the  sensitive  volume.  Walled  detectors  compress  all 
of  these  low  energy  deposition  events  into  considerably  fewer  events  with  higher  energy 
deposition  per  event.  For  this  reason,  application  of  the  ICRU  40  recommendations  for 
determination  of  mean  quality  factor  may  be  incorrect  for  walled  detectors. 

Mean  O  from  ICRP  26  and  60  Recommendations 

Currently,  NASA  is  applying  ICRP  recommendations  for  the  determination  of  Q  from 
microdosimetric  spectra  under  the  assumption  that  y  =  L.  On  an  event-by-event  basis,  y 
does  not  equal  L.  Based  on  Fe  particles  with  the  same  L  (see  Figs.  4-13, 4- 17a,  and  4- 
17b),  detector  response  (y)  varied  considerably  dependent  on  particle  trajectory.  For  Fe 
particles  through  the  detector  that  just  grazed  the  cavity/wall  interface,  lineal  energy  values 
were  higher  than  300  keV  pim'\  though  the  particles  responsible  for  these  events  had  an  L 
of  150  keV  ft m'1. 

For  the  experimental  data,  differences  in  calculated  Q  values  (using  ICRP  26  and 
y  =  L)  from  the  actual  Q  are  negative  but  less  than  6.2  %  (in  magnitude).  For  the  same 
ICRP  recommendations,  but  with  y  =  8/9  L,  the  difference  in  calculated  Q  from  the  actual 
Q  were  all  negative,  but  greater  in  magnitude  than  10  %  for  all  of  the  ions  but  He.  Thus,  at 
least  for  these  ions,  using  y  =  L  is  more  appropriate  for  calculation  of  Q  than  using  y  =  8/9 
L.  The  same  general  conclusion  was  reached  in  application  of  ICRP  60  recommendations. 
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For  the  particles  without  radial  energy  loss,  differences  between  calculated  Q  values 
and  the  actual  Q  had  magnitudes  as  high  as  9  %  for  the  ions  considered  under  ICRP  26 
recommendations  with  y  =  L.  For  the  same  ICRP  recommendations,  but  with  y  =  8/9  L, 
discrepancies  were  less  than  4.1%  for  all  of  the  ions  considered.  For  these  ions,  using 

y  =  8/9  L  is  more  appropriate  for  calculation  of  Q  than  using  y  =  L,  except  for  Fe.  The 
same  general  conclusion  was  reached  in  application  of  ICRP  60  recommendations,  again 
except  for  Fe. 

The  recommendation  to  use  y  =  8/9L  by  Ricourt  et  al.  (1981)  was  for  lineal  energy 
distributions  of  detectors  irradiated  with  neutrons  of  energies  less  than  a  few  MeV.  Recoil 
nuclei  from  neutrons  of  this  energy  will  have  low  velocity  and  have  little  radial  energy  loss 
for  site  diameters  in  the  micron  range.  Distortions  from  the  wall  effect  for  high  velocity 

ions  make  the  assumption  of  y  =  L  more  appropriate  for  calculation  of  Q  in  walled 
detectors. 

Fragments  Produced  in  the  Detector  Wall 

General 

The  summary  of  detector  response  to  the  charge-changing  fragmentation’s  produced  in 
the  detector  wall  is  displayed  in  Fig.  4-24.  As  noted  earlier,  caution  should  be  applied  in 
interpreting  the  data  for  intervals  of  impact  parameter  near  the  cavity/wall  interface  because 
of  uncertainties  in  impact  parameter  and  a  non-uniform  particle  fluence.  For  intervals  of 
impact  parameter  near  the  center  of  the  detector,  these  effects  introduce  less  uncertainty. 

Due  to  these  shortcomings,  assessment  of  the  contribution  of  recoil  nuclei  to  the  total 
energy  deposition  cannot  be  made  in  a  concise  manner.  The  discussion  and  observations 
will  be  qualitative  in  nature  only. 

Data  Analysis 

For  events  with  impact  parameter  less  than  2  mm,  the  mean  response  of  the  detector  to 
the  fragments  produced  in  the  wall  was  only  72  %  of  the  response  to  Fe  particles  through 
the  detector.  For  events  in  this  interval  of  impact  parameter  (see  Fig.  4-21),  it  appeared  that 
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detector  response  was  independent  of  any  recoil  nuclei  that  were  produced,  but  dependent 
only  on  whether  the  nuclear  interaction  occurred  in  the  front  or  rear  wall  of  the  detector  and 
the  degree  of  charge-reduction  in  the  projectile.  Of  the  697  events  in  Fig.  4-21,  only  one 
had  a  detector  response  greater  than  250  keV  pa n'1  (lineal  energy)  and  was  likely  due  to  a 
wire  hit.  Recoil  nuclei  are  generally  scattered  in  a  forward  direction.  If  recoil  nuclei  had  a 
significant  contribution  to  the  energy  deposition  of  these  fragmentation  events,  one  would 
expect  that  for  the  fragmentation  events  that 

1)  had  Mn  or  Cr  products  and 

2)  occurred  in  the  front  wall  of  the  detector, 

energy  deposition  would  be  enhanced  over  Fe  particles  through  the  detector.  But,  it  is  not. 
And  the  fraction  of  total  events  greater  than  190  keV  pan1  is  significantly  less  for  the  events 
where  fragmentation  occurred  as  compared  to  Fe  through  the  detector. 

Mean  detector  response  for  fragmentation  events  that  occurred  in  the  detector  wall  was 
lower  than  that  for  the  Fe  events  through  the  detector  for  all  impact  parameter  intervals  in 
the  cavity  (Fig.  4-24).  It  is  unknown  what  fraction  of  the  detector  response  is  due  to 
electron  ionization  and  recoil  nuclei  for  these  events.  But,  the  reduction  in  detector 
response  due  to  a  net  reduction  in  energy  deposition  from  electrons  is  not  fully 
compensated  by  energy  deposition  from  recoil  nuclei. 

Recommendations  for  Future  Work 

A  more  definitive  analysis  could  be  accomplished  if  there  were  a  two  PSD  pairs 
upstream  of  the  detector.  This  was  the  intention  for  these  experiments,  but  difficulties  were 
encountered  in  getting  one  upstream  PSD  pair  to  work.  With  two  pairs  of  operational 
upstream  PSDs,  impact  parameter  uncertainties  could  be  as  low  as  0.23  mm.  However, 
detailed  plots  like  those  obtained  for  the  Fe  events  through  the  detector  could  not  be 
created.  In  the  case  of  the  Fe  data,  it  was  reasonable  to  assume  that  the  response  of  the 
detector  was  homogenous  for  particles  having  die  same  impact  parameter.  But,  for 
distributions  of  detector  response  where  every  events  is  associated  with  a  nuclear 
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interaction  in  the  wall  of  the  detector,  this  assumption  would  be  unreasonable  because 
detector  response  is  influenced  by  the  site  of  the  fragmentation. 

If  two  pairs  of  upstream  and  downstream  PSDs  were  used,  nuclear  interaction  location 
could  be  determined.  Due  to  the  uncertainties  in  the  position  signals  from  PSDs,  there 
would  be  some  uncertainty  in  position,  but  the  ability  to  determine  whether  the  scattering 
event  occurred  in  a  PSD  rather  than  the  TEPC  may  be  valuable.  Discrimination  of 
interactions  in  the  PSDs  could  be  made  from  those  in  the  TEPC.  Coincidental  scattering  in 
different  detectors  would  complicate  this  analysis. 

Detector  Response  to  Ions  of  Other  Z.  but  the  Same  Velocity  as  Fe 

General 

The  stopping  power  of  high  velocity  ions  of  the  same  velocity  is  nearly  proportional  to 
the  Z2  of  the  ion  because  the  doubly  differential  secondary  electron  spectrum  is  nearly 
proportional  to  Z2  (ICRU  1996).  Subsequently,  for  ions  of  different  charge  but  the  same 
velocity,  the  spatial  pattern  of  energy  deposition  should  be  the  same  and  likewise  for  the 
fraction  of  L  captured  in  discrete  volumes. 

Separate  experiments  with  ions  of  different  charge,  but  the  same  velocity,  could  verify 
this  effect.  But,  because  experiments  may  be  conducted  on  different  days,  calibration 
uncertainties  between  individual  experiments  may  limit  the  ability  to  assess  any  differences 
that  may  exist.  The  experiments  described  here  were  not  conducted  primarily  for  the 
purpose  of  investigating  this  effect.  However,  for  the  experiment  conducted  with  the  Cu 
target,  numerous  projectile  fragmentation  products  were  created  and  could  be  separated 
based  on  the  energy  deposited  in  the  silicon  detectors.  Evaluation  of  ions  of  different 
charge  within  a  single  experiment  eliminates  calibration  uncertainties.  However,  other 
uncertainties  in  the  technique  exist. 
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Uncertainties 


1.  Trajectory 

One  important  consideration  in  evaluating  the  fragment  data  is  the  potential  for 
introduction  of  error  based  on  the  relative  trajectory  of  projectile  particles  from  the  same 
incident  projectile.  Figure  5-2  illustrates  the  trajectory  of  three  fragments  of  a  projectile. 
Because  the  evaluation  of  detector  response  was  for  particles  through  the  center  of  the 
detector  and  the  heaviest  projectile  fragment  dominates  the  PSD  position  signal,  particle  (a) 
corresponds  to  the  heaviest  fragment.  Particles  (b)  and  (c)  refer  to  lighter  fragments. 

Particle  (b)  illustrates  a  fragment  that  traverses  PSD1,  but  not  PSD2.  If  the  energy 
deposited  by  fragment  (b)  comprised  a  significant  fraction  of  the  energy  deposited  in  the 

PSD1,  the  event  would  have  been  rejected  from  the  analysis  because  the  quotient  of  AE  y2 

by  AE  Xj  would  have  been  outside  the  acceptance  region  (see  Fig.  4-26).  However,  for 
particle  (c),  it  traverses  PSD  1  and  PSD  2,  but  doesn’t  traverse  the  TEPC.  This  event 
would  not  be  rejected  through  analysis  of  the  quotient  of  AE  y2  by  AE  xp  and  would  have 

a  low  quotient  of  detector  response  (TEPC)  by  AE  y2  as  compared  to  the  case  where  both 

the  heavy  and  lighter  fragments  all  traverse  the  TEPC. 

The  bias  introduced  by  this  effect  is  dependent  on  the  combination  of  projectile 
fragments.  Figure  5-3  is  a  bar  graph  of  the  stopping  power  of  various  heavy  projectile 
fragments  and  that  of  a  single  lighter  fragment  complement.  All  values  are  normalized  to 
the  stopping  power  of  Fe.  Mn  is  not  displayed,  but  the  difference  in  ionization  between 
Mn  and  proton  vs.  Mn-alone  is  less  than  0. 16  %.  In  the  case  of  a  heavy  fragment  of  sulfur 
(Z  =  16),  the  largest  lighter  fragment  has  a  Z  of  10  (neon)  and  comprises  28  %  of  the 
combined  stopping  power  of  the  two.  If  the  sulfur  ion  traversed  the  TEPC  and  PSD2,  but 
the  neon  only  traverses  PSD2,  then  significant  bias  would  be  introduced  in  evaluation  of 
the  detector  response. 
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Figure  5-2.  Trajectories  of  Projectile  Fragments 
(from  the  Same  Incident  Projectile)  through  the 


105 


PSD2  TEPC  PSD1 


Fraction 


Figure  5-3.  Normalized  Stopping  Power  of  Heavy 
and  Light  Fragment  in  Select  Fragmentation  Events 


Z  of  Heavy  Fragment 
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Most  fragmentation  events  with  significant  reduction  in  the  Z  of  the  projectile, 
however,  will  not  have  only  one  lighter  fragment.  Rather,  these  events  are  characterized  by 
high  multiplicity  in  fragment  number,  often  numerous  doubly  and  singly  charged  nuclei. 
Appendix  C  contains  experimental  data  for  S6Fe  fragmentation. 

Calculations  were  performed  to  determine  the  angular  acceptance  of  detectors  in  the 
experimental  set-up  for  straight  projectile  trajectories  on  the  centerline  of  the  TEPC  and  for 
fragmentation  events  that  occur  in  the  Cu  target.  The  acceptance  angles  are  as  follows: 

a.  TEPC  -  0.8°, 

b.  PSD2  - 1.7°,  and 

c.  PSD1-8.30 

The  mean  emission  angle  for  a  proton  produced  in  Fe  fragmentation’s  is  2.78  +  0.04°, 
with  that  of  higher  Z  fragments  being  lower.  The  scattering  of  a  proton  from  the  projectile 
at  this  angle  in  the  experimental  set-up  would  introduce  a  lateral  separation  from  the  original 
trajectory  of  about  22  and  37  mm,  respectively  at  the  TEPC  and  PSD2.  This  proton  would 
miss  both  the  TEPC  and  PSD2.  Alternatively,  the  proton  could  be  scattered  at  an  angle  of 
1.3°,  allowing  passage  through  the  PSD2,  but  outside  the  TEPC. 

Even  for  scattering  angles  less  than  0.8°,  where  the  trajectory  of  the  lighter  fragment 
crosses  both  the  TEPC  and  PSD2,  differences  in  the  trajectory  of  the  heavy  and  lighter 
fragments  across  the  cavity  will  introduce  subtle  differences  in  detector  response  as 
compared  to  the  particles  having  the  same  trajectory. 

2.  Position  Errors 

The  PSDs  determine  the  location  of  a  particle  trajectory  based  on  the  relative  charge 
collected  on  either  side  of  the  PSD.  For  a  single  charged  particle,  position  uncertainties  are 
due  to  the  finite  resolution  of  the  PSD.  However,  when  a  combination  of  charged  particles 
traverse  a  PSD,  the  position  signal  may  be  ambiguous.  The  cases  for  three  groups  of 
projectile  fragments  are  illustrated  in  Fig.  5-4.  In  Fig.  5-4a,  the  heavy  fragment  is  Cr  (Z2  = 
576)  and  is  accompanied  by  two  protons.  For  this  case,  the  ionization  produced  by  the 
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Figure  5-4.  Combinations  of  Projectile  Fragments  from 
the  Same  Incident  Projectile  through  a  PSD  Pair 


a 

Heavy  Fragment:  Cr(Z  =  24) 
Light  Fragments:  2  -  Protons 


Heavy  Fragment:  Ca(Z  =  20) 

Light  Fragments:  Li  (Z  =  3), 
He,  &  Proton 


Heavy  Fragment:  A1  (Z  =  13) 

Light  Fragments:  O  (Z  =  8), 
2  -  He,  Proton 


protons  are  insignificant  compared  to  the  Cr  and  the  position  determination  for  the  Cr 
should  be  unaffected.  In  Fig.  5-4b,  the  heavy  fragment  is  Ca  (Z2  =  400)  and  is 
accompanied  by  Li,  He,  and  a  proton.  For  this  case,  the  ionization  produced  by  the  three 
lighter  fragments  is  3.5  %  of  that  of  the  heavy  fragment  and  a  minor  difference  in  the 
assigned  and  actual  position  of  the  heavy  particle  will  be  introduced.  In  Fig.  5-4c,  the 
heavy  fragment  is  A1  (Z2  =  169)  and  is  accompanied  by  O  (Z2  =  64),  two  He,  and  a  proton. 
The  ionization  produced  for  the  lighter  fragments  is  over  40  %  of  that  of  the  Al.  For  this 
case,  the  signal  from  the  device  will  be  ambiguous. 

Data 

Detector  response  was  analyzed  for  Fe  and  projectile  fragments  down  to  a  heavy 
fragment  Z  of  18.  Energy  deposition  in  the  TEPC  and  silicon  detectors  was  comprised  of 
that  from  all  projectile  fragments  traversing  the  active  volumes.  Energy  deposition  for 
events  that  had  the  same  heavy  fragment  may  be  vary  because  of  the  varied  number  and 
charge  of  the  lighter  fragments.  The  energy  deposited  in  PSD  y2  should  be  nearly  equal  to 
the  combined  stopping  power  of  all  projectile  fragments  traversing  it  because  the  detector  is 
thick,  an  approximation  to  CPE  conditions.  Therefore,  the  energy  deposited  in  PSD  y2 
provides  a  reasonable  estimate  of  L  for  projectile  fragments  traversing  the  TEPC.  Variance 
in  the  trajectory  among  projectile  fragments  will  bias  the  TEPC  response  to  lower  energy 
deposition  as  compared  to  that  in  PSD  y2. 

The  ratio  of  energy  deposited  in  the  detector  compared  to  AE  y2  was  consistent  among 

Fe  and  projectile  fragments  with  heavy  fragment  of  Z:  21  -  25.  For  the  fragments  with 
heavy  fragment  of  Z:  18  -  20,  the  mean  ratio  of  energy  deposited  in  the  detector  compared 
to  that  in  PSD  y2  compared  to  that  of  Fe  was  0.974  +  0.013.  For  an  event  with  this  ratio 
and  an  Ar  (Z  =  18)  heavy  fragment,  the  approximate  charge  that  traverses  PSD  y2,  but 

misses  the  TEPC  corresponds  to  Li  (Z  =  3)  or  two  ex-particles. 
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Recommendations  for  Future  Work 

A  more  definitive  analysis  could  be  accomplished  if  there  were  two  PSD  pairs  upstream 
of  the  detector,  but  downstream  of  the  Cu  target.  To  reduce  bias  observed  in  the  present 
work,  future  experiments  should  have  a  silicon  device  downstream  of  the  detector.  This 
device  should  be  place  in  close  proximity  of  the  detector  and  have  an  active  area  of  the  same 
dimensions  as  that  of  the  detector.  This  would  better  insure  that  the  projectile  fragments 
that  traverse  the  downstream  silicon  detector  also  traversed  the  detector. 
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Chapter  6:  Conclusions 

Experiments  were  performed  to  measure  the  response  of  spherical  TEPCs  to  Fe 
particles  at  1  GeV/nucleon.  The  objective  was  to  characterize  effects  caused  by  the  abrupt 
change  in  density  between  the  wall  of  the  detector  and  the  cavity  gas  used  to  simulate  tissue 
volumes  1, 2,  and  3  jim.  in  diameter.  The  experimental  arrangement  included  a  particle 
spectrometer  that  determined  the  species  and  trajectory  of  the  particles  passing  through  the 
detector. 

The  characterization  of  the  wall  effect  in  this  work  is  original  and  unique.  The 
technique  was  superior  to  previous  characterizations  because  the  analysis  could  be 
performed  for  a  spatially  uniform  particle  fluence  and  the  effect  could  be  observed  with 
respect  to  particle  trajectory. 

The  data  indicated  that  for  Fe  particles  passing  through  the  center  of  the  detector,  only 
about  75  to  80  %  of  the  energy  transferred  to  the  cavity  gas  based  on  L  was  actually 
measured  in  the  cavity  gas  ( 1  pim  simulated  diameter).  A  calculation  using  a  model  for  the 
radial  distribution  of  energy  along  a  heavy  ion  track  indicated  that  30  %  of  the  energy 

carried  by  6-rays  would  escape  the  spherical  volume.  The  difference  between  data  and  the 
model  was  a  slight  increase  in  energy  deposition  by  6-rays  produced  in  the  front  wall  of  the 

detector.  However,  for  this  trajectory,  the  enhanced  detector  response  from  6-rays 
produced  in  the  wall  of  the  thick-walled  detector  was  not  enough  to  compensate  for  the 
radial  energy  loss  from  6-rays  created  in  and  escaping  the  cavity. 

The  measured  value  of  lineal  energy  was  less  than  that  predicted  by  L  alone  for  all 
trajectories  with  impact  parameters  less  than  5.5  mm  from  the  center  of  the  cavity.  At  about 
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5.5  mm,  the  enhancement  of  detector  response  from  6-rays  created  in  the  wall  were  about 


equal  to  the  radial  energy  loss  from  6-rays  created  in  the  cavity  gas.  For  particle  trajectories 

approaching  the  cavity/wall  interface,  the  measured  energy  deposition  increased 
dramatically.  Events  with  the  largest  energy  deposition  (i.e.,  y  >  250  keV  pim'1)  were 
associated  with  a  very  narrow  band  of  impact  parameters  located  at  the  cavity/wall 
interface.  An  analysis  of  nuclear  cross-sections  and  variations  in  the  trajectory  through  the 
detector  indicated  that  the  very  large  energy  deposition  events  were  due  to  enhanced  energy 
deposition  by  soft  collision  electrons  in  the  cavity  rather  than  that  of  recoil  nuclei. 

The  data  for  energy  deposition  as  a  function  of  impact  parameter  for  the  thick-walled 
TEPC  were  integrated  to  determine  the  absorbed  dose  measured  in  the  simulated  volume  of 
tissue.  For  a  /<- random  Fe  particles,  the  absorbed  dose  corresponded  to  particles  with  an  L 
=  154  keV  pim~l  in  1  and  2  fim  diameters  of  tissue.  A  value  of  150  keV  pm1  was  obtained 
for  the  3  pirn  simulated  diameter.  The  minor  difference  in  the  experimental  data  from  theory 
was  believed  to  arise  from  calibration  or  position  uncertainties.  Thus,  the  enhancement  of 

6-rays  to  the  detector  response  was  sufficient  to  create  CPE  in  the  cavity. 

Walled  spherical  TEPCs  can  provide  accurate  determination  of  absorbed  dose, 
however,  the  distribution  of  energy  deposition  on  an  event  by  event  basis  does  not 
represent  that  that  would  be  observed  in  homogenous  media. 

Q  was  determined  for  a  /^-random  irradiation  to  Fe  particles  where  the  detector 
simulated  a  1  ptm  diameter  site.  Close  approximation  of  the  actual  Q  was  obtained  under 
the  assumption  that  y  =  L  for  both  ICRP  26  and  60  recommended  relationships  between  Q 
and  L.  Similar  results  were  obtained  when  the  data  was  modified  for  ions  of  the  same 
velocity,  but  lower  Z. 

For  projectile  fragments  with  a  heavy  fragment  of  Z  between  21  and  25  traversing  the 
center  of  the  detector,  the  fraction  of  L  absorbed  in  the  detector  was  the  same  as  that  of  Fe. 
For  heavy  fragments  of  Z  between  18  and  20,  the  fraction  of  L  absorbed  in  the  detector 
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was  about  97  %  of  that  of  Fe.  Reduced  energy  deposition  was  due  to  inadequacies  in  the 
experimental  design. 

For  charge-changing  nuclear  interactions  in  the  wall  of  the  detector,  detector  response  is 
lower  than  that  of  Fe  particles  through  the  detector.  Energy  deposition  from  recoil  nuclei 
do  not  appear  to  compensate  for  decreased  energy  deposition  from  electrons. 
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Appendix  A:  Chatterjee  Track  Structure  Model 

The  Chatteijee  model  (1973, 1976,  and  1993)  defines  energy  density  around  the 
trajectory  of  heavy  charged  particles.  The  energy  deposition  is  defined  for  two  distinct 
regions:  the  core  and  penumbra.  The  energy  deposited  in  the  core  is  defined  by  equation 
A-l: 


L[2  +  21n(%)] 
2jtrc2  [1  +  21n(%)]  ’ 


(A-l) 


where  rp,  rc,  L  are  the  penumbra  radius,  core  radius,  and  the  stopping  power,  respectively. 
The  energy  deposited  in  the  penumbra  is  defined  by  equation  A-2: 


L 

Ppenumbra  ~  ^  r — 7 - »  (A-2) 

2nrz  [1  +  21n(%)] 

where  r  is  the  radial  distance  from  the  particle  trajectory.  The  plot  in  Fig.  A-l  provides 
core  and  penumbra  radii. 

Figure  A-2  provides  a  plot  of  dose  (keV  ptm'3)  versus  radial  distance  from  the  trajectory 
of  a  1  GeV/nucleon  56Fe  in  homogenous  unit  density  media.  Although  it  is  agreed  that 
there  is  no  clear  demarcation  between  a  core  region  and  the  penumbra  portion  of  the  track, 
it  is  treated  as  such  in  the  model.  For  both  regions  of  the  track,  it  is  assumed  that  energy 
deposition  is  homogenous.  From  the  graph,  it  is  apparent  that  beyond  the  core  region, 
dose  is  a  function  of  the  reciprocal  of  the  radial  distance  squared. 
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Figure  A-l.  Core  and  Penumbra  Radii  for  Chatterjee  Model  (1993) 


Figure  A-2.  Dose  vs.  Radial  Distance  from  Trajectory  For  1  GeV/nucleon 
S6Fe  Particle,  Chatteijee  Model  (1993) 


The  model  can  be  used  to  calculate  energy  absorbed  is  volumes  of  any  size  in  the 
vicinity  of  heavy  ion  trajectories  as  illustrated  in  Figure  A-3  for  a  sphere.  In  the  Figure,  the 
trajectory  is  perpendicular  to  the  page,  dense  ionization  in  vicinity  of  the  core  is  simulated 

by  black  and  dark  shading,  while  6-rays  are  extending  radially  from  the  trajectory.  Three 

trajectories  are  illustrated:  one  where  the  trajectory  is  penetrating  the  sphere  in  the  center 


Figure  A-3.  Trajectories  of  Heavy  Ion  in  Vicinity  of  Spherical  Volume. 


where  the  chord  length  through  the  volume  is  maximum  (a),  another  where  the  trajectory  is 
outside  of  the  sphere  (b),  and  the  last  where  the  trajectory  is  through  the  sphere,  but  near 
the  edge  (c).  For  volumes  with  diameters  on  the  order  of  microns  that  are  applicable  to 
microdosimetry  (and  illustrated  in  Fig.  A-3),  virtually  all  soft  collision  electrons  produced 
along  the  ion  trajectory  in  the  sensitive  volume  will  be  absorbed.  For  trajectories  where 
the  ion  does  not  intersect  the  volume  (Fig.  A-3b),  soft  collision  electrons  will  not  contribute 

to  absorbed  energy  in  the  volume  and  energy  deposition  will  be  from  energetic  6-rays 
generated  outside  the  sensitive  volume.  For  the  model,  provided  the  radial  distance  of  the 
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ion  trajectory  to  the  volume  is  less  the  penumbra  radius,  ionization  from  6-rays  will 

contribute  to  the  absorbed  energy  in  the  volume. 

Figure  A-4  illustrates  the  approach  taken  to  calculate  the  energy  absorbed  in  sensitive 
volumes  in  the  vicinity  of  ion  trajectories.  The  volume  of  the  detector  was  divided  into 


Figure  A-4.  Voxelized  Detector  Volume  in  Vicinity  of  Ion  Trajectory. 


many  smaller  discrete  volumes  (voxels)  as  illustrated  in  the  Figure.  The  energy  deposited 
in  each  volume  was  the  product  of  the  radial  dose  and  the  volume  of  the  voxel.  The  voxel 
volume  was  the  product  of  the  cross-sectional  area  of  each  voxel  and  the  voxel  length. 
Voxel  length  was  calculated  from  geometrical  dimensions  of  a  sphere.  For  calculations, 
detector  volumes  were  broken  into  4280  voxels  of  cross-sectional  voxel  area  dependent  on 
the  radius  of  the  sphere  simulated.  Absorbed  energy  in  the  vicinity  of  1  GeV/amu  S6Fe  was 
calculated  for  spheres  of  varying  radii.  Figure  A-5  provides  a  summary  of  the  calculations 
for  a  0.5  jm l  radius  sphere  in  terms  of  the  mean  lineal  energy. 
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The  fraction  of  the  L  absorbed  in  volumes  directly  traversed  by  ions  is  a  parameter 
of  interest.  Chatteijee’s  model  ( 1993)  provides  a  formula  that  relates  the  fraction  of  L 
absorbed  in  an  infinite  cylinder  around  the  trajectory  of  an  ion  for  cylinders  of  various 
radii.  The  formula  is  a  simple  integration  of  Eqs.  A-l  and  -2  from  the  center  of  the  cylinder 
to  the  radius  of  interest.  Figure  A-6  provides  a  plot  of  the  fraction  of  L  absorbed  in  infinite 
cylindrical  and  bound  spherical  volumes  for  1  GeV/nucleon  ions  at  various  radii.  As 
expected,  the  infinite  cylinders  will  have  higher  absorbed  fraction  of  L  than  the  spheres. 
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Fraction  of  L 


Appendix  B:  Calibration 


General 

The  objective  of  microdosimetry  is  to  measure  the  energy  deposited  in  a  sensitive 
volume,  while  TEPCs  measures  the  charge  collected  per  event  that  is  proportional  to  the 
number  of  ion  pairs  formed  within  the  simulated  volume  (cavity).  The  two  parameters  are 
related  to  each  other  by  W,  the  mean  energy  expended  per  ion  pair  formed.  W-values  vary 
considerably  among  radiation  types,  proportional  counter  fill  gases,  energy  of  the  incident 
radiation,  and  the  energy  deposited  by  the  radiation. 

Energy  calibrations  must  be  carried  out  by  irradiating  the  counter  with  a  radiation  that 

deposits  a  known  amount  of  energy.  Generally,  a-particle,  neutron,  or  x-ray  sources  are 

used.  K-line  aluminum  (1.49  keV),  55Fe,  and  37Ar  x-ray  sources  are  used,  with  the 
aluminum  and  argon  more  appropriate  for  counters  that  simulate  diameters  less  than  a 
micron.  Additionally,  x-ray  calibration  sources  are  more  appropriate  for  low-LET  radiation 
measurements  since  the  energy  range  and  W-values  of  the  two  are  more  closely  matched. 
For  x-ray  sources,  the  range  of  the  secondaries  are  considerably  lower  than  a  few  microns 
(ICRU  1983),  dimensions  typical  of  traditional  microdosimeters.  Thus,  the  geometry  does 
not  limit  the  ionization  produced.  For  walled  proportional  counters  constructed  of  A-150 
plastic,  over  time  and  under  continuous  vacuum,  the  plastic  shrinks  and  becomes 
dimensionally  unstable  leading  to  errors  in  the  actual  size  of  the  cavity  (Kliauga  1990). 
These  errors  can  invalidate  a  calibration.  Another  problem  with  some  x-ray  sources  is  that 
they  do  not  sample  the  entire  detector  volume;  37Ar  sources  avoid  this  problem  because  they 

are  mixed  with  the  counting  gas  (Anachkova  etal.  1991).  Built-in  a-particle  calibration 
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sources  of  241Am  and  244Cm  are  commonly  used  in  commercially  available  microdosimeters 
(Schrewe  et al.  1988).  The  sources  are  commonly  mounted  on  the  outside  of  the  wall  of 
the  detector.  A  small  hole  drilled  in  the  wall  allows  penetration  and  collimation  of  the 

a-particles.  A  dominant  peak  in  the  spectrum  is  caused  by  a-particles  that  cross  the  cavity 

along  the  full  chord  length  for  spherical  counters  and  normally  to  the  height  for  cylindrical 
counters.  The  spherical  TEPCs  used  in  these  experiments  have  built-in  244Cm  sources  that 
are  gravity  controlled.  External  fast  neutron  sources  are  another  method  of  calibration  and 

use  recoil  protons  as  the  radiation  source.  In  contrast  to  built-in  a-calibration  sources, 

these  sources  do  not  provide  recoil  protons  of  a  single  energy  or  single  pathlength  in  the 
cavity.  Calibration  is  made  by  evaluation  of  the  proton-edge  of  the  collected  spectrum  -  a 
position  in  the  spectrum  where  protons  have  both  maximum  stopping  power  and  maximum 
chord  length.  Calibration  of  cylindrical  counters  with  neutron  sources  requires  special 
procedures  due  to  the  chord  length  distribution  (Anachkova  etcd.  1994). 

W-values 

ICRU  Report  31  (1979)  provided  a  comprehensive  review  of  W-values.  ICRU  defines 
W  as  the  quotient  of  E  by  N,  where  N  is  the  mean  number  of  ion  pairs  formed  when  the 
initial  kinetic  energy,  E,  of  a  charged  particle  is  completely  dissipated  in  the  gas: 


The  differential  value,  w,  of  the  mean  energy  required  to  produce  an  ion  pair  is  the  quotient 
of  dE  by  dN,  where  dE  is  the  mean  energy  lost  by  a  charged  particle  of  energy  E  in 
traversing  an  absorber  of  thickness  dx,  and  dN  is  the  mean  number  of  ion  pairs  produced 
when  dE  is  completely  dissipated  in  the  gas: 


w 


dE 

dN 


(B-2) 
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The  relationship  between  W  and  w  is  defined  by: 


W(E)  = 


E 

jf  de/w(e)’ 


(B-3) 


where  I  is  the  lowest  ionization  potential  of  the  medium  and  e  is  the  instantaneous  energy  of 

the  particle  while  slowing  down.  Numerous  measurements  of  W  have  been  performed, 
many  of  which  are  summarized  in  ICRU  Report  31,  but  very  few  have  been  made  on  w 
(Bichsel  1972). 

ICRU  (1979)  and  Goodman  and  Coyne  (1980)  have  summarized  W  values  for 
protons,  a-particles,  and  heavy  ions  in  tissue-equivalent  (TE)  gas.  Goodman  and  Coyne 

performed  least-squares  evaluations  of  experimental  data  Empirical  functions  for  protons 
and  a-particles  from  their  work  are: 


W(protons)  =  29.28  ( In  E  )'2  +  29.99  (eV  ion  pair'1)  and  (B-4) 

W(a-particles)  =  49.96  ( In  E  )'02418  (eV  ion  pair'1) ,  (B-5) 

where  E  is  initial  particle  energy  in  keV/nucleon.  Figure  B-l  provides  a  plot  of  the 
Goodman  and  Coyne  empirical  functions  for  protons,  a-particles,  and  oxygen  ion.  Three 

noticeable  features  of  the  functions  are  apparent.  First,  W  values  for  all  three  particles 
increase  significantly  for  low  initial  energy  particles  and  is  due  to  a  larger  fraction  of  the 
energy  transferred  to  atoms  or  molecules  in  the  form  of  excitation  rather  than  ionization. 
Second,  for  low  energy  particles,  W  values  are  higher  for  particles  of  higher  Z-number. 
And  lastly,  at  higher  velocities,  W  values  become  relatively  independent  of  energy.  Myers 
( 1968)  noted  that  W  values,  in  general,  should  exhibit  little  dependence  on  energy  above: 

Eo=|mvJ,  (B-6) 
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where  M  is  the  projectile  mass  and  v0  is  the  velocity  of  an  electron  in  the  first  Bohr  orbit. 


Figure  B-l.  Empirical  Functions  of  W  for  Protons,  a-Particles, 
and  Oxygen  Ions  (Goodman  and  Coyne  1980). 


Energy  (keV/nucIeon) 

Very  little  experimental  work  on  W-  and  w -values  determination  have  been 
accomplished  for  high-energy  heavy  ions.  Goodman  and  Colvett  ( 1977)  reported  that  the 

value  of  w(e)  in  TE  gas  for  429  MeV/amu  Ar18+  ions  to  be  similar  to  that  of  137Cs  photons, 

under  the  assumption  that  the  w(e)  of  60Co,  29.2  ±  0.6  eV  ion  pair1  from  Leonard  and 

Boring  (1973),  is  the  same  as  that  for  137Cs.  This  value  is  reasonably  close  to  the  W  of 
3 1.0  eV  ion  pair'1  for  5.3  MeV  using  the  Goodman  and  Coyne  empirical  formula.  Thomas 

etal.  (1980)  measured  w(e)  values  for  250  MeV/nucleon  C6+,  375  MeV/nucleon  Ne10+,  and 

479  MeV/nucleon  Ar18+  ions  in  nitrogen  gas.  Respectively,  the  reported  values  along  with 
percent  fractional  errors  (%  FE)  were:  36.4  ±  0.6  (%  FE  =  1.7  %),  35.4  ±  0.8  (%FE  = 

2.3  %),  34.7  +  0.5  (%  FE  =  1.4  %)  eV  ion  pair'1.  From  the  experimental  data,  the  authors 
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could  not  conclude  that  there  was  a  Z  dependence  of  w(e)  nor  any  significant  difference 

between  the  three  values.  Schimmerling  etal.  (1982)  measured  w(e)  for  Ne10+,  Al13+,  Si14+ 

in  nitrogen  and  P-7  gas  over  particle  energy  range  of  100  -  700  MeV/nucleon.  The  mean 
value  reported  for  nitrogen  gas  was  33.0  +  2.3  eV  ion  pair1,  with  %  FE  of  7.0  %.  The 
values  of  Thomas  etal.  are  close  to  the  recommended  W  value,  36.39  +  0.23  eV  ion  pair1, 

reported  in  ICRU  Report  31  for  a-particles  of  approximate  energy  5.3  MeV  in  nitrogen. 

The  value  reported  by  Schimmerling  etal.  is  lower  than  the  ICRU  value,  but  may  not  be 
different  considering  the  relatively  large  %  FE. 

Calibration  Comparisons  with  Far-West  Technologies  1/2”  LET  TEPC 

Calibration  was  performed  by  three  separate  techniques  for  the  thin-walled  TEPC  used 
in  these  experiments: 

1.  single-event  peak  of  unscattered  a-particles  from  the  built-in  244Cm  source, 

2.  proton-edge  from  irradiation  with  239Pu:Be  unmoderated  neutron  source,  and 

3.  a-particle-edge  from  a-particles  that  are  scattered  in  the  borehole  of  the  detector 
wall  as  described  by  Schrewe  et  al.  (1988). 

The  first  and  third  technique  were  performed  for  the  thick-walled  detector . 

In  the  comparison,  stopping  power  values  for  protons  and  a-particles  were  taken  from 

ICRU  Report  49  (1993)  for  propane-based  TE  gas.  Mean  energy  deposited  by  particles 
was  calculated  with  these  stopping  power  values  and  the  continuous  stopping  distance 
approximation  (CSDA).  Because  the  particles  used  for  the  calibrations  do  not  completely 
stop  in  the  cavity  of  the  detector,  W-values  were  determined  from  the  empirical  formulas  of 
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Goodman  and  Coyne  and  Equation  B-2,  where  dE  is  difference  in  the  initial,  E,,  and  final, 
Ej,  energy  of  the  particle,  and 

dN  =  (Ej/Wj)  -  (Ef  / Wf  ),  (B-7) 

where  W;  and  Wf  are  the  W- values  respectively  for  particles  of  energy  equal  to  the  initial 
and  final  energy  of  the  particle.  The  classical  definition  of  W-value  infers  an  integration  of 

w(e)  from  the  initial  particle  energy  to  that  of  the  lowest  ionization  potential.  To  avoid 

confusion,  W-values  calculated  in  this  manner  will  be  referred  to  as  integral  w -values. 
Initial  a-particle  energy  from  the  244Cm  source  was  assumed  to  be  5.789  MeV,  the 

frequency-weighted  mean  of  the  energy  of  the  two  alpha  particles  emitted.  Also,  separate 
calculations  of  energy  transferred  were  performed  for  the  single-event  peak  of  unscattered 

a-particles  from  the  built-in  244Cm  source  for  the  thin-  and  thick-walled  detectors.  This  is 
necessary  since  the  initial  and  final  energy  of  the  a-particle  will  be  different  for  the  two 

detectors  because  there  is  a  difference  in  the  energy  lost  by  a-particles  in  the  wall  region. 

Table  B-l  summarizes  factors  used  in  the  calibrations  and  comparison  among  techniques. 
From  the  energy  transferred  values  listed  in  the  Table,  it  is  apparent  that  there  are 

significant  differences  dependent  on  simulated  diameter.  For  the  a-peak  values,  there  is  a 

1.7  and  1.9  percent  difference  between  the  energy  transferred  for  the  0.5  and  3  pa n 
simulated  diameters,  respectively  for  the  thin  and  thick  walled  cavities.  However,  for  the 

a-  and  proton-edge  values,  the  differences  are  greater,  3.3  %  and  15.8  %,  respectively. 

In  comparison  of  the  integral  w-values  from  0.5  to  3  ptm  for  each  of  the  techniques,  the 
maximum  difference  was  0.7  %  and  deemed  trivial.  However,  the  difference  between  the 

integral  w-values  of  the  a-edge  technique  and  the  other  two  is  as  high  as  7.4  %.  This 
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difference  is  very  close  to  that  reported  by  Schrewe  etal.  (1988)  for  a  similar  comparison 
of  calibration  techniques. 


Table  B-l.  Calibration  Factors 


Technique 

Simulated 
Diameter  (/mi) 

Energy  Transferred 
(keV  //m1) 

Initial  Particle 
Energy  (MeV) 

Final  Particle 
Energy  (MeV) 

Integral  w- Value 
(eV  ion  pair'1) 

a -peak 

0.5 

83.79 

5.785 

5.743 

29.93 

thin 

1 

84.06 

5.781 

5.697 

29.93 

walled 

2 

84.62 

5.772 

5.603 

29.94 

detector 

3 

85.19 

5.764 

5.508 

29.95 

a -peak 

0.5 

83.83 

5.781 

5.739 

29.93 

thick 

1 

84.15 

5.772 

5.688 

29.93 

walled 

2 

84.80 

5.756 

5.568 

29.94 

detector 

3 

85.48 

5.739 

5.482 

29.95 

proton 

0.5 

100.8 

0.107 

0.057 

30.83 

edge 

1 

97.75 

0.139 

0.041 

30.82 

(max. 

2 

91.21 

0.206 

0.024 

30.79 

energy) 

3 

86.06 

0.270 

0.015 

30.76 

alpha 

0.5 

259.3 

0.694 

0.564 

32.22 

edge 

1 

257.8 

0.765 

0.507 

32.21 

(max. 

2 

255.2 

0.926 

0.416 

32.17 

energy) 

3 

250.9 

1.111 

0.357 

32.10 

Comparison  of  Thin-Walled  Detector 

The  spectrum  collected  for  the  internal  a-calibration  source  is  displayed  in  Figure  B-2. 


The  prominent  peak  near  channel  60  corresponds  to  unscattered  a-particles,  while  the  point 
of  maximum  energy  deposition  of  scattered  a-particles  occurs  near  channel  150  and  is 
called  the  a-edge.  Events  with  magnitude  significantly  higher  than  that  of  the  a-edge  are 
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few  in  number  and  are  likely  attributable  to  a-particle  coincidences.  The  peak  to  the  left  of 

that  for  unscattered  a-particles,  has  a  peak  frequency  about  5  %  of  the  peak  for  unscattered 

a-particles.  The  lower  energy  deposition  peak  has  a  mean  energy  deposition  about  half 

that  of  the  prominent  peak  and  is  due  to  a-particles  that  are  absorbed  by  the  helix  or  anode 
wires  after  passing  part  of  the  way  through  the  cavity. 

Figure  B-2.  Internal  a-Particle  Calibration  Source  Spectrum 
(Serial  Number  =  1373,  HV  =  600  volts.  Simulated  Diameter  =  1  pirn) 
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The  peak  corresponding  to  unscattered  a-particles  appears  asymmetric  and  would  result 

in  a  difference  between  the  mean  and  most  probable  value  (Schrewe  etal.  1988).  Schrewe 
et  al.  (1988)  calculated  the  distribution  functions  for  various  simulated  diameters  of 
spherical  TEPCs  and  determined  that  energy  loss  straggling  follows  a  Vavilov  distribution 
for  diameters  up  to  2  ^m  and  that  the  relative  difference  in  the  mean  and  most  probable 
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energy  deposition  are  small  at  1  and  0.6  %,  respectively  for  1  and  2  ym.  Because  the 
difference  is  small  for  simulated  diameters  greater  1  ym,  and  1  ym  was  the  smallest 
diameter  used  in  our  experiments,  energy  deposition  corresponding  to  the  most  probable 
value  will  be  assumed  to  be  equal  to  the  mean.  To  evaluate  the  magnitudes  corresponding 

to  the  most  probable  energy  and  a-edge,  the  first  derivative  of  the  frequency  distribution 

was  calculated  by  a  numerical  technique  of  Savitzky  and  Golay  (1964). 

Figure  B-3  is  the  first  derivative  of  the  calibration  data  of  Fig.  B-2,  truncated  for  data 
channels  in  the  vicinity  of  the  most  probable  value.  Interpolation  of  the  data  points 
provides  an  approximate  “zero”  crossing-point  at  channel  56% . 


Figure  B-3.  first  Derivative  of  Data  from  Fig.  B-2  (Most  Probable  Value  Region). 


Figure  B-4  is  the  first  derivative  of  the  same  data  in  the  vicinity  of  the  a-edge.  Precise 

evaluation  of  the  a-edge  is  complicated  and  would  require  a  detailed  simulation  of  the 
energy  loss  and  straggling  outside  and  inside  of  the  detector  according  to  Schrewe  etal. 
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(1988).  From  inspection  of  the  plot,  there  is  a  broad  minimum  in  the  distribution  with  the 
center  near  channel  157. 


Figure  B-4.  First  Derivative  of  Data  from  Fig.  B-2  (a-Edge  Region) 


Channel  Number 

For  neutron  spectra  collection,  a-particle  calibrations  were  performed  before  and  after 
the  neutron  spectrum  was  collected  to  insure  constant  gas  pressure  in  the  TEPC  existed 
throughout  the  measurements.  Before  and  after  a-calibration  spectra  are  displayed  in  Fig. 

B-5  truncated  in  the  region  of  the  most  probable  value.  Figure  B-6  contains  the  first 
derivative  spectra  of  the  data  in  Fig.  B-5.  Most  probable  value  constancy,  ~  channel  80.6, 
indicates  constant  gas  pressure  throughout  the  measurements. 

Figure  B-7  is  a  plot  of  the  neutron  spectrum  collected  over  a  14-hour  period.  Apparent 
from  the  spectrum  is  the  proton-edge  near  channel  100.  Figure  B-8  contains  the  first 
derivative  spectrum  of  the  neutron  data  truncated  for  the  channels  in  the  vicinity  of  the 

proton-edge.  As  used  for  evaluation  of  the  a-edge,  the  local  minimum  in  the  derivative  is 
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Figure  B-5.  Before  and  After  a-Calibration  Spectra  for  Neutron  Measurements. 
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Figure  B-6.  First  Derivative  of  Data  from  Fig.  B-5. 
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Figure  B-7.  Neutron  Spectrum  (Serial  Number  =  1373,  HV  =  600  volts. 
Simulated  Diameter  =  1  /<m) 
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the  point  where  the  frequency  spectrum  has  its  steepest  negative  slope  and  is  designated  the 
proton-edge.  For  this  data  set,  channel  98  possesses  the  lowest  value  in  the  derivative 
spectrum  in  this  region.  From  inspection  of  Fig.  B-7,  this  criterion  seems  reasonable. 

To  provide  accurate  comparison  of  the  three  techniques,  a  research  pulser  calibration  of 
the  multichannel  analyzer  was  performed.  Table  B-2  provides  pulser  settings  for  analyzer 
channels  5  through  254. 

The  relative  location  of  the  three  features  are  compared.  Pulser  settings  are  related  to 
relative  energy  deposited  from  Table  2. 

1.  The  ratio  of  the  pulser  settings  determines  the  relative  ratio  of  ion  pairs  created 
by  two  techniques  being  compared. 

2.  The  ratio  of  the  number  of  ion  pairs  expected  between  two  techniques  being 
compared  is  calculated  from  energy  transferred  values  and  integral  w -values  from  Table  1. 
Finally,  the  two  ratios  from  1.  and  2.  above  are  compared.  Equation  B-8  is  a  summary  of 

the  calculations  for  comparison  of  the  a-peak  and  proton-edge: 

Ratio  =  -  RS'(test)  *  Et(«-Peak)  *  w(test) 

P.S.(a-peak)  Et(test)  w(a  -  peak)  ’ 

where  P.S.  is  the  relative  pulser  setting  and  Et  is  the  energy  transferred.  For  comparison  of 
the  a-edge  to  a-peak,  the  ratio  is: 


Ratio  = 


605.88  *  84.06  #  32.21 
207.66  257.8  29.93 


1.024. 


(B-9) 


Equation  B-10  summarizes  the  calculation  of  the  ratio  of  the  expected  number  of  ion  pairs 
between  the  proton-edge  to  those  produced  by  the  mean  of  the  a-peak. 


367.44  *  84.06  *  30.82 
325.38  97.75  29.93 


Ratio  = 


(B-10) 


Table  2.  Relationship  of  Research  Pulser  Settings  to  Multichannel 
Analyzer  Channel  Number  for  Conversion  Gain  =  256 
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Channel  Pu 
Number  Set 


From  the  ratios,  it  is  apparent  that,  within  the  uncertainties  of: 

1.  selection  of  an  appropriate  position  of  the  a-  and  proton-edges, 

2.  integral  w -values,  and 

3.  energy  transferred  calculations, 

that  the  three  techniques  provide  similar  results  for  this  detector. 

Comparison  of  Thick- Walled  Detector 

The  spectrum  collected  for  the  internal  a-calibration  source  is  displayed  in  Fig.  B-9. 

The  prominent  peak  of  the  spectrum  is  similar  to  that  of  the  thin-walled  detector.  But, 
some  notable  differences  exist.  The  peak  to  the  left  of  the  prominent  peak,  that  is  caused 

by  absorption  of  a-particles  in  the  helix  and  anode  wires,  is  of  greater  magnitude  than  that 

of  the  thin-walled  detector.  This  is  due  to  a  higher  degree  of  collimation  and  a  higher 
fraction  of  unscattered  particles  being  absorbed  by  the  wires.  Also,  the  scattered 
contribution  to  the  spectrum  is  greater  with  a  more  prominent  plateau  than  for  that  of  the 
thin- walled  detector.  Figure  B-10  is  the  first  derivative  of  the  spectrum  contained  in  Fig. 
B-9  truncated  for  the  region  of  the  most  probable  value  that  is  located  at  channel  65.  Figure 

B-ll  is  the  first  derivative  of  the  spectrum  in  the  vicinity  of  the  a-edge.  The  location  of  the 

a-edge  (channel  181)  is  more  apparent  than  was  the  case  of  the  thin-walled  detector  that 

had  a  broader  distribution  at  the  minimum.  This  is  reasonable  since  the  thick-walled 
detector  had  a  higher  proportion  of  scattered  particles.  Using  Equation  B-8,  the  ratio  of  the 

expected  number  of  ion  pairs  produced  by  the  a-edge  and  a-peak  is  calculated: 

Z2LZ 1  ,  8415  ,  3i21_ 

240.14  257.8  29.93 


Ratio  = 


(B-ll) 


d  (Counts)  Counts 


Figure  B-9.  Internal  a-Particle  Calibration  Source  Spectrum 
(Serial  Number  =  1376,  HV  =  600  volts.  Simulated  Diameter  =  1  pm) 
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'igure  B-10.  First  Derivative  of  Data  from  Fig.  B-9  (Most  Probable  Value  Region) 
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Figure  B- 1 1  First  Derivative  of  Data  from  Fig.  B-9  (a-Edge  Region) 


Thus,  for  the  thick-walled  detector,  as  was  the  case  for  the  thin,  there  is  close  agreement 
between  the  two  techniques. 

Calibration  for  S6Ions 

Based  on  the  comparison  of  calibration  techniques,  best  available  literature  on  up¬ 
values,  it  is  reasonable  to  apply  the  most  probable  value  of  the  internal  a-calibration  source 
spectra  directly  to  S6Fe  ions  without  any  correction  for  vv-values. 
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Appendix  C:  Nuclear  Interactions 


General 

Considerable  interest  and  effort  to  model  nuclear  interaction  cross-sections  of  high- 
velocity  ions  have  been  performed  by  NASA  and  others  over  the  last  40  years.  Details  on 
current  models  used  by  NASA  are  in  NASA  Reference  Publication  1257  (Wilson  etal 
1991).  Heavy  charged  particles  traversing  matter  lose  energy  primarily  through  ionization 
and  excitation  of  atoms,  except  at  very  low  velocities  where  nuclear  interactions  become  an 
important  form  of  energy  loss.  While  nuclear  interactions  are  responsible  for  relatively 
small  energy  transfer  in  high-velocity  ions,  the  interactions  alters  energy  deposition  at  the 
interaction  site  and  the  characteristics  of  the  radiation  field  from  fragmentation  interactions 
(Wilson  etal.  1991).  Nuclear  interactions  between  projectile  and  target  nuclei  involve 
electromagnetic  forces,  also  referred  to  as  the  Coulombic  force.  Nuclear  forces  are  also 
involved  when  the  impact  parameter  between  the  two  nuclei  is  less  than  the  sum  of  the 
respective  nuclear  radii.  In  interactions  of  this  type,  the  nuclear  forces  typically  dominate 
the  interaction.  For  nuclear  interactions  that  occur  at  impact  parameters  greater  than  the 
sum  of  the  nuclear  radii,  electromagnetic  interactions  are  only  important.  These  interactions 
contribute  to  a  phenomenon  called  multiple  Coulomb  scattering,  have  a  negligible 
contribution  to  energy  loss,  have  negligible  interaction  cross-sections  as  compared  to  that 
of  collisions  (Wilson  etal.  1991),  and  are  ignored  in  discussion  of  nuclear  cross-sections 

in  the  remainder  of  this  appendix.  The  total  nuclear  interaction  cross-section,  o^.,  is  the 
sum  of  the  elastic  interaction  cross-section,  as,  and  the  absorption  (inelastic)  interaction 
cross-section,  aABS. 
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Elastic  Scattering 


Elastic  scattering  is  a  nuclear  interaction  process  where  there  is  conservation  of 
momentum  and  kinetic  energy  and  the  two  nuclei  remain  intact.  Figure  C-l  provides  an 
example  of  elastic  scattering.  The  example  illustrates  an  interaction  between  a  high-velocity 
56Fe  projectile  nucleus  and  a  target  12C  nucleus.  In  the  interaction,  the  path  and  velocity  of 
the  iron  nucleus  remains  relatively  unchanged  and  neither  of  the  nuclei  have  a  change  in 
nucleon  number,  A  or  charge,  Z.  Particles  from  these  interactions  are  scattered  primarily  in 
the  forward  direction  (ICRU  1993). 

Figure  C-l.  Elastic  Scattering. 

Before  After 


Wilson  and  Costner  (1975)  applied  the  coherent  wave  as  an  approximation  to  the  elastic 
channel  and  found  good  agreement  with  available  experimental  data.  NASA  uses  these 
approximations  in  their  nuclear  transport  calculations  (Wilson  etal.  1991).  For 
measurements  collected  in  this  experiment,  production  of  scattered  nuclei  in  the  wall  of  the 
TEPC  could  influence  the  response  in  the  cavity  if  these  particles  have  sufficient  penetration 
range.  Equation  B-l  defines  the  maximum  energy  transfer  from  a  projectile  to  target 
particle  at  rest  (Turner  1992). 
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2  max  — 


(C-l) 


2y  2  Mi  V2 


where  Qmax  is  the  maximum  kinetic  energy  transfer,  is  the  rest  mass  of  the  target  nuclei, 

is  rest  mass  of  the  projectile,  y  =  V  ^1-p2  ,  P  =  V/c,  and  c  is  the  speed  of  light. 

For  1  GeV/amu  S6Fe  ions,  the  maximum  kinetic  energy  transfer  to  is  5800  and  3200 
MeV/nucleon  to  proton  and  carbon  nuclei,  respectively,  the  two  primary  components  of  the 
TEPC  wall.  Figure  C-2  is  a  plot  of  particle  penetration  range  in  water  versus  kinetic  energy 
based  on  ICRU  49  (ICRU  1993)  and  Equation  4.25  of  Turner  (1992).  From  the  plot,  it  is 
apparent  that  carbon  or  hydrogen  recoil  nuclei  of  maximum  energy  from  1  GeV/nucleon  Fe 
ions  will  have  a  penetration  range  in  excess  of  106  ja m  ( 1  m).  But,  the  majority  of  recoils 
will  have  considerably  less  energy,  where  the  penetration  range  is  greatly  reduced. 

Absorption  (Inelastic)  Interactions 

In  general,  products  of  absorption  nuclear  interactions,  are  fragments  of  the  projectile 
nucleus,  fragments  of  the  target  nucleus,  and  nucleons  (or  groups  of  them  like  a-particles) 

produced  in  the  interaction  region  where  the  two  nuclei  overlap  in  the  collision  (Curtis  and 
Schimmerling  1977).  Projectile  fragments  emerge  from  the  interaction  in  a  direction  very 
close  to  the  initial  particle  direction  and  with  little  or  no  change  in  velocity.  Target 
fragments,  however,  have  a  broad  angular  distribution,  low-energies,  and  high  multiplicity 
in  particle  number.  Central  or  “head  on”  collisions  are  unique  where  products  are  formed 
from  projectile  and  target  nucleons,  and  have  intermediate  particle  kinetic  energies  between 
that  described  above  for  projectile  and  target  fragments.  According  to  Curtis  and 
Schimmerling  (1977),  projectile  fragmentation,  target  fragmentation,  and  central  collisions 
have  the  following  additional  characteristics. 
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Range  in  Water  (/*m) 


Figure  C-2.  Range  of  Proton  and  Carbon  Nuclei  in  Water  (unit  density). 
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Projectile  Fragmentation 

1.  The  reactions  resulting  in  projectile  fragments  are  due  to  grazing  collisions. 

2.  Total  reaction  cross-sections  are  generally  independent  of  energy  for  projectiles  of 


energy  greater  than  300  MeV/nucleon. 
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3.  Fragment  velocity  and  trajectory  is  close  to  that  incident  projectile. 

4.  The  probability  of  producing  a  fragment  is  related  to  the  target  mass  only  through  a 
geometric  scale  factor. 

5.  The  interaction  cross-sections  have  been  approximated  by  considering  the  overlap 
of  two  black  discs,  that  represent  the  area  of  the  projectile  and  fragment  nuclei. 

Examples  of  projectile-only  fragmentation  are  provided  in  Figure  C-3.  By  a  convention 
used  by  Heckman  etal.  (1978),  these  events  are  referred  to  as  “Type  1.”  In  example  a, 
S6Fe  is  fragmented  into  52Cr  and  4He.  While  in  example  s  b  and  c,  only  one  nucleon  is 
stripped  off,  respectively  a  neutron  and  proton.  Reactions  of  the  type  described  in 
examples  b  and  c  present  difficulties  in  experiments  designed  to  evaluate  inelastic  cross- 
sections.  For  both  solid-state  detectors  and  photographic  nuclear  emulsions,  lack  of  charge 
in  neutrons  prevents  detection  of  neutron  loss-only  fragments.  While  in  the  case  of  proton 
loss-only  fragmentation,  the  ionization  can  be  too  low  to  be  detected  in  some  photographic 
nuclear  emulsion  experiments  (Heckman  etal.  1978). 

Target  Fragmentation 

1.  The  distribution  of  reaction  products  is  not  significantly  dependent  on  heavy-ion 
energy  between  0.3  and  2  GeV/nucleon. 

2.  Target  fragmentation  reaction  products  are  nearly  independent  of  projectile  mass 
except  for  interaction  cross-section  scaling  factors.  Examples  of  target-only  fragmentation 
are  provided  in  Figure  C-4.  Heckman  etal.  (1978)  referred  to  these  as  ‘Type  4.”  Example 
a  is  typical  in  that  it  contains  multiplicity,  with  three  fragments.  The  other  examples  are 
less  typical  but  possible,  and  contain  single  nucleon-loss  fragmentation  of  a  proton  and 
neutron,  that  normally  are  difficult  to  interpret  by  photographic  film  techniques. 
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Figure  C-3.  Projectile-Only  Fragmentation  (Type  1). 

Before  After 
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Figure  C-4.  Target-Only  Fragmentation  (Type  4) 


Before  After 
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Nuclear  absorption  interactions  of  heavy  projectiles  more  commonly  are  involved  in 
fragmentation  of  both  the  target  and  projectile.  Interactions  of  this  type  are  referred  to  as 
“Type  3.”  Figure  C-5  provides  examples.  The  characteristics  of  the  reaction  products, 
respectively,  are  similar  to  those  noted  for  projectile-  and  target-only  fragmentation. 

Figure  C-5.  Target  and  Projectile  Fragmentation  (Type  2). 

Before  After 


153 


Absent  from  the  two  cases  in  Figure  C-5  are  target  or  projectile  fragments  with  a  single 
nucleon-only  removal  fragment.  Single  nucleon-only  removal  is  more  common  with 
grazing  collisions  that  have  small  target-projectile  nuclear  overlap  whereas  collisions  of  this 
type  have  greater  nuclear  overlap  and  thus  higher  multiplicity  in  projectile  fragments. 

Central  Collisions 

Central  collisions  involve  head  on  collisions  resulting  in  both  target  and  projectile 
fragmentation  that  result  in  large  amounts  of  energy  transferred  to  the  target.  Some 
researchers  have  postulated  that  nuclear  shock  waves  are  produced  by  these  interactions 
(Curtis  and  Schimmerling  1977).  The  end  result  of  these  interactions  is  illustrated  in  Fig. 
C-6  where  there  is  a  catastrophic  destruction  of  both  target  and  projectile,  and  no  forward 
fragments  from  the  projectile. 

Figure  C-6.  Central  Collision  (Type  3). 

Before  After 


Model  Predictions  and  Experimental  Data 

Wilson  and  Costner  Theoretical  Data  on  Cross-Sections 

Wilson  and  Costner  (1975)  presented  a  high-energy  heavy-ion  reaction  model  that 
provided  both  total  and  absorption  nuclear  interaction  cross-sections.  The  model  was 
compared  to  experimental  data  with  close  agreement.  Figure  C-7  provides  cross-sections 
for  iron  projectiles  on  carbon  from  the  model.  A  recent  measurement  from  Dudkin  etal. 
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(1990)  of  iron  on  carbon  is  noted  on  the  graph  with  reasonable  agreement  with  the  model. 
The  ratio  of  absorption  cross-section  to  total  cross-section,  about  55  %,  is  reasonably 
independent  of  energy. 


Figure  C-7.  Total  and  Absorption  Nuclear  Interaction 
Cross-Sections  for  Fe  Projectiles  on  Carbon 


Dudkin  etal.  (1990  and  1994)  respectively  found  for  0.48  and  1.8  GeV/nucleon  Fe 
projectiles,  close  agreement  in  measured  and  theoretically  calculated  inelastic  cross- 
sections  using  the  Brant-Peters  formula: 

°abs  =  3i  Tq  [  Ap3  +  Aj3  -  B  ( Ap1/3  +  Aj13  )  ]2,  (C-2) 

where  Ap  and  Ax  are  the  respective  masses  of  projectile  and  target  nuclei,  r0  =  1.32  fm,  and 
B  =  0.85. 
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Westfall  Experimental  Data  on  Fragmentation 

Westfall  etal.  (1979)  measured  the  element  production  cross-section  of  56Fe  on  many 
targets  with  a  solid-state  detector  transmission  system.  The  results  of  the  carbon  target  are 

shown  in  Fig.  C-8.  From  the  graph,  it  is  obvious  that  small  number  nucleon  losses  (AZ:  1 

-  4)  account  for  over  half  of  the  projectile  fragmentation  events,  excluding  central  collisions 
that  could  not  be  evaluated.  Total  charge-changing  cross-section,  o&z>l,  and  mass¬ 
changing  cross-section,  cr^,  for  the  carbon  target  were  1.56  +  0.05  and  1.66  ±  0.06 

bams,  respectively.  Neutron-only  removal  cross-sections  could  not  be  measured  and  were 
estimated  at  100  millibams  (mb). 


Figure  C-8.  Element  Production  Cross-Sections  from  56Fe 
Projectiles  on  Carbon  at  1.88  GeV/nucleon  (Westfall  etal.  1979). 
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Dudkin  Experimental  Data  on  Fragmentation 

Dudkin  etal.  (1990, 1991, 1993  and  1994)  evaluated  target  and  projectile  fragment 
characteristics  from  1.8  GeV/nucleon  56Fe  incident  on  nuclei  of  photographic  nuclear 
emulsion.  They  separated  interactions  into  one  of  three  categories  based  on  the  target 
nuclei:  Ag  or  Br;  H;  and  C,  N,  or  O.  Only  a  summary  of  the  CNO  data  will  be  presented 
here  since  it  is  most  applicable  to  secondaries  created  in  A- 150  plastic. 

In  the  1990  work,  fragments  of  interactions  were  separated  into  one  of  five  types: 

1.  black (b)-Ep< 26 MeV, 

2.  gray  (g)  -  26  <  Ep  <  400  MeV, 

3.  stream  (s)  -  Ep  >  400  MeV, 

4.  stream  (s’)  -  Ep  >  400  MeV  (singly  charged  projectile  particles), 

5.  heavily  ionizing  (h)  particles  from  target  breakup,  and 

6.  fragments  of  the  projectile  nuclei. 

Total  charge  of  the  non-interacting  fragments  of  the  projectile  nucleus  Q  =  2.  was 

determined  for  each  interaction,  where  nt  is  the  number  of  fragments  and  Z;  is  the  charge  of 
all  relativistic  single-charged  fragments  within  5°  of  the  projectile  trajectory. 

They  determined  the  fragmentation  cross-section  for  different  projectile  products  as 
follows:  0.49 ±  0.06  bams  (20  < Z<  26),  0.50  ± 0.08  bams  (10  <  Z <  19), ),  0.38  +  0.06 

bams  (6  <  Z  <  9), ),  0.36  +  0.06  bams  (3  <  Z  <  5), ),  3.38  ±  0.31  bams  (a-particles), 

7.02  +  0.57  bams  (protons).  They  found  the  mean  emission  angle  of  the  fragments  in  each 
group  to  be:  1.03  ±  0.08°  (20  <  Z  <  26),  1.16  +  0.08°  (10  <  Z  <  19),  1.44  ±  0.13°  (6  <  Z 

<  9),  2.00  +  0. 15°  (3  <  Z  <  5),  2.50  +  0.08°  (a-particles),  and  2.78  +  0.04°  (protons). 

Thus,  for  increasing  fragment  size,  the  mean  angle  of  emission  decreases.  They 
determined  the  relative  percent  probabilities  of  definite  number  of  projectile  fragments  of 
S6Fe  from  interactions  with  CNO  target  nuclei.  Figure  C-9  is  adapted  from  the  data. 
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Graph  C-9.  Relative  Percent  Probabilities  of  Definite  Number 
of  56Fe  Fragments  [Adapted  From  Dudkin  etal  (1990)]. 


From  the  Figure,  about  2  %  of  the  interactions  did  not  have  any  relativistic  charged 
particles  emitted  along  the  initial  trajectory  of  the  projectile.  These  interactions  are 
categorized  as  central  collisions  (type  3).  The  complexity  of  the  projectile  fragmentation 
interactions  is  apparent  from  the  data.  For  example,  about  63  %  of  the  interactions  have 
one  projectile  fragment  with  Z  >  3  and  10.9  %  have  two  projectiles  with  Z  >  3.  Combined, 
these  interactions  are  predominately  grazing  collisions  that  leave  a  large  portion  of  the 
projectile  intact.  On  the  other  hand,  about  25  %  of  the  interactions  had  no  projectile 
fragments  of  Z  >  3  and  were  dominated  by  collisions  with  large  overlap  of  target  and 
projectile  nuclei.  One-third  of  the  interactions  produced  five  or  more  proton  fragments,  at 

least  70  %  produced  at  least  two  protons,  and  over  one-third  had  two  or  more  a-particles. 

Studies  were  conducted  on  CNO  target  fragments  produced  by  S6Fe  projectiles.  Per 
the  classification  system  given  above,  the  means  were  as  follows:  1.77  +  0. 10  (nj. 
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2.87  ±  0.14  (it ),  4.64  +  0.13  (r^),  8.12  +  0.50  (ns),  3.90  +  0.18  (« ,),  and 


20.05  +  0.35  (Q).  Figure  C-10  provides  correlations  of  the  mean  Q  “( Q  )”,  (r^),  (n h  ), 
and  (  ns )  with  respect  to  ng,  the  number  of  gray  particles. 


Figure  C-10.  Correlations  of  ( Q ),  (  ),  { ),  and  (ns)  with  ng 

from  S6Fe  on  CNO  [Adapted  from  Dudkin  etal.  (1990)]. 


From  the  figure,  ( Q )  and  ( nb )  are  negatively  correlated  to  ng,  while  { tih )  and  ( n3 )  are 
positively  correlated.  Independent  of  ng,  the  mean  number  of  black  particles,  <  >,  is  less 

than  the  shower  particles,  <  ns ),  or  the  heavily  ionizing  particles,  <  ).  Dudkin  etal. 

(1993)  later  determined  the  correlation  between  (  Q  ),  (  ns ),  and  (  ng  )  with  impact 


159 


parameter  for  various  projectile  and  target  nuclei  combinations.  Figure  C-l  1  provides 
values  of  <  Q )  and  ( )  and  vs.  impact  parameters  for  56Fe  projectiles  on  CNO. 


Figure  C-l  1.  ( Q )  and  (  ng )  vs.  Impact  Parameter  Range 
of  S6Fe  on  CNO  (Dudkin  etal.  1993). 


As  expected,  for  grazing  interactions  that  have  little  nuclei  overlap,  high  conservation  of 
mean  projectile  charge  and  a  decrease  in  the  number  of  gray  particles.  Another  point  of 
interest  from  the  mean  values  of  Q  and  ns,  is  that  the  sum  is  less  than  23.95,  less  than  the  Z 
of  iron.  Thus,  for  some  interactions,  conservation  of  charge  would  not  exist  along  the 
particle  trajectory. 


Dudkin  et  al.  (1991)  analyzed  characteristics  of  target  fragment  protons  and  a-particles. 

Table  C-l  lists  characteristics  S6Fe  projectiles  interacting  on  CNO  for  events  that  had 
almost  complete  target  breakup.  Particles  were  separated  into:  protons  of  energy  30  -  400 

MeV,  protons  <  30  MeV ,  and  a-particles  less  than  40  MeV.  For  each  category,  mean 
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values  of  multiplicity  (n),  emission  angle  (&),  energy  ( E),  transverse  momentum  ( P±  ), 

and  longitudinal  momentum  (  Py  )  or  forward/backward  ratio  of  the  secondaries  (n  /  are 
provided. 


Table  C-l.  Characteristics  of  Proton  and  a-Particle  Secondaries  from  Fragmentation 
of  1.8  GeV/nucleon  S6Fe  on  CNO  of  Photographic  Nuclear  Emulsion  (Dudkin  etal.  1991). 


Secondary 

Particle 

Range 

(MeV) 

<»> 

<0  > 
(degree) 

<E> 

(MeV) 

<Pi> 

(MeV/c) 

<P«> 

(MeV/c) 

n  In 

proton 

30-400 

2.46  +  0.13 

48+1 

198  +  5 

359+8 

394+  14 

r 

1 

proton 

<30 

1.66  +  0.10 

84  +  2 

9.1  +  0.4 

91  +  3 

1.3  +  0.2 

a -particle 

<40 

0.18  +  0.03 

79+6 

17.2  +  1.2 

264+  19 

— 

2.2  +  0.5 

Simple  observations  can  be  drawn  from  the  data.  First,  for  the  high-energy  protons:  30  - 
400  MeV,  the  transverse  and  longitudinal  momenta  are  similar,  and  the  particles  are 

directed  primarily  in  the  forward  direction.  For  the  low-energy  protons  and  a-particles, 

there  is  a  considerably  higher  fraction  of  the  secondaries  emitted  in  the  backward  direction 
with  about  43  and  24  %,  respectively. 

Dudkin  etal.  (1994)  analyzed  0.1  -  0.5  GeV/nucleon  S6Fe  interactions  in  photographic 
emulsion.  The  range  of  velocity  of  these  projectiles:  0.43  -  0.76c  was  significantly  less 
than  the  0.94c  velocity  of  the  1.8  GeV/nucleon  projectiles  of  their  other  work.  The 
following  mean  multiplicities  were  reported  for  CNO  target  nuclei:  2.43  ±  0. 14  (nj,  2.26 
±  0.17  (ng),  4.64  ±  0.13  (nj,  3.04  ±  0.24  (ns),  3.46  ±  0.21  (ns.),  22.18  ±  0.29  (Q), .  In 
comparison  to  the  1.8  GeV/amu  data: 

1.  the  sum  of  the  mean  Q  and  ns.  is  25.64, 1%  higher  than  the  1.8  GeV/amu  data, 

2.  there  was  a  decrease  in  the  total  inelastic  cross-section,  and 

3.  there  was  a  lower  degree  of  destruction  of  the  colliding  nuclei. 
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Antonchik  etal.  (1982)  also  evaluated  mean  multiplicities  of  target  fragments  from 
high-velocity  heavy  ions.  They  compared  1.2  GeV/nucleon  14N  to  3.6  GeV/nucleon  12C 
and  found  no  appreciable  change  in  multiplicities  of  any  category  of  target  fragment.  Thus, 
the  energy  dependence  of  target  fragment  multiplicities  is  apparent  only  in  the  low-energy 
region. 

Heckman  Experimental  Data  on  Fragmentation 

Using  photographic  nuclear  emulsion  detectors,  Heckman  etal.  (1978)  measured  the 
absorption  interaction  cross-sections  of  2.1  GeV/nucleon  4He,  12C,  14N,  and  160 
projectiles.  By  examination  of  the  films,  they  were  able  to  categorize  fragmentation  events 
by  type  as  used  here  already  in  this  appendix.  Figure  C-12  lists  the  percent  interaction  of 
each  type  of  the  total  for  each  projectile.  One  notable  feature  from  initial  inspection  of  the 
data  is  that  for  each  projectile,  about  50  %  of  the  interactions  involved  a  change  in  mass  of 
both  the  projectile  and  target.  Also,  for  heavy  ion  projectiles  involved  in  interactions  of  this 
type,  rarely  are  mass  changes  of  the  neutron-only  type.  Another  notable  feature  of  the  data 
is  the  relatively  constant  proportion  of  interaction  type,  independent  of  projectile.  One 
difficulty  in  interpretation  of  nuclear  films  is  proper  correction  of  the  data  for  neutron(s)- 
only  removals.  The  authors  estimated  that  about  5.3, 6.0, 5.4,  and  4.8  %  of  the  interaction 
events  had  neutron(s)-only  removal,  respectively  in  order  from  lower  to  higher  mass 

projectile.  Also,  for  the  three  heavy  projectiles,  AZ  =  1  removals  from  the  target  and 

projectiles  were  difficult  to  assess  because  the  ionization  produced  by  singly  charged  ions 

was  too  low.  Corrections  for  AZ  =  0, 1  were  calculated  and  determined  to  belong  to  type  1 

and  4  interactions.  The  maximum  correction  possible  for  each  type  is  represented  by  a 
small  line  bar  above  the  data. 
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Percent  of  Total  Interactions 


Figure  C-12.  Classification  of  Interaction  Types  of  2. 1  GeV/nucleon  4He,  12C, 
14N,  and  160  on  Photographic  Nuclear  Emulsion  (Heckman  etal  1978). 
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Appendix  D:  Position  Determinations 

General  Description 

Figure  D-l  illustrates  the  arrangement  of  the  position-sensitive  detectors  (PSD1  and 
PSD2)  with  respect  to  the  TEPC  in  a  three  coordinate  system.  Each  PSD  consists  of  two 
devices,  one  oriented  to  determine  relative  position  in  the  x-direction  and  the  other  in  the  y- 
direction. 


Figure  D-l.  Arrangement  of  Position-Sensitive  Detectors  in 
Relationship  to  TEPC  in  a  Three-Dimension  Frame 


PSD  1 


TEPC 


PSD  2 


z 


Each  device  in  the  PSD  has  an  active  silicon  thickness  of  about  1000  pirn  (1  mm).  The  two 
devices  were  placed  in  contact  and  assumed  to  have  the  same  z-position.  The  z-position  for 
each  PSD  pair  relative  to  the  center  of  the  TEPC  was  measured.  All  x-  and  y-position 
values  were  normalized  to  the  line  traversing  the  three  detectors,  with  x,  y,  and  z  =  0  for 
the  TEPC  center.  This  was  necessary  since  the  position-sensitive  devices  had  a  small  x  and 
y  offset  from  the  center  of  the  TEPC. 
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Position  Calculation  Equations 


FORTRAN  code  was  written  to  calculate  the  perpendicular  distance  of  the  heavy  ion 
trajectory  from  the  detector  center  and  is  referred  to  as  the  impact  parameter.  For  particles 
traversing  the  cavity  of  the  detector,  chord  length,  /,  of  the  trajectory  in  the  cavity  was 
calculated,  with  the  maximum  being  12.7  mm,  the  diameter. 

For  the  calculations,  particle  trajectories  were  assumed  to  be  on  a  straight  line  and 
defined  by  Equation  D-l: 


x-x!  _  y-yt  _  z-Z! 
x2-xi  y2-yi  z2  - zx  5 


(D-l) 


where  x15  x2,  yp  and  y2  are  the  normalized  position  values  determined  by  the  PSDs  and  z, 
and  Zj  are  the  measured  locations  of  the  PSDs  relative  to  the  center  of  the  TEPC.  The 
direction  cosines  are  defined  as  «,  v,  and  w: 


X2  -  Xi 

u  =  cosa  =  — - - ,  (D-2) 

d 

v  =  cos(3  =  — — — ,  and  (D-3) 

d 

Z2  -  Zi 

w  =  cosy  =  — - - ,  (D-4) 

d 

where  d  is  the  distance  between  the  two  points  and  a,  p,  and  y  are  the  angles  formed 

between  the  particle  trajectory  and  the  respective  axes  x,  y,  and  z.  The  following 
relationship  exists  among  the  direction  cosines: 

u2+v2+w2=  1.  (D-5) 

For  simplicity  of  presentation,  /,  m,  and  n  are  defined  as: 
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/=X2  —  Xj, 

(D-6) 

m  =  y2-yi,and 

(D-7) 

n  =Z2  -Zj. 

(D-8) 

Rearranging  terms,  the  following  relationships  exist: 


(D-9) 


(D-10) 


(D-ll) 


Similar  equations  can  be  formulated  for  a  line  from  the  detector  center  perpendicular  to  the 
particle  trajectory.  The  direction  cosines  of  this  line  are  defined  as  u\  v\  and  w’: 


U  =  cos  a!  = 


x-0 
b  ’ 


(D-12) 


v'  =  cosp' 


(D-13) 


w'  =  cosy'  = 


z  -  0 
b 


(D-14) 


where  b  is  the  impact  parameter  and  a’,  p\  and  y’  are  the  angles  formed  between  the  line 

and  the  respective  axes  x,  y,  and  z.  Like  the  other  set  of  direction  cosines,  the  sum  of  the 
direction  cosines-squared  equals  one.  Because  the  two  lines  are  perpendicular,  the 
following  relationship  exists  among  the  six  direction  cosines: 

utf  +vV+wW  =  0.  (D-15) 
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Then  by  substitution. 


x  /  y  m 
bd  +b  d  + 


Zffl  ^ 

—  —  =  0  and 
bd 


h.  +  my  +  nz  =  0 . 


(D-16) 

(D-17) 


Substitution  of  the  line  equation  (Eq.  D-l)  into  this  equation  for  each  of  the  three  axes,  die 
coordinates  of  the  intercept  points  xp  y„  and  Zj  can  be  determined: 


XI 


(m2  + 


n 


‘)  xt  -m/y1  -n/zj 


(D-18) 


(  r  +  n2)yi  -m/xi-nmzi 

yi - '-*±-2 - 1 - L,  and 

dz 

+  m2)z1-n/X|  -  n  m  y1 


(D-l  9) 


(D-20) 


where 

d2  =  Z2  +  m2  +  n2  ■  (D-21) 

It  follows  that  the  impact  parameter,  b,  and  the  chord  length,  /,  can  be  calculated: 


b  =  Jx?  +y2  +z2  311(1  (D-22) 

/  =  2t/ r2  -  b2  ,  (D-23) 

where  r  is  the  detector  radius  (6.35  mm).  To  provide  an  ability  to  determine  the  trajectory 
of  individual  particles  in  the  detector  plane,  the  angle  of  incidence,  0,  was  calculated  for 

each  trajectory.  Because,  the  particle  trajectories  had  a  limited  angle  of  acceptance  in  the  z- 
dimension,  0  was  calculated  from  the  x:  and  y,  intercept  values: 
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0  =  tan’1!  — 

Ui, 


(D-24) 


Sign  corrections  were  made  with  the  final  output  angle  conforming  to  the  convention 
illustrated  in  Fig.  D-2,  where  the  cross-section  of  the  detector  is  shown  in  the  x-y  plane. 

Figure  D-2.  Angle  Convention  Used  for  Data  Analysis. 


Position  Uncertainty  Estimations 

Figure  D-3  is  a  histogram  of  impact  parameter  for  events  with  energy  deposition  greater 
than  170  keV  (y  =  255  keV  pm1)  in  the  thin-walled  detector.  Over  85  %  of  the  events  are 
distributed  about  the  cavity/wall  interface  (impact  parameter  =  6.35  mm).  For  events  with 
impact  parameter  between  5.7  and  7.0  mm  (140  events),  the  mean  was  6.356  +  0.017  mm 
with  a  standard  deviation  of  the  data  about  the  mean  of  0.23  mm.  The  standard  deviation 
value  for  this  distribution  of  events  is  a  reasonable  estimate  of  the  uncertainty  in  impact 
parameter  at  the  cavity/wall  interface,  provided  the  impact  parameter  interval  that  gives  rise 
to  these  events  is  small  compared  to  the  uncertainties  in  impact  parameter. 
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Number  of  Events 


Figure  D-3.  Distribution  of  Impact  Parameter 
for  Fe  Events  with  High  Lineal  Energy  in  the 
Thin-Walled  Detector  [Simulated  Diameter  =  1  /im] 


Impact  Parameter  (mm) 
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At  the  cavity/wall  interface,  about  12,000  particles  per  mm  of  impact  parameter  interval 
were  observed  from  plots  of  detector  response  from  selection  on  impact  parameter.  For  the 
140  particles  (the  number  of  particles  between  5.7  and  7.0  mm  given  in  Fig.  D-3),  the 
impact  parameter  interval  would  be  0.0125  mm,  based  on  12,000  particles  per  mm.  This 
value  is  much  lower  than  the  estimated  uncertainty  in  impact  parameter,  0.23  mm.  Thus, 
0.23  mm  will  be  used  as  an  estimate  of  uncertainty  in  impact  parameter  at  the  cavity/wall 
interface. 

Many  factors  may  contribute  to  the  uncertainty  in  impact  parameter.  Most  important  is 
the  uncertainty  in  the  signals  from  the  PSD.  Other  contributions  may  be  from  distortions  in 
the  shape  of  the  detector  wall,  multiple  coulomb  scattering,  and  errors  in  measured 
distances  between  components  in  the  experimental  set-up.  Errors  in  measured  distances 
between  components  in  the  experimental  set-up  are  assumed  to  be  negligible.  The 
contribution  from  wall  distortions  is  difficult  to  assess,  but  no  distortions  were  observed  in 
x-ray  photographs  of  the  wall.  Calculations  were  made  to  assess  the  impact  of  multiple 
coulomb  scattering  on  position  uncertainties.  The  effect  had  a  small  contribution  to 
position  uncertainties  compared  to  that  estimated  for  the  PSD  devices  alone. 

Estimation  of  PSD  Uncertainty 

The  estimated  uncertainty  in  impact  parameter  at  the  cavity/wall  interface  is  based  on  the 
combined  uncertainty  of  all  four  PSD  devices.  An  analysis  was  performed  to  estimate  the 
uncertainty  in  the  individual  PSD  under  the  assumptions  that: 

1.  uncertainties  were  based  only  on  uncertainties  in  PSD  devices  and 

2.  the  uncertainty  among  the  devices  was  equal. 

A  computer  analysis  was  performed  on  an  Excel l^1  spreadsheet  with  Crystal  Ball.^2 
Monte  Carlo-based  software.  The  analysis  was  based  on  parallel  particle  trajectories  with 
true  impact  parameter  of  6.35  mm  randomly  located  around  the  detector. 

1  Microsoft  Corporation,  1993-1994. 

2  Decisioneering  Inc.,  1988-1996. 
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An  initial  guess  of  the  PSD  uncertainty  was  input  into  the  spreadsheet.  A  Monte  Carlo 
simulation  was  performed  with  an  output  distribution  of  impact  parameter.  The  standard 
deviation  of  the  calculated  distribution  of  impact  parameter  was  the  estimate  of  uncertainty 
in  impact  parameter.  The  input  value  of  PSD  uncertainty  was  modified  until  the  simulated 
distributions  of  impact  parameter  yielded  an  estimate  of  uncertainty  in  impact  parameter 
reasonably  close  to  0.23  mm.  A  PSD  uncertainty  of  0.325  mm  in 30,000  Monte  Carlo 
simulations  yielded  a  standard  deviation  in  the  distribution  of  impact  parameter  equal  to 
0.228  mm,  within  one  percent  of  the  target  value. 

Monte  Carlo  simulations  were  performed  for  other  impact  parameters  using  a  PSD 
uncertainty  of  0.325  mm.  Estimated  uncertainties  in  impact  parameter  ranged  from 
a  low  of  0.214  mm  (impact  parameter  =  0.5  mm)  to  a  high  of  0.231  mm  (impact  parameter 
=  4.25  mm)  for  30,000  trials  each.  Thus,  uncertainties  in  impact  parameter  are  relatively 
constant,  independent  of  impact  parameter. 

Impact  Parameter  Bias  Estimates 

Impact  parameters  were  calculated  from  the  combined  signals  from  the  four  PSD 
devices.  In  the  calculations,  position  signals  were  squared.  As  a  result,  uncertainties  in 
PSD  signals  will  introduce  a  bias  in  the  calculated  impact  parameter. 

Analyses  were  carried  out  to  determine  bias  in  calculation  of  impact  parameter.  The 
calculations  were  performed  with  the  same  spreadsheet  described  above.  PSD  uncertainty 
was  set  equal  to  0.325  mm,  10  different  true  values  of  impact  parameter  were  evaluated, 
and  30,000  simulations  were  performed.  Fig.  D-4  contains  the  results  expressed  in  terms 
of  percent  bias  vs.  impact  parameter.  For  all  impact  parameters  evaluated,  impact 
parameter  was  positively  biased  compared  to  true  impact  parameter.  Bias  was  highest  for 
small  impact  parameters  and  very  small  for  large  values.  A  power  function  fit  to  the  data  is 

shown  as  a  solid  line  in  the  Figure.  One  standard  deviation  (la)  errors  were  not  displayed 

because  they  were  on  the  order  of  the  size  of  the  data  symbols  or  less.  The  data  was 
replotted  with  the  abscissa  on  a  power  function  scale  to  weight  the  errors  with  respect  to  a 
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/^-random  irradiated  detector.  From  the  Fig.  D-5,  it  is  apparent  that  a  relatively  low  fraction 
of  particles  have  significant  bias. 

Impact  Parameter  Determinations:  PSD1  vs.  Both  Pairs 

In  analysis  of  the  data  where  charge-changing  fragmentation  events  occurred  in  the  wall 
of  the  detector,  impact  parameters  would  be  incorrect  if  calculated  with  both  pairs  of  PSDs. 
This  data  was  evaluated  with  impact  parameter  estimated  from  PSDl-only.  Impact 
parameter  estimates  from  PSDl-only  will  have  greater  uncertainty  than  those  determined 
with  both  set  of  PSDs  because  the  beam  did  not  contain  particles  that  were  perfectly  parallel 
in  trajectory.  This  analysis  estimated  the  uncertainty  in  impact  parameter  for  PSDl-only. 

The  analysis  of  uncertainty  was  based  in  part  on  experimental  data  and  Monte  Carlo 
simulations.  Uncertainty  in  impact  parameter  is  based  on  the  uncertainty  in  PSD  x15  PSD 
yP  and  that  from  the  uncertainty  in  particle  trajectory.  The  latter  effect  is  due  to  the  non¬ 
parallel  beam.  Earlier  in  this  appendix,  uncertainty  in  PSD  response  was  estimated  at 
0.325  mm.  Particle  trajectory  can  be  evaluated  from  the  distributions  of  (Xj  -  x,)  and 
(y2  -  yi)-  Figure  D-6  contains  the  distribution  of  (Xj  -  x,)  for  Fe  particles  through  the 
detector  with  the  copper  target  The  data  is  normally  distributed  with  a  mean  and  standard 
deviation  of  -3.47  and  1.84  mm,  respectively.  One  standard  deviation  in  the  trajectory 
angle  of  these  particles  in  the  x-dimension  is  0. 17°,  based  on  a  distance  between  the  PSDs 
of  622  mm  and  the  standard  deviation  of  1.84  mm  in  the  distribution  of  (Xj  -  xx). 

1.  True  Impact  Parameter  =  5.5  mm 

A  Monte  Carlo  simulation  was  performed  for  an  impact  parameter  of  5.5  mm  in  PSD1. 
The  standard  deviation  of  1.8  mm  was  used  for  (x2  -  xj  and  (y2  -  y,).  PSD  uncertainties 
of  0.325  mm  were  used.  The  output  of  the  Monte  Carlo  simulation  was  impact  parameter. 

For  40,000  trials,  the  mean  of  the  distribution  was  5.567  +  0.004,  with  a  CV  of 
14.8  %.  The  ratio  of  the  mean  to  the  true  impact  parameter  was  1.01.  The  kurtosis  and 
skewness  of  the  distribution  were  3.03  and  -0.002,  respectively,  indicating  normality. 
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Figure  D-6.  Distribution  of  (X2  -  XI)  for 
Fe  Particles  through  the  Detector  with 
Copper  Target  [Simulated  Diameter  =  1  /#m] 


x2  -  xl  (mm) 
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An  evaluation  of  experimental  data  was  performed  to  confirm  the  results  and 
assumption  used  in  the  Monte  Carlo  simulation.  For  Fe  events  through  the  detector,  impact 
parameter  was  calculated  with  both  PSD  pairs  and  PSDl-only.  Fig.  D-7  contains  the 
distribution  of  the  ratio  of  impact  parameter  from  PSDl-only  to  that  calculated  from  both 
PSDs.  The  distribution  was  created  for  Fe  events  through  the  detector,  the  copper  target  in 
place,  and  impact  parameters  between  5  and  6  mm  (as  determined  by  both  PSDs).  The 
data  follows  a  normal  distribution  and  has  a  mean  of  1.007  and  CV  of  14.6%. 

The  data  is  closely  approximated  by  the  Monte  Carlo  simulated  distribution.  At  5.5 
mm,  impact  parameter  based  on  PSDl-only  had  a  slight  positive  bias  compared  to  impact 

parameter  based  on  both  PSDs.  The  la  uncertainty  in  impact  parameter  at  5.5  mm  was 

estimated  at  0.84  mm  by  the  Monte  Carlo  calculations. 

2.  Uncertainty  in  Impact  Parameter 

Uncertainties  in  impact  parameter  were  estimated  by  Monte  Carlo  calculations  for  other 
impact  parameters.  Distributions  of  the  ratio  of  the  impact  parameter  from  PSDl-only  to 
that  from  both  pairs  were  created.  A  summary  of  the  analysis  of  the  experimental  data  and 
paired  estimates  of  uncertainties  by  Monte  Carlo  calculations  are  in  Table  C-l. 
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Figure  D-7.  Impact  Parameter  Ratio: 
(PSDl-Only)  to  (PSD1  &  PSD2)  for  Fe 
Particles  through  the  Detector  with 
Impact  Parameters  between  5  and  6  mm 


Impact  Parameter  (PSD1  &  PSD2) 
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Table  D-l.  Ratios  of  Impact  Parameter  from  (PSDl-Only)  to  that  with  Both  PSD  Pairs, 
and  Estimated  Uncertainties  in  Impact  Parameter  by  Monte  Carlo  Calculations 
[Experimental  Data  with  Fe  Particles  through  the  Detector  with  Cu  Target  in  Place] 


Impact 

Parameter 

Range  (mm) 

Experimental  Data 
fPSDU  . 

(PSD1  &  PSD2) 

Estimated 
Uncertainty 
Monte  Carlo 

(la)  (mm) 

Mean 

Standard  Deviation 

Event  Number 

2-3 

1.067  ±  0.003 

0.337  ±  0.003 

10,467 

0.81 

3-4 

1.054  ±  0.002 

0.239  ±  0.002 

15,321 

0.84 

4-5 

1.029  ±0.001 

0.180  ±0.001 

20,519 

0.82 

5-6 

1.007  ±  0.001 

0.146  +  0.001 

24,786 

0.84 

6-7 

0.995  +  0.001 

0.124  +  0.001 

26,275 

0.85 

7-8 

0.988  ±  0.001 

0.1 12  ±0.001 

26,429 

0.89 

8-9 

0.983  ±  0.001 

0.100  ±0.001 

24,977 

0.90 

177 


Appendix  E:  Calculations  for  Random  Variables 


General 

The  energy  deposition,  AE,  in  the  PSDs  was  dependent  on  many  factors.  As  noted  in 

Chapter  4,  these  factors  are  stochastic  in  nature  and  can  be  treated  as  random  variables,  Xr . 
The  “average”  or  expected  value,  E[X],  for  the  random  variable,  Xr,  is  one  measure  of 
central  tendency.  Another  property  of  random  variables  is  the  variance,  Var[Xf],  and 
provides  some  indication  of  how  the  distribution  is  spread.  The  variance  of  a  random 
variable  is  defined  as: 

Var[X]  =  E[{X-E[X]}2],  (E-l) 

where  X  is  the  random  variable  and  E[X]  is  the  expected  value  of  the  random  variable. 

More  commonly  the  expectation  value  is  given  the  symbol,  pi,  and  the  variance,  o2. 

The  square  root  of  the  variance,  a,  is  the  standard  deviation  and  is  more  commonly  used  to 

describe  dispersion.  The  coefficient  of  variation  (CV)  is  the  quotient:  . 

Quotient  of  Two  Random  Variables 

A  function  of  two  random  variables  is  their  quotient.  Let  X  and  Y  be  two  random 

variables  with  respective  expectation  values,  pix  and  piY,  and  variances,  a^.and  o2^  The 

covariance  between  the  two  random  variables  is  defined  as  Cov  [X,  Y]  or  simply  oLx  r  In 

general,  no  exact  solutions  exist  for  the  expectation  and  variance  of  the  quotient  of  two 
random  variables  (Mood  etal.  1974).  Approximate  solutions  given  by  Mood  etal.  (1974) 
are  given  in  Eqs.  E-2  and  -3. 
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1  rr2 

T°x.Y 
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Px  2 
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(E-2) 


/£x. 

Hx 


2ax,rN 

PxMt^ 


(E-3) 


Equation  E-4  was  produced  by  substitution,  rearrangement,  and  under  the  simplification: 
p,rx-i  «  from  Eq.  E-2. 


o\Y  [CV2(X)  +  CV2(y)  -  CV2(% 
PxM'F  2 


(E-4) 


Estimation  of  the  CV  for  the  Covariant  Term  between  AE  vt  and  AE  xt 


Distributions  of  AE  y,,  AE  x1}  and 


AEyj 


AE  Xj 


were  displayed  in  Figs.  4-2  and  -3. 


Using  the  statistics  of  the  distributions  tabulated  in  Table  E-l,  an  estimation  of  the 
covariance  between  AE  yt  and  AE  \  was  made  using  Eq.  E-4.  The  covariance  was 


estimated  at  3. 1  MeV2  (a  =  1.8  MeV).  Percent  CV’s  for  the  distributions  are  listed  in  Table 
E-l  and  the  estimated  CV  of  the  covariant  term  between  AE  yt  and  AE  x,  is  listed  as  well. 
For  AE  y,  and  AE  x,,  the  covariant  term  represents  22  and  16  %  of  the  total  variance. 


respectively. 
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The  mean  of  the  quotient  distribution, 


AEjql’ 

AExi 


is 


very  close  to  the  quotient  of  the 


mean  of  the  distribution  of  energy  in  each  PSD, 


PAE  yt  N 
M’AE  xi  ) 


Thus,  the  assumption  made  in 


forming  equation  E-4  was  appropriate  in  this  case. 


Table  E-l.  Summary  Statistics  and  Estimated  Percent  CV 
of  the  Covariant  Term  between  AE  yj  and  AE 


Distribution 

(MeV) 

Percent  Coefficient  of  Variation 

a 

Total 

Covariant  Term 

AEy, 

258 

3.75 

1.45 

0.68 

AEx! 

279 

4.40 

1.58 

0.63 

mM1 

UiJ 

0.925 

0.0179 

NA 

NA 
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Appendix  F:  Summary  of  Miscellaneous  Calculations 

Charged  Particle  Equilibrium 

Calculations  were  performed  for  the  detector  based  on  response  to  a  ^-random  incident 
beam  of  Fe  particles.  The  response  was  compared  to  that  of  an  equivalent  homogenous  site 
irradiated  to  an  equivalent  spatially  infinite  beam.  For  the  homogenous  case,  dose  is 
defined  by  Eq.  F-l: 

DOSE  =  LLL,  (F-l) 

P 


where  <j>  is  the  particle  fluence,  L  is  the  stopping  power,  and  p  is  the  density  of  the 
medium. 

Dose  to  the  cavity  of  the  detector  is  based  on  the  integral  of  Eq.  F-2: 

kfnax  A  •  £(b) 

DOSE  =  2  - — —  AA , 

bmin  m 


(F-2) 


where  as  shown  in  Fig.  F-l,  b^  is  the  center  of  the  detector,  bmM  is  the  outer  boundary  of 
the  wall  of  the  detector,  AA  is  the  incremental  area  as  shown  by  the  annulus,  and  m  is  the 

mass  of  die  cavity  fill  gas.  In  the  figure,  the  cavity  defines  the  detector  cross-section  area, 

£ 

Ad.  Based  on  dA  =  2k  •  b  •  Ab;  y  =  -=• ;  and  m  =  V  •  p,  where  V  is  the  cavity  volume; 


Eq.  F-2  can  be  redefined. 


DOSE  =  <|)  2 


1  •  y(b) 


b  .  p  •  v 

Dmin  r 


2ji  •  b  •  Ab. 


(F-3) 
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With  V 


it  r^3 

6D'Ad  = 


and  substituting  the  expression  for  dose  from 


Eq.  F-l,  the  following  equivalence  should  exist  if  CPE  has  been  achieved. 


<|>  *  L 
P 


<t> 

P  •  ad 


'max 


2  y0>) 


'min 


2ji  •  b  •  Ab. 


Simplifying,  Eq.  F-4  is  solved  for  L. 


7  max 


2  y(b)  •  2ji  •  b  •  Ab 


L  = 


LD 


For  the  calculations,  y(b)  was  substituted  for  y(b). 


(F-4) 


(F-5) 


Figure  F-l.  Illustration  of  Parameters  Used  to 
Determine  Absorbed  Dose  in  the  Detector 
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Frequency  Mean  and  Dose  Mean  Lineal  Eners 


Frequency  mean  lineal  energy,  yF  was  calculated  with  Eq.  F-6. 


2  y(b)  •  b  •  Ab 


UlOA 

2  b  •  Ab 


Dose  mean  lineal  energy,  yD ,  was  calculated  with  Eq.  F-7. 


(F-6) 


_  bmjn 


max  __  ~ 

2  y(b)2  •  b  •  Ab 


Mean  Quality  Factor 


Mean  quality  factor,  Q ,  was  determined  for  the  distributions  of  y  vs.  impact 
parameter,  b,  with  Eq.  F-8. 


2  Q(y)  •  y(b)  •  b  •  Ab 


max  _ 

2  y(b)  •  b  •  Ab 
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